The Role of Capillarity in Two-Phase Flow through Porous Media

When determining experimentally relative permeability and capillary pressure as a function of saturation, a self-consistent system of macroscopic equations, that includes Leverett's equation for capillary pressure, is required. In this technical note, such a system of equations, together with the conditions under which the equations apply, is formulated. With the aid of this system of equations, it is shown that, at the inlet boundary of a vertically oriented porous medium, static conditions pertain, and that potentials, because of the definition of potential, are equal in magnitude to pressures. Consequently, Leverett's equation is valid at the inlet boundary of the porous medium, provided cocurrent flow, or gravity-driven, countercurrent flow is taking place, and provided the porous medium is homogeneous. Moreover, it is demonstrated that Leverett's equation is valid for flow along the length of a vertically oriented porous medium, provided cocurrent flow, or gravity-driven, countercurrent flow is taking place, and provided the porous medium is homogeneous and there are no hydrodynamic effects. However, Leverett's equation is invalid for horizontal, steady-state, forced, countercurrent flow. When such flow is taking place, it is the sum of the pressures, and not the difference in pressures, which is related to capillary pressure.     

基于CFD-DEM的高炉风口回旋区形状和传热特性数值模拟研究

采用计算流体力学-离散单元法(CFD-DEM)耦合方法,对高炉风口回旋区进行了数值模拟研究。首先通过与实验结果对比,验证了CFD-DEM模型的正确性;然后考察了不同气速对风口回旋区形状和传热特性及颗粒接触的影响。数值模拟结果表明:风口回旋区的大小和形状均受气速影响较大,在较大进气速度下,颗粒受到的曳力大于颗粒间的摩擦阻力并破坏颗粒间的桥力,形成较大尺寸的回旋区;且颗粒之间接触力较小,形成较大的空隙结构,更有利于热气体向周围扩散以强化传热。目前考察的三种气速结果表明:当气速为11m/s时,热量向下方传递速度最快;当气速为13m/s时,热量向上方传递速度最快;而当气速为15m/s时,热量向右方传递速度最快;此外,气速越大流态化越明显,颗粒间接触越少,接触力也越小。     

Hydraulic analysis of free-surface flows with solidification

Summary The paper is concerned with a one-dimensional analysis of plane open-channel flow with continuous solidification. The process is of relevance for recent developments in the casting of steel and other metals. The bottom of the channel consists of a rotating casting roll and a horizontal cooling table, where the solidified material is withdrawn with given velocity. The study is restricted to the region downstream of the top of the casting roll. Surface tension is neglected. In the main part of the analysis inviscid fluid flow is considered since the Reynolds number is very large in the applications. It is found that the steady-state solutions are nonunique in a certain parameter range. In addition to a continuous solution, there are two solutions including hydraulic jumps, with one hydraulic jump being located on the casting roll, the other one on the cooling table. Regarding the stability of the non unique solutions, the evolution of disturbances is investigated numerically as an initial-value problem. It is concluded that the hydraulic jump on the cooling table is unstable, while the other discontinuous solution as well as the continuous solution are stable for sufficiently small disturbances. Which stable solution is attained in the steady state, depends on the history of the process. Friction at the liquid/solid interface is taken into account in the last part of the analysis. A constant friction coefficient is assumed. It is found that the history of the process determines the steady-state solution if, and only if, the friction coefficient is sufficiently small. For larger values of the friction coefficient, the steady-state solution is unique and independent of the history of the transient process. Furthermore, for sufficiently large friction coefficients, stable hydraulic jumps are found, in contrast to the inviscid case, also on the cooling table. Received 19 March 1999; accepted for publication 3 May 1999     

Non-linear argumental oscillators: Stability criterion and approximate implicit analytic solution

The behaviour of a space-modulated, so-called “argumental” oscillator, is studied. The oscillator is submitted to an external harmonic force, which is amplitude-modulated by the oscillator's position in space. An analytic expression of a stability criterion is given. Using the averaging method, an integrating factor and a Van der Pol representation in the (amplitude, phase)-space, an exact implicit analytic solution is given when there is no damping, and an approximate implicit analytic solution is given when there is damping, allowing the plotting of the separatrix curve. An attractor is identified.     

小波分析用于陀螺仪漂移测试信号分析的研究

为建立陀螺仪漂移特性模型，需对其进行测试。由测试得到的数据是含有噪声的，且一般是非平稳的，用小波分析测试数据是一种很有效的方法。在简述小波分析中的多分辨分析理论的基础上，用其对某型陀螺仪的实测数据进行了分析－预处理，其结果证实了所研究方法的有效性。     

A VIEWPOINT OF THE TENSORAL PROPERTY OF COORDINATE TRANSFORMATION COEFFICIENT AND ANALYSISON THE CONCEPT OF HYBRID TENSOR1)

