二维振荡叶栅非定常粘性流动数值模拟

采用显式四步Runge-Kutta格式，结合Baldwon-Lomax紊流模型求解Navier-Stokes方程，借助运动网格技术，完成了对二维振荡叶栅非定常粘性流动的数值模拟。为了加速求解过程，引入了变系数隐式残差光顺方法，取得了较好效果。数值结果与已公布的数据有很好的一致性。     

高黏性流动有限元模拟的迭代稳定分步算法

分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于 某个临界值方能使算法稳定. 然而在高黏性流动模拟中，已有的显式和半隐式分步算法由于 其显式本质，必须采用小时间步长计算，不但降低了计算效率，同时也常与为使分步算法稳 分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于 某个临界值方能使算法稳定. 然而在高黏性流动模拟中，已有的显式和半隐式分步算法由于 其显式本质，必须采用小时间步长计算，不但降低了计算效率，同时也常与为使分步算法稳 定必须满足的最小时间步长要求冲突. 本文目的是构造一种含迭代格式的分步算法，它能在 保证精度的前提下大幅度地增大时间步长. 方腔流和平面Poisseuille流数值计算结果证实 了此特点，该方法被有效应用于充填流动过程的数值模拟.     

射流元件控制通道流动的三维数值模拟

针对一种射流元件控制通道的复杂结构 ,采用分块对接技术和网格“融合”技术生成计算网格 ,并运用五步显式格式的 Runge-Kutta法和多重网格法对含全 N-S方程、RNG k-ε湍流模型和两层分区壁面模型的流动模型进行数值求解。通过对控制通道内部流动的数值模拟和流场特性分析 ,提出了改进方案     

RTM充模过程数值模拟的隐式有限元算法

建立了基于欧拉方法描述树脂传递模塑(RTM)工艺充模过程的基本数学方程,并采用有限元隐式时间积分方法对基本方程进行了数值求解．编制了基于隐式有限元算法及传统有限元控制体算法的程序,通过具体算例比较了这两种算法的优缺点．与传统的有限元控制体法相比,该文提出的隐式有限元算法能节省计算时间,特别适合于单元、节点数目多的情况．隐式有限元算法是一种纯有限元方法,不需要使用控制体积技术,采用该算法计算出的流动前沿与时间步长无关。     

绕Apollo飞船的高超声速化学非平衡流动的数值模拟

利用混合通量分裂方法，建立了很方便求解的隐式ＮＮＤ格式，求解了完全气体和化学非平衡空气绕Ａｐｏｌｌｏ飞船的流动，计算结果和实验值作了比较，应用拓扑分析方法，研究了背风区和尾迹内的流动结构。     

一种基于块雅可比迭代的高阶FR格式隐式方法

最近, 基于非结构网格的高阶通量重构格式(flux reconstruction, FR)因其构造简单且通用性强而受到越来越多人的关注. 但将FR格式应用于大规模复杂流动的模拟时仍面临计算开销大、求解时间长等问题. 因此, 亟需发展与之相适应的高效隐式求解方法和并行计算技术. 本文提出了一种基于块Jacobi迭代的高阶FR格式求解定常二维欧拉方程的单GPU隐式时间推进方法. 由于直接求解FR格式空间和隐式时间离散后的全局线性方程组效率低下并且内存占用很大. 而通过块雅可比迭代的方式, 能够改变全局线性方程组左端矩阵的特征, 克服影响求解并行性的相邻单元依赖问题, 使得只需要存储和计算对角块矩阵. 最终将求解全局线性方程组转化为求解一系列局部单元线性方程组, 进而又可利用LU分解法在GPU上并行求解这些小型局部线性方程组. 通过二维无黏Bump流动和NACA0012无黏绕流两个数值实验表明, 该隐式方法计算收敛所用的迭代步数和计算时间均远小于使用多重网格加速的显式Runge-Kutta格式, 且在计算效率方面至少有一个量级的提升.      

高雷诺数下求解NS方程的无网格算法

