Simulating the hydraulic characteristics of the lower Yellow River by the finite‐volume technique

The finite‐volume technique is used to solve the two‐dimensional shallow‐water equations on unstructured mesh consisting of quadrilateral elements. In this paper the algorithm of the finite‐volume method is discussed in detail and particular attention is paid to accurately representing the complex irregular computational domain. The lower Yellow River reach from Huayuankou to Jiahetan is a typical meandering river. The generation of the computational mesh, which is used to simulate the flood, is affected by the distribution of water works in the river channel. The spatial information about the two Yellow River levee, the protecting dykes, and those roads that are obviously higher than the ground, need to be used to generate the computational mesh. As a result these dykes and roads locate the element interfaces of the computational mesh. In the model the finite‐volume method is used to solve the shallow‐wave equations, and the Osher scheme of the empirical function is used to calculate the flux through the interface between the neighbouring elements. The finite‐volume method has the advantage of using computational domain with complex geometry, and the Osher scheme is a method based on characteristic theory and is a monotone upwind numerical scheme with high resolution. The flood event with peak discharge of 15 300 m³/s, occurring in the period from 30 July to 10 August 1982, is simulated. The estimated result indicates that the simulation method is good for routing the flood in a region with complex geometry. Copyright © 2002 John Wiley & Sons, Ltd.     

Unsteady source flow in weakly heterogeneous porous media

Average nonuniform flows in heterogeneous formations are modeled with the aid of the nonlocal effective Darcy's law. The mean head for flow toward source of instantaneous discharge in a heterogeneous medium of given statistics represents the fundamental solution of the average flow equation and is called the Mean Green Function (MGF). The general representation of the MGF is obtained for weakly heterogeneous formations as a functional of the logconductivity correlation function. For Gaussian logconductivity correlation, the MGF is derived in terms of one quadrature in time t and it is analyzed for isotropic media of any dimensionality d and for 3D axisymmetric formations. The MGF is further applied to determining the mean head distribution for flow driven by a continuous source of constant discharge. The large time asymptotic of the mean head is analyzed in details.     

正交各向异性岩体孔径变形法的二维地应力分析方法改进

根据正交各项异性材料的弹性理论,推导出偏轴坐标系下的三维、二维本构方程;利用应力变分法对钻孔周围的应力分布作出修正,进一步得出孔壁径向位移的计算公式。结合孔径变形法测量地应力的特点,建立正交各向异性岩体地应力的分析理论,通过一典型工程与传统的各向同性理论进行了对比。     

海浪场非平衡信息论熵的研究初探

鉴于传统熵与时间无关 ,不能反映海 -气耦合机制下海浪场的基本属性———不可逆性 ,文章将非平衡信息论熵的概念引入海浪场的研究。依据海浪场相空间的Fokker Planck方程 ,首次推导出了非平衡信息论熵的熵平衡方程、熵密度、熵流密度、熵产生密度、广义力和广义力密度的具体表达式。     

温度效应对盐岩力学特性影响的试验研究

对经历不同温度后盐岩的力学性能进行了试验研究,分析了盐岩应力-应变曲线、峰值应力、峰值应变、弹性模量等的变化情况。分析结果表明,温度对盐岩的物理力学特性影响较大,随着温度的升高,盐岩的力学性能劣化明显,峰值应力和弹性模量降低,峰值应变增加;峰值应力可以表示成围压与温度的函数,并从盐岩弹性模量随温度的变化规律入手,导出了热损伤演化方程和一维TM耦合损伤本构方程。     

Axisymmetric interaction of a rigid disc with a transversely isotropic half‐space

A theoretical formulation is presented for the determination of the interaction of a vertically loaded disc embedded in a transversely isotropic half‐space. By means of a complete representation using a displacement potential, it is shown that the governing equations of motion for this class of problems can be uncoupled into a fourth‐order partial differential equation. With the aid of Hankel transforms, a relaxed treatment of the mixed‐boundary value problem is formulated as dual integral equations, which, in turn, are reduced to a Fredholm equation of the second kind. In addition to furnishing a unified view of existing solutions for zero and infinite embedments, the present treatment reveals a severe boundary‐layer phenomenon, which is apt to be of interest to this class of problems in general. The present solutions are analytically in exact agreement with the existing solutions for a half‐space with isotropic material properties. To confirm the accuracy of the numerical evaluation of the integrals involved, numerical results are included for cases of different degrees of the material anisotropy and compared with existing solutions. Further numerical examples are also presented to elucidate the influence of the degree of the material anisotropy on the response. Copyright © 2009 John Wiley & Sons, Ltd.     

提高全球导航卫星定位系统精度和稳定性的算法研究

针对地面用户接收机的性能易受各种复杂环境影响的问题,研究了在遮挡、多径等复杂环境下,提高和保持定位精度及稳定性的定位算法。基于RAIM思想,对量测信息进行故障诊断和处理,给出了一种高性能、低复杂度的EKF算法。仿真和试验表明,该算法不但提高了定位的精度,而且具有较强的鲁棒性和较强的工程实用性。     

Estimating evapotranspiration fraction by modeling two-dimensional space of NDVI/albedo and day-night land surface temperature difference: A comparative study

Three simple models estimating the evapotranspiration fraction (EF) were tested, validated and compared with available data from NDVI/albedo and day-night land surface temperature difference (DT). They were analyzed using ground based measurements collected by the energy balance Bowen ratio system at the 11 enhanced facilities located at the Southern Great Planes of the USA between April 2001 and May 2005. One model was superior in predicting the EF. It used an extension of the Priestly-Taylor equation and a relation between the normalized difference vegetation index (NDVI) and DT.     

地震波数值模拟中差分阶数的边界效应

提出了声波正演数值模拟中计算网格间差分阶数(精度)的不衔接而引起的边界反射效应问题,采用不同中心网格有限差分法求解声波波动方程来验证.数值实例分析表明,同差分阶数间不存在任何边界效应,而当差分阶数较低且网格间差分阶数递变较大时,边界效应显著,通过缩小差分网格间的递变阶数并提高相应的离散阶数,可以有效压制该边界效应,并保...