本文致力于澄清一个十分基本的问题：坐标变换系数是否为张量?传统观念认为,坐标变换系数不是张量。为了揭示坐标变换系数的本质,本文采用“从一般到特殊”的研究策略,重塑了张量的内涵和外延,引入了杂交张量概念,进而颠覆了坐标变换系数不是张量的传统观念,确切地讲,它就是度量张量的杂交分量。这一结果扩张了张量概念的集合,提升了张量分析学内在的统一性、对称性和不变性,减少了连续介质力学的运算量。     

Investigation of calculation method for anti-penetration of reactive power steel fiber concrete (RPC)

为了解决评估超高强活性粉末钢纤维混凝土的抗侵彻能力问题，根据混凝土在弹体侵彻过程中的介质状态及能量的分配关系，揭示了侵彻近区耗散能量结构与破碎区之间的比例尺度关系。由于近区介质的运动特性变形接近一维应变状态，受力类似于流体动力学状态，据此推导并简化了具有宽广的适用范围的侵彻计算公式，认为侵彻的深度主要由弹体的质量、形状、速度及靶体的强度、密度、断裂韧度、变形波速决定。通过多次的实弹试验，验证了公式的可靠性，得到了RPC的抗侵彻能力是普通混凝土的3倍的结论。     

浮筏的主动吸振控制实验研究

浮筏在低频减振方面存在较大不足,而这恰是主动减振的优势.将两者相结合,研究了浮筏的主动吸振控制问题.研制了一种主动吸振器,建立了以DSP F2812为核心的控制系统;提出了符号LMS算法,与广泛使用的滤波x-LMS等算法相比,该算法不需要精确的模型信息;在此基础上设计了自适应前馈主动吸振控制律,并在浮筏实验平台上进行主动吸振控制实验.实验结果表明,通过主动吸振控制,浮筏的低频减振能力得到了显著的改善.     

介质参数反演的广义射线近似方法

在对无粘性介质参数反演问题进行的研究中，引入一种全波场广义射线近似形式，提出一种新的反演参数的方法，文中，首先对由弹性波动方程演变成的声波方程进行分析，引入背景场量和扰动量，并结合Ｇｒｅｅｎ函数理论，得到了介质参数的积分方程；然后结合前人对非均匀介质中波函数局部理论的定性分析，引入一种全波场广度射线近似形式，把问题归结为一个第一类Ｆｒｅｄｈｏｌｍ积分方程；最后对半空间问题层状介质模型进行了反演，算     

Steady sliding frictional contact problem for a 2d elastic half-space with a discontinuous friction coefficient and related stress singularities

The steady sliding frictional contact problem between a moving rigid indentor of arbitrary shape and an isotropic homogeneous elastic half-space in plane strain is extensively analysed. The case where the friction coefficient is a step function (with respect to the space variable), that is, where there are jumps in the friction coefficient, is considered. The problem is put under the form of a variational inequality which is proved to always have a solution which, in addition, is unique in some cases. The solutions exhibit different kinds of universal singularities that are explicitly given. In particular, it is shown that the nature of the universal stress singularity at a jump of the friction coefficient is different depending on the sign of the jump.     

自适应卡尔曼滤波在惯性测量组合误差补偿中的应用

惯性元件误差是捷联惯导系统的主要误差源，必须在导航过程中加以补偿。根据机动目标跟踪理论和惯性测量组合动态模型，分别建立状态方程和观测方程，利用机动频率自适应的算法进行卡尔曼滤波，以此达到惯性测量组合动态误差和随机误差补偿的目的。仿真结果说明该方法可行有效，优于传统的误差补偿算法，能较好地提高系统导航精度。     

基于伪距宽巷组合的GPS模糊度固定方法

中长基线模糊度快速解算是GPS网络差分技术的核心内容。根据网络差分模糊度固定只用于基线解算而不用于定位的特点,通过对常用GPS原始观测数据组合方式的分析,提出伪距宽相组合的数据处理方法,有效消除了电离层、对流层传播误差,形成抗差性强、大气误差自由的GPS组合观测值;在此基础上引入双差伪距宽相组合进行模糊度浮点解并建立法方程,应用高度角与大气误差的关系构造权阵,对常规LAMBDA算法进行了改进,形成一种适合GPS网络差分的中长基线模糊度解算方法。对三个参考站同步观测数据的实际测试结果表明：使用该方法网络模糊度解算时间小于300s,基线长超过60km,并满足闭合性原则。     

伪可变体系的几何构造分析

伪可变体系几何可变性的研究，对轻型结构的设计分析已变得十分重要。本文先分析能量与平衡之间的普遍关系，进而得出判定体系可变性的能量准则。通过拉格朗日乘子的引入，建立能量泛函，得出判定极值的二次型。然后证明了乘积力法与能量法的一致性，并讨论了宜于计算机分析实现的矩阵表示方法。结果表明，若二次型确定，则体系伪可变；当半确定时，体系部分伪可变部分可变；否则体系含二阶以上的无穷小机构。     