提出了一种适合高雷诺数NS方程求解的隐式无网格算法。针对高雷诺数粘性流动的特点,在附面层内的粘性影响区域采用法向层次推进布点的方法形成离散点云,在附面层外的计算区域内实行填充式布点的方法形成离散点云。根据附面层内外点云的不同构造特点,推导出运用格林公式和最小二乘曲面拟合方法求取空间导数的统一形式,在此基础上运用AUSM _up格式求得数值通量,并引入BL湍流模型对雷诺平均NS方程的湍流应力项进行封闭。时间推进格式方面,采用了计算效率较高的隐式高斯-赛德尔迭代算法。为了验证本文方法的计算精度和鲁棒性,对NACA0012翼型低速流动、RAE2822翼型跨音速绕流和二维圆柱的分离流动进行了数值模拟。     

非稳态流动的隐式最小二乘等几何方法

对非稳态粘性不可压流动问题提出了一种隐式最小二乘等几何计算方法。该方法先用隐式的向后多步差分格式对Navier-Stokes方程进行时间离散,再用Newton法线性化对流项,最后在每个时间步上用最小二乘等几何方法进行求解。根据该算法编制了计算程序,通过构造解析解的方法验证了程序的正确性,用该程序求解了雷诺数为5000时的非稳态二维顶盖驱动流问题,计算结果捕捉到了流动过程中涡的演化过程,表明本文方法可用于非稳态流动的求解。     

基于隐式扩散的直接力格式浸没边界格子Boltzmann方法

采用浸没边界格子Boltzmann(immersed boundary-lattice Boltzmann,IB-LB)模型执行动边界绕流数值模拟时,信息交互界面和边界力计算格式直接影响流动求解器的数值精度和计算效率.基于隐式扩散界面,一种改进的直接力格式IB-LB模型被提出.边界力表达式基于欧拉/拉格朗日变量同一性准...     

侧向多喷口干扰复杂流动数值模拟研究

采用具有高分辨率的NND格式,通过数值求解N-S方程对典型外形多喷口侧向喷流复杂干扰流动进行了数值模拟. 为了提高计算效率,采用了LU-SGS隐式算法. 采用分块对接网格技术,生成高质量的贴体计算网格,精确模拟喷口截面. 对比分析了不同计算格式、限制器形式、网格拓扑及流动形态(层流与湍流)对喷流干扰流场结构和压力分布特性的影响,研究和分析了喷口附近流场的涡系结构、波系结构和喷流干扰引起的气动力特性. 在上述研究的基础上,针对典型飞行器外形的侧向喷流干扰特性进行了详细的数值模拟,得到了喷口参数(喷口位置、数目等)及来流条件对喷流干扰流场结构、气动力特性的影响规律,并对其流动机理进行了相应的分析. 研究表明,发展的针对多喷口侧喷干扰的数值计算方法是成功的,可以应用于飞行器侧向喷流干扰的流场结构分析及气动力特性数值预测.      

基于Surfer平台的FLAC3D 复杂三维地质建模研究

复杂三维地质体模型的建立一直是FLAC^3D软件前处理过程中的一个难题。因此,提出以Surfer软件为过渡平台,通过对地表及岩层分界面三维地质信息的提取、转换,并将Surfer输出的数据经由FLAC^3D软件内嵌的Fish语言编制的程序的二次转换,生成FLAC^3D软件可以直接读取的模型数据文件,从而实现了该软件前处理过程中复杂三维地质模型的快速、准确建立。同时结合工程实际,检验了该方法的可行性和仿真效果,结果表明,建模过程简捷、实用,可操作性强,是FLAC^3D软件前处理的一种新方法。     

抗滑桩设计推力计算方法研究

滑坡推力是抗滑桩设计推力的重要依据,工程中通常使用剩余推力法来计算滑坡推力。为了比较准确合理地计算出抗滑桩的设计推力,提出了在剩余推力法基础上的一种新的抗滑桩设计推力取值方法。分析极限平衡状态下和工程安全要求下各条块的剩余下滑力,选择合理的桩位,使得抗滑桩位于滑坡体的抗滑段;然后调整桩后滑体的剩余下滑力,其调整量一方面保证桩前桩后滑体满足工程安全要求,另一方面则作为抗滑桩设计推力的修正依据,其过程由程序实现。通过实例证明了该方法的合理性。同时对今后抗滑桩设计推力的合理取值提供了参考作用。     

基于混沌神经网络理论的城市地面沉降量预测模型

