强旋湍流气－固两相流动的颗粒随机轨道法模拟

应用颗粒随机轨道模型，并与一种新的代数Ｒｅｙｎｏｌｄｓ应力模型相结合，对新型煤粉涡旋燃烧炉内强旋湍流气。固两相流动进行了数值模拟。得到了与实验相符合的颗粒相密度分布和质量流分布。计算结果表明，在涡旋燃烧炉内的强旋湍流流场中，外壁附近颗粒浓度最高，颗粒停留时间加长，气－固两相间滑移速度增大。     

强旋湍流气－固两相流动的颗粒随机轨道法模拟

应用颗粒随机轨道模型，并与一种新的代数Ｒｅｙｎｏｌｄｓ应力模型相结合，对新型煤粉涡旋燃烧炉内强旋湍流气。固两相流动进行了数值模拟。得到了与实验相符合的颗粒相密度分布和质量流分布。计算结果表明，在涡旋燃烧炉内的强旋湍流流场中，外壁附近颗粒浓度最高，颗粒停留时间加长，气－固两相间滑移速度增大。     

用PDF方程法分析悬沙垂线浓度分布

从颗粒运动的PDF(概率密度分布函数)输运方程出发。建立颗粒相的质量、动量和脉动速度二阶矩方程．对于明渠二维恒定均匀流。利用垂向动量方程导出了新的泥沙扩散方程。表明颗粒脉动强度梯度、升力、重力沉降和紊动扩散都影响悬沙运动。说明了传统扩散方程的不足．理论分析了水沙两相物理属性和水流条件对泥沙扩散系数和浓度分布特征的影响。并通过细颗粒试验资料的分析进行了定量研究．     

二阶矩两相湍流模型中的改进算法

建立了融合SIMPLEC算法在内的考虑了浓度修正影响的颗粒压力修正方程。提出了二阶矩湍流模型下考虑浓度修正值影响的两相湍流流动的算法,并将它和没有考虑浓度修正值影响的二阶矩湍流模型进行了对照。结果表明在二阶矩模型中是否考虑浓度修正影响会不同程度地影响流场的速度以及浓度等参数分布,考虑浓度修正影响的二阶矩湍流模型更能有效地预测稠密两相湍流流动。     

采用精确化理论求解厚板任意形开孔动应力集中

摘要：本文基于复变函数与保角映射法,采用平板弯曲振动精确化方程[9],对含任意形开孔平板中弹性波散射与动应力集中问题进行了研究。利用正交函数展开的方法将待解的问题归结为对一组无穷代数方程组的求解。作为算例,计算了自由边界条件下圆孔和椭圆孔的动弯矩集中系数的数值结果,并对板厚与孔径比对动弯矩分布的影响做了分析研究。结果表明：入射波数、平板厚度和椭圆偏心率等参数对动弯矩的分布都有很大的影响。在较低频率和平板较薄的情况下,基于文献[9]的方程与基于Mindlin板的动弯矩结果在数值分布上是基本一致的;在较高频率和平板较厚的情况下,基于文献[9]的方程与基于Mindlin板的动弯矩结果在数值分布计算结果相差较大。由于文献[9]给出的平板振动精确化方程是在没有任何工程假设条件下得到的,因此本文的分析计算结果更精确一些。     

极性材料的热弹性力学

本文在文献[1]、[2]、[3]、[5]、[7]的基础上,讨论了线性耦合下,极性材料三维热弹性力学的能量方程、熵生成定律、耗散函数、Fourier热传导方程和它的变分原理。     

最小二乘边界配点法解二维问题探讨

1.引言薄板弯曲问题和弹性力学平面问题的边界配点法,通常都是选用双调和方程的各类特解序列来构造试函数。在弹性力学教本中出现过的三类特解序列都可用来构造试函数。文献[2]以多项式级数为试函数,使用混合法分析过薄板强度。文[3]在边界配点法中使用了两种试函数——多项式级数和双曲三角函数。文[4]提出了一种产生双调和多项式的方法。文[5,6]使用的是双曲三角函数。文[7]综合了极坐标系中各种     

旋转直管内气固或液固两相流动浓度分布

本文针对旋转直管内气固或液固两相流动，建立了经过适当简经处理的颗粒无量纲运动方程，得到了稳定状态下直管内颗粒的浓度分布公式，从而为气固两相流风机或液固两相流泵的颗粒浓度分布研究打下了基础。     

弹性薄圆板的后屈曲分析

本文对薄圆板的后屈曲进行了研究。采用Galerkin法,试函数选为Legendre多项式,控制方程是Von-karman大挠度方程。考虑了简支,夹支两种边界条件。计算结果与有关文献[1]进行了比较,表明以Legendre多项式为试函数收敛快,精度高,且计算工作量较文献[1]为小。     

雾层两相耦合流场及气溶胶颗粒物的动力学性质