Entropy of microstructure

Two points are made in this paper: first, energy of random structures is not determined uniquely by any finite set of the characteristics of microstructure. The information lost is characterized by entropy of microstructure; it describes the scattering of the values of energy. Therefore, entropy of microstructure is a key thermodynamic parameter in phenomenological modeling of the behavior of random structures. Second, mathematical modeling of a random structure is based on the construction of its probabilistic measure; a way to select the probabilistic measure from the experimental data is outlined. The corresponding probabilistic measure is remarkably similar to that of classical statistical mechanics, though the underlying physics is quite different. After the probabilistic measure is chosen, the entropy of microstructure can be found from the analysis of the homogenization problem. Entropy of microstructure is computed in two example problems. Applications to phenomenological modeling of work hardening are discussed.     

基于最小搜索超椭球的GNSS模糊度固定及检验方法

为进一步简化GNSS模糊度解算流程,降低计算复杂度,重点提高LAMBDA算法的搜索效率,对模糊度解算作出以下改进:1)模糊度检验采用后验概率检验方法,并对其目标函数进行适当简化;2)把简化后的目标函数嵌入模糊度搜索过程,省去了单独的模糊度检验步骤;3)推导了模糊度空间最小搜索超椭球,把搜索区域限制在该超椭球中,缩小了搜索范围,从而大大降低了搜索复杂度。用三组实测数据实验比较了新方法和传统的LAMBDA方法,结果显示新方法搜索复杂度降低普遍在30%左右,最高可接近60%。理论推导和实验结果均证明了新方法的高效性。     

某软土深基坑工程时间效应有限元分析

软土地区的深基坑因土体的固结作用和流变性而具有了时间效应。本文以Biot固结有限元法为基础,用三元件模型中的第一个线性弹簧模拟固结作用,弹性模量考虑了开挖应力路径和应力历史的影响;另外一部分(KELVIN模型)来模拟土体的流变性,以实际变形的反演来得到两个参数的大致取值,再对基坑的变形情况以及进一步开挖进行分析。假定为正常固结饱和粘土,平面应变问题。通过对某饱和软粘土地基深基坑开挖工程实例的分析,得到的挡墙水平位移曲线与实测曲线很吻合,表明程序较好地反映出基坑的时间效应。     

单步辛算法的相位误差分析及修正

若一个算法的幅值误差和相位误差都不累加，则该算法就是最理想的算法， 但这样的算法难以构造. 辛几何算法解决了幅值误差的累加问题，但相位误差累加问题仍然 存在. 给出了单步隐式辛算法相位误差的精确估计公式，提出了简单而实用的修正方法. 以 Euler中点隐式辛差分格式为例，针对几个线性动力学系统，对相位误差进行了数值分析和 修正.     

GLONASS星历电文特征及其解算方法

针对飞行试验中由于GLONASS星历解算错误导致的定位结果异常问题,研究了GLONASS星历电文下传的基本特征,即通过第1~4串电文下传,并且一个更新周期内下传60组数据完全相同的星历。基于此,提出了基于电文串标识的GLONASS星历解算基本算法,但发现该算法在电文串丢失且发生星历更新时解算出错误星历的问题。为确保星历来源于连续的1~4串电文,提出了基于时间比较的星历解算改进算法,发现GLONASS星历更新时,不保证从第1串开始,也不保证在连续的1~4串电文中更新完毕,改进算法依然无法确保获取到正确星历。最后分析了星历电文误码时的特征,提出了基于星历合法性检测的星历解算可靠算法,该算法综合考虑卫星不健康、星历更新以及电文误码等异常情况,采用轨道特性检测法和原码比对检测法验证星历合法性,采用电文串标识法充分利用有效电文数据。试验结果表明,该算法的星历误码识别率达到100%,星历更新异常识别率达到100%,获取的GLONASS星历数据正确率到100%。     

Variational formulation of the equivalent eigenstrain method with an application to a problem with radial eigenstrains

In this paper a variational formulation of the equivalent eigenstrain method is established. A functional of the Hashin–Shtrikman type is proposed such that the solution of the equivalent eigenstrain equation is a unique minimizer of the functional. Moreover, it is also shown that the equivalent eigenstrain equation is the Euler–Lagrange equation of the potential energy of the inclusions. An approximate solution of the equivalent eigenstrain equation is then found as a minimizer of the functional on a finite dimensional span of basic eigenstrains. Special attention is paid to possible symmetries of the problem. The variational formulation is illustrated by determination of effective linear elastic properties. In particular, material with a simple cubic microstructure is considered in detail. A solution for the polynomial radial basic eigenstrains approximation is found. In particular, for the homogeneous eigenstrain approximation, the effective moduli are derived in an exact closed form.