通过分析城市地面沉降量时间序列的非线性动力学系统,认为该时间序列具有混沌特性。在此基础上,通过相空间重构的方法建立了用于城市地面沉降量预测的混沌神经网络模型;并利用此模型对高桥地面沉降量进行了预测,并和实际监测沉降量进行了比较,最大绝对预测误差为1. 7,预测的平均误差为0. 0833,研究结果表明,应用混沌神经网络模型进行城市沉降预测是可行、精确的。     

Analysis of the USBM wettability test

In this paper, we analyse the capillary pressure curves obtained by the centrifuge method in order to perform the USBM wettability test. The physical displacement mechanisms present both in the porous plate and in the centrifuge method, are described for different cases of wettability of the pore surface.The wetting fluid is defined as the fluid being at the lower pressure while displacing the other fluid, and this displacement is defined as imbibition. On the other hand, the process in which the fluid under the lower pressure is the displaced fluid is defined as drainage. The capillary pressure is defined as the positive pressure difference between the two fluids. By adhering to these definitions, there is a unique and consistent terminology for the same physical process: the displacement of oil by water in an oil wet system and the displacement of water by oil in an water wet system are both designated as drainage.An important result is that the centrifuge method is limited to the determination of drainage capillary pressure curves for strongly oil or water wet samples. There is no capillary equilibrium possible when a water wet sample is centrifuged under water because the wetting phase is under higher pressure than the nonwetting phase; the resulting forced imbibition curve should not be called a capillary pressure curve. For samples with bicontinuous fractional wettability, the curves obtained by the centrifuge method correspond to combination displacement, i.e. a combination of equilibrium drainage and forced imbibition coupled with blob mobilisation.     

冲击渐进振动系统相邻基本振动的转迁规律

冲击振动现象广泛存在于动力机械系统中,使得系统表现出复杂的动力学响应.目前对冲击振动系统的p/1类基本振动的稳定性及分岔研究报道较少,而且已有的对冲击振动系统动力学的研究基本都是基于单参数分岔进行分析的.研究以小型振动冲击式打桩机为工程背景,建立了冲击渐进振动系统的力学模型.分析了激振器和缓冲垫发生碰撞的类型,以及滑块渐进运动的条件.给出了系统可能呈现的四种运动状态的判断条件和运动微分方程.通过二维参数分岔分析得到系统在(ω,l)参数平面内存在的各类周期振动的参数域和分布规律.详细分析了相邻p/1类基本振动的转迁规律.在5/1基本振动的参数域的右边区域,相邻p/1基本振动的参数域临界线上存在一个奇异点X_p,相邻p/1类基本振动的分岔特点以奇异点X_p为临界点.在l小于l_X_p的区域内,相邻p/1基本振动经实擦边分岔和鞍结分岔相互转迁,实擦边分岔线和鞍结分岔线之间存在迟滞域,迟滞域内,系统存在两个周期吸引子共存的现象.在l大于l_X_p的区域内,相邻p/l类基本振动的参数域之间存在一个中间过渡区域.中间过渡区域内,系统呈现(2p+2)/2和(2p+1)/2周期振动等.在5/1基本振动的参数域的左边区域,p/1基本振动经多重滑移分岔产生(P+1)/1基本振动.     

Upscaling in Subsurface Transport Using Cluster Statistics of Percolation