雾层气溶胶系统涉及复杂的动力学演变过程:碰撞、凝并、破碎、冷凝/蒸发、成核、沉积、表面化学反应等.因此,发展雾层与气相流场耦合的Eulerian-Lagrangian两相流模型、颗粒动力学及随机轨道模型,考虑重力、曳力、布朗力、Basset力等对颗粒相的作用.基于SIMPLE和多重Monte Carlo算法求解颗粒群平衡方程,自行开发了FAD程序首先对室内燃烧源细微颗粒物的扩散实验展开数值模拟,计算结果与实验数据吻合较好.将建立的模型和方法数值研究气溶胶污染物在雾环境中的输运过程,分析雾消散阶段颗粒相浓度、平均尺度的时空分布.结果显示:当时间演化至60 min,雾滴的平均尺度减小到初始的65.67％,而气溶胶颗粒最大数目对应的尺度为0.006 μm.     

The distribution of solid particles suspended in a turbulent flow: A stochastic approach

Basic fluid mechanics and stochastic theories are applied to show that the concentration distribution of suspended solid particles in a direction normal to the mean streamlines of a two-dimensional turbulent flow is greatly influenced by the lift force exerted on them in the vicinity of the wall. Analytic solution shows that, when the direction of the mean flow is horizontal, the probability density functionp (y, t) for random displacements of the particles will have a maximum value at a point from the wall where the perpendicular component of the lift force precisely balances particle gravity. Interpretation of experimental observations is presented using this theory.     

风沙两相流跃移层中沙粒相的速度分布

从单个跃移沙粒在气流中的运动方程出发导出了风沙两相流中沙粒相速度分布函数的Boltzmann方程，对风沙流研究中几种不同的分布函数及其相应的统计平均值等基本概念给出了严密的数学定义，指出了不同分布函数之间的区别和联系，在略去铅垂方向空气阻力的情况下，给出了沙粒相速度分布函数沿铅垂方向的边缘分布，作为风沙流中跃移理论的主要基础之一。利用结果对前人在风沙流研究中发现的某些重要规律和现象进行了解释。     

Particle collisions in a turbulent flow

An equation for the two-point probability density function of the two-particle the coordinate and velocity distribution is obtained. A closed system of equations for the first and second two-point moments of the velocity fluctuations of a pair of particles with allowance for the turbulent flow inhomogeneity is given. Boundary conditions for the equations of the particle concentration and the intensity of the relative random velocity during particle collision are obtained. A unified formula describing the interparticle collision process as a result of turbulent motion and the average relative particle velocity slip is obtained for the kernel of the coagulation equation. The effect of the average velocity slip of the particles and the carrier phase on the parameters of motion of the dispersed admixture and its coagulation is investigated on the basis of a two-point two-time velocity fluctuation autocorrelation function with two time and space scales representing the energy-bearing and small-scale motion of the fluid phase.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 104–116, March–April, 1996.     

Particle deposition on a channel wall in a turbulent disperse flow with gradient

Kinetic ideas about the motion of a set of particles (droplets) in a turbulent gas flow with gradient are used to derive a Fokker-Planck equation for the case of sufficiently large particles (more than few microns). This equation describes the process in which they are deposited on the wall of a channel. Satisfactory agreement has been obtained between the numerical solution to this equation for the deposition rate and the experimental data published in the literature. Under the assumption that the parameters of the carrier gaseous flow vary fairly slowly, a generalized equation is derived for particle diffusion in turbulent flow. This takes into account the intensity gradient of transverse pulsations in the velocity of the carrier gaseous flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 57–63, July–August, 1985.     

Mean and fluctuating components of drag and lift forces on an isolated finite-sized particle in turbulence

We perform fully resolved direct numerical simulations of an isolated particle subjected to free-stream turbulence in order to investigate the effect of turbulence on the drag and lift forces at the level of a single particle, following Bagchi and Balachandar’s work (Bagchi and Balachandar in Phys Fluids 15:3496–3513, 2003). The particle Reynolds numbers based on the mean relative particle velocity and the particle diameter are Re?=?100, 250 and 350, which covers three different regimes of wake evolution in a uniform flow: steady axisymmetric wake, steady planar symmetric wake, and unsteady planar symmetric vortex shedding. At each particle Reynolds number, the turbulent intensity is 5–10% of the mean relative particle velocity, and the corresponding diameter of the particle is comparable to or larger than the Kolmogorov scale. The simulation results show that standard drag values determined from uniform flow simulations can accurately predict the drag force if the turbulence intensity is sufficiently weak (5% or less compared to the mean relative velocity). However, it is shown that for finite-sized particles, flow non-uniformity, which is usually neglected in the case of the small particles, can play an important role in determining the forces as the relative turbulence intensity becomes large. The influence of flow non-uniformity on drag force could be qualitatively similar to the Faxen correction. In addition, finite-sized particles at sufficient Reynolds number are inherently subjected to stochastic forces arising from their self-induced vortex shedding in addition to lift force arising from the local ambient flow properties (vorticity and strain rate). The effect of rotational and strain rate of the ambient turbulence seen by the particle on the lift force is explored based on the conditional averaging using the generalized representation of the quasi-steady force proposed by Bagchi and Balachandar (J Fluid Mech 481:105–148, 2003). From the present study, it is shown that at Re?=?100, the lift force is mainly influenced by the surrounding turbulence, but at Re = 250 and 350, the lift force is affected by the wake structure as well as the surrounding turbulence. Thus, for a finite-sized particle of sufficient Reynolds number supporting self-induced vortex shedding, the lift force will not be completely correlated with the ambient flow. Therefore, it appears that in order to reliably predict the motion of a finite-sized particle in turbulence, it is important to incorporate both a deterministic component and a stochastic component in the force model. The best deterministic contribution is given by the conditional average. The influence of ambient turbulence at the scale of the particle, which are not accounted for in the deterministic contribution, can be considered in stochastic manner. In the modeling of lift force, additional stochastic contribution arising from self-induced vortex shedding must also be included.     

