三维非稳定渗流自由面边界积分项的精确数值计算

利用固定网格法分析三维非稳定渗流问题时,将要面对两项积分难题:以自由面及单元表面为边界的空间积分及以自由面为边界的曲面积分。针对常用的任意8结点6平面三维普通单元,提出采用坐标变换及等参变换技术求取空间积分项的精确数值解;至于曲面积分项,建议改用单元非饱和区部分表面作为积分边界,经过坐标变换及等参变换处理积分边界后,利用高斯数值积分可求出曲面积分项的精确数值解。通过一个普通单元及一项均质半无限边界堤坝的实例分析,表明此方法的精确性和稳定性良好。     

波尔兹曼数字岩芯致密砂岩渗透率研究

致密砂岩渗透率在油气勘探开发、地应力测量及水库地质灾害等领域具有重要意义,但鉴于目前实验技术局限性,温压耦合渗透率测量尚无法通过实验手段实现. 在已有研究工作基础上,首次建立了基于D3Q27数字岩芯模型,并计算了高温压耦合低渗砂岩渗透率问题. 首先,以鄂尔多斯盆地某油田延长组致密砂岩为例,利用X射线CT断层成像技术岩芯获取10μm, 5μm, 2μm分辨率致密砂岩内部结构数据,应用基于量子力学第一性原理的D3Q27格子波尔兹曼数字岩芯模型建立数值模型. 进而,利用自编3DLBM程序分别计算了不同分辨率渗透率随围压(0～200MPa)、孔隙压(0～65MPa)和温度(25℃～180℃)变化规律,通过与Inc AUTOLAB2 000C岩石测试分析系统实验结果对比,验证了程序的可靠性,得到低渗砂岩断层最佳分辨率;最后,在并行CPU-GPU平台上计算了高温高压耦合(0℃～400℃, 0～1.4 GPa)下致密砂岩渗透率值及其各向异性随温压变化规律,并讨论了致密砂岩中水在达到超临界状态后对致密岩石内部结构的影响.     

带有全铰节点的三维杆系结构静力分析的回传波射矩阵法研究

推导出含有全铰节点的三维杆系结构的回传波射矩阵表达式,完善了具有任意连接和约束的空间杆系结构静力分析的回传波矩阵法.基于节点平衡方程和协调方程,推导出表达杆件近端位移和远端位移关系的传递分配矩阵及载荷源向量,并通过由对偶坐标系下近端位移和远端位移的关系获得结构的总体相位矩阵,再引入转列矩阵,进而推导出结构的回传波射矩阵,在此基础上求解以杆端位移为基本未知量的线性方程组,最终得到精确确定所有杆件的杆端位移及杆端内力的矩阵列式.给出了空间杆系结构算例分析,与有限元结果比较,验证了回传波射矩阵法的计算精度.     

On a 3D moving load problem for an elastic half space

The paper deals with 3D dynamic response of an elastic half-space loaded by a point force moving at a constant speed along a straight line on the surface. The problem is formulated within the framework of the asymptotic hyperbolic–elliptic model developed earlier by two of the authors. The validity of the model is restricted to the range of speeds close to the Rayleigh wave speed. Steady-state near-field solutions are derived in terms of elementary functions. Transient analysis of surface motion illustrates peculiarities of the resonance associated with the Rayleigh wave.     

Meshing Requirements for Tidal Modeling in the Western North Atlantic

This paper presents an a posteriori approach to unstructured mesh generation via a localized truncation error analysis and applies it to the Western North Atlantic Tidal (WNAT) model domain. The WNAT model domain encompasses the Gulf of Mexico, the Caribbean Sea, and the North Atlantic Ocean east to the 60°W Meridian. Herein, we pay particular attention to the area surrounding the Bahamas.

A bathymetric data set with fine resolution is employed in seven separate linear, harmonic simulations of shallow water tidal flow for seven different tidal-forcing constituents. Each set of simulation results is used to perform a truncation error analysis of a linear, harmonic form of the depth-averaged momentum equations for each of the seven different tidal-forcing frequencies. Truncation error is then driven to a more uniform, domain-wide value by solving for local node spacing requirements. The process is built upon successful research aiming to produce unstructured grids for large-scale domains that can be used in the accurate and efficient modeling of shallow water flow. The methodology described herein can also be transferred to other modeling applications.     

一类常微分方程组的紧支撑样条小波方法探讨

应用三次紧支撑样条小波插值函数得到了求一类常微分方程组数值解的隐式公式,并求得到其局部截断误差为O(-h5).在此基础上给出1个显式校正求解公式,并讨论得到其局部截断误差为O(-h4).     

一维细胞神经网络的完全稳定性分析

运用Lyapunov函数方法,讨论了一维细胞神经网络模型的完全稳定性问题,给出了四组使模型具有完全稳定的充分条件.     

Law of the iterated logarithm for the periodogram

We consider the almost sure asymptotic behavior of the periodogram of stationary and ergodic sequences. Under mild conditions we establish that the limsup of the periodogram properly normalized identifies almost surely the spectral density function associated with the stationary process. Results for a specified frequency are also given. Our results also lead to the law of the iterated logarithm for the real and imaginary parts of the discrete Fourier transform. The proofs rely on martingale approximations combined with results from harmonic analysis and techniques from ergodic theory. Several applications to linear processes and their functionals, iterated random functions, mixing structures and Markov chains are also presented.     

Square-mean weighted pseudo almost automorphic solutions for non-autonomous stochastic evolution equations

In this paper, we first introduce the concepts and properties of the square-mean weighted pseudo almost automorphy and the square-mean bi-almost automorphy for a stochastic process. With these preliminary settings and by virtue of the theory of the semigroups of the operators, the Banach fixed point theorem and the stochastic analysis techniques, we investigate the well-posedness of the square-mean weighted pseudo almost automorphic solutions for a general class of non-autonomous stochastic evolution equations that satisfy either global or only local Lipschitz condition. Moreover, we estimate the boundedness of attractive domain for the case where the only local Lipschitz condition is taken into account. Finally, we provide two illustrative examples to show the practical usefulness of the analytical results that we establish in the paper.     

Strong solutions for two-dimensional nonlocal Cahn–Hilliard–Navier–Stokes systems

A well-known diffuse interface model for incompressible isothermal mixtures of two immiscible fluids consists of the Navier–Stokes system coupled with a convective Cahn–Hilliard equation. In some recent contributions the standard Cahn–Hilliard equation has been replaced by its nonlocal version. The corresponding system is physically more relevant and mathematically more challenging. Indeed, the only known results are essentially the existence of a global weak solution and the existence of a suitable notion of global attractor for the corresponding dynamical system defined without uniqueness. In fact, even in the two-dimensional case, uniqueness of weak solutions is still an open problem. Here we take a step forward in the case of regular potentials. First we prove the existence of a (unique) strong solution in two dimensions. Then we show that any weak solution regularizes in finite time uniformly with respect to bounded sets of initial data. This result allows us to deduce that the global attractor is the union of all the bounded complete trajectories which are strong solutions. We also demonstrate that each trajectory converges to a single equilibrium, provided that the potential is real analytic and the external forces vanish.