Transport/flow problems in soils have been treated in random resistor network representations (RRNs). Two lines of argument can be used to justify such a representation. Solute transport at the pore-space level may probably be treated using a system of linear, first-order differential equations describing inter-pore probability fluxes. This equation is equivalent to a random impedance network representation. Alternatively, Darcys law with spatially variable hydraulic conductivity is equivalent to an RRN. Darcys law for the hydraulic conductivity is applicable at sufficiently low pressure head in saturated soils, but only for steady-state flow in unsaturated soils. The result given here will have two contributions, one of which is universal to any linear conductance problem, i.e., requires only the applicability of Darcys (or Ohms) law. The second contribution depends on the actual distribution of linear conductances appropriate. Although nonlinear effects in RRNs (including changes in resistance values resulting from current, analogous to changes in matric potential resulting from flow) have been treated within the framework of percolation theory, the theoretical development lags the corresponding development of the linear theory, which is, in principle, on a solid foundation. In practice, calculations of the nonlinear conductivity in relatively (compared with soils) well characterized solid-state systems such as amorphous or impure semiconductors, do not agree with each other or with experiment. In semiconductors, however, experiments do at least appear consistent with each other.In the limit of infinite system size the transport properties of a sufficiently inhomogeneous medium are best calculated through application of critical rate analysis with the system resistivity related to the critical (percolating) resistance value, R_c. Here well-known cluster statistics of percolation theory are used to derive the variability, W (R,x) in the smallest maximal resistance, R of a path spanning a volume x³ as well as to find the dependence of the mean value of the conductivity, (x). The functional form of the cluster statistics is a product of a power of cluster size, and a scaling function, either exponential or Gaussian, but which, in either case, cuts off cluster sizes at a finite value for any maximal resistance other than R_c. Either form leads to a maximum in W (R,x) at R=R_c. When the exponential form of the cluster statistics is used, and when individual resistors are exponential functions of random variables (as in stochastic treatments of the unsaturated zone by the McLaughlin group [see Graham and MacLaughlin (1991), or the series of papers by Yeh et al. (1985, 1995), etc.], or as is known for hopping conduction in condensed matter physics), then W (R,x) has a power law decay in R/R_c (or R_c/R, the power being an increasing function of x. If the statistics of the individual resistors are given by power law functions of random variables (as in Poiseiulles Law), then an exponential decay in R for W (R,x) is obtained with decay constant an increasing function of x. Use, instead, of the Gaussian cluster statistics alters the case of power law decay in R to an approximate power, with the value of the power a function of both R and x.     

The Singular Nature of Emergent Levels: Suggestions for a Theory of Emergence

Suggestions are offered for a theory of emergence based on a clarification and new interpretation of the singular nature of emergent levels. These suggestions cover formalisms, formulations, and measurements. In contrast to mere collectivities, as well as the rendering of macro- and micro-levels in entropy formulations, order parameters, and distinctions in temporal dynamics, emergent levels are described as privileged and confounded. A discussion of the insufficiency of previous formalisms in dealing with the structural novelty of emergent levels sets the stage for the introduction of a new formal construct, that of self-transcending constructions. This construct is linked to the idea of logical depth as a complexity measure. The advantages of a semantic rather than information–theoretic perspective are discussed. In addition, the tendency to confuse levels in models with levels in emergent phenomena themselves is described. Finally, conclusions about emergent levels as a new natural kind construct are offered.     

The Three Dimensional Viscous Camassa–Holm Equations, and Their Relation to the Navier–Stokes Equations and Turbulence Theory

We show here the global, in time, regularity of the three dimensional viscous Camassa–Holm (Navier–Stokes-alpha) (NS-) equations. We also provide estimates, in terms of the physical parameters of the equations, for the Hausdorff and fractal dimensions of their global attractor. In analogy with the Kolmogorov theory of turbulence, we define a small spatial scale, , as the scale at which the balance occurs in the mean rates of nonlinear transport of energy and viscous dissipation of energy. Furthermore, we show that the number of degrees of freedom in the long-time behavior of the solutions to these equations is bounded from above by (L/ )³, where L is a typical large spatial scale (e.g., the size of the domain). This estimate suggests that the Landau–Lifshitz classical theory of turbulence is suitable for interpreting the solutions of the NS- equations. Hence, one may consider these equations as a closure model for the Reynolds averaged Navier–Stokes equations (NSE). We study this approach, further, in other related papers. Finally, we discuss the relation of the NS- model to the NSE by proving a convergence theorem, that as the length scale ₁ tends to zero a subsequence of solutions of the NS- equations converges to a weak solution of the three dimensional NSE.     

紊动流场中悬浮颗粒分布的随机理论

通过分析固体颗粒在紊动流场中的随机运动,建立了二维流场中垂直于时均流动的方向上颗粒随机位移的概率密度分布函数所满足的方程。由该方程解出的分布函数在一定条件下即相当于颗粒浓度分布函数。运用这一方法研究了[1]、[2]中报道的壁面附近颗粒浓度降低的现象。     

On the approximate solution of elastic contact problems for a circular annulus

Paper concerns the Boussinesq and Reissner-Sagoci contact problems for a rigid punch in contact with the plane face of an isotropic elastic half-space over a circular annulusbra. It is shown that a simple combination of the interior and exterior Dirichlet and Neumann solutions for a circle, when expressed in terms of oblate spheriodal coordinates, yields extremely accurate approximate solutions provided that 0b/a0.8.