Physical interrelation between Fokker-Planck and random walk models with application to Coulomb interactions

A model of the random walk is formulated to allow a simple computing procedure to replace the difficult problem of solution of the Fokker-Planck equation. The step sizes and probabilities of taking steps in the various directions are expressed in terms of Fokker-Planck coefficients. Application is made to many particle systems with Coulomb interactions. The relaxation of a highly peaked velocity distribution of particles to equilibrium conditions is illustrated. Part I of a thesis submitted in partial fulfillment of Doctor of Philosophy requirements at Colorado State University.     

分布式运动约束下悬臂输液管的参数共振研究

输液管道结构在航空、航天、机械、海洋、水利和核电等工程领域都有广泛应用,其稳定性、振动与安全评估备受关注.针对具有分布式运动约束悬臂输液管的非线性动力学模型,分别采用立方非线性弹簧和修正三线性弹簧来模拟运动约束的作用力,研究了管道在脉动内流激励下的参数共振行为.首先,从输液管系统的非线性控制方程出发,利用Galerkin方法进行离散化;然后,由Floquet理论得出线性系统在失稳前两个不同平均流速下脉动幅值和脉动频率变化时的共振参数区域;最后,考虑系统的几何非线性项和分布式非线性运动约束力的影响,求解了管道的非线性动力学响应,讨论了非线性项及运动约束力对管道参数共振行为的影响.研究结果表明,系统非线性共振响应的参数区域与线性系统的共振参数区域是一致的,分布式运动约束力对发生参数共振时管道的位移响应有显著影响;立方非线性弹簧和修正三线性弹簧模型所预测的分岔路径存有较大差异,但都可诱发管道在一定的参数激励下出现混沌运动.      

The Dilute Rheology of Swimming Suspensions: A Simple Kinetic Model

A simple kinetic model is presented for the shear rheology of a dilute suspension of particles swimming at low Reynolds number. If interparticle hydrodynamic interactions are neglected, the configuration of the suspension is characterized by the particle orientation distribution, which satisfies a Fokker-Planck equation including the effects of the external shear flow, rotary diffusion, and particle tumbling. The orientation distribution then determines the leading-order term in the particle extra stress in the suspension, which can be evaluated based on the classic theory of Hinch and Leal (J Fluid Mech 52(4):683–712, 1972), and involves an additional contribution arising from the permanent force dipole exerted by the particles as they propel themselves through the fluid. Numerical solutions of the steady-state Fokker-Planck equation were obtained using a spectral method, and results are reported for the shear viscosity and normal stress difference coefficients in terms of flow strength, rotary diffusivity, and correlation time for tumbling. It is found that the rheology is characterized by much stronger normal stress differences than for passive suspensions, and that tail-actuated swimmers result in a strong decrease in the effective shear viscosity of the fluid.     

A Novel Framework for Integration of Random-Walk Particle-Tracking Simulation in Subsurface Multi-phase Flow Modeling

Coarse-scale models are generally preferred in the numerical simulation of multi-phase flow due to computational constraints. However, capturing the effects of fine-scale heterogeneity on flow and isolating the impacts of numerical (artificial) dispersion, which increases with scale, are not trivial. In this paper, a particle-tracking method is devised and integrated in a scale-up workflow to estimate the conditional probability distributions of multi-phase flow functions, which can be considered as inputs in coarse-scale simulations with existing commercial packages. First, a novel particle-tracking method is developed to solve the saturation transport equation. The transport calculation is coupled with a velocity update, following the implicit pressure, explicit saturation framework, to solve the governing equations of two-phase immiscible flow. Each phase particle is advanced in a deterministic convection step according to the phase velocity, as well as in a stochastic dispersion step based on the random Brownian motion. A kernel-based formulation is proposed for computation of fluid saturation in accordance with the phase particle distribution. A novel aspect is that this method employs the kernel approach to construct saturation from phase particle distribution, which is an important improvement to the conventional box method that necessitates a large number of particles per grid cell for consistent saturation interpolation. The model is validated against various analytical solutions. Finally, the validated model is integrated in a statistical scale-up procedure to calibrate effective, or “pseudo,” multi-phase flow functions (e.g., relative permeability functions). The proposed scale-up framework does not impose any length scale requirement regarding the distribution of sub-grid heterogeneities.