带电沙粒的跃移云数值模型

给出吹过无限大平坦沙面的稳定风场中风、跃移沙粒及静电力相互耦合的方程 ,并用简化的沙粒——床面碰撞模型求解了方程 .计算结果表明静电力对沙粒跃移轨道及贴地风速有强烈影响 ,当跃移沙粒的带电量为 - 6 0 μC/ kg时 ,计算所得的跃移层中风速沿高程分布与实验结果相吻合     

形状因素对风沙跃移影响的数值研究

目的研究沙粒形状因素对沙粒跃移的影响。方法采用两类形状系数对阻力系数加以修正,进行了沙粒跃移的数值模拟。结果获得了形状系数修正后的模拟结果,并对比分析了其与采用球形假设的模拟结果以及真实沙粒跃移轨迹的高速摄影实验结果之间的差异。结论沙粒形状因素对沙粒跃移的影响是不能忽略的。当粒径>160μm时,两种修正方法得到的计算结果与实验轨迹基本一致,而当粒径≤160μm时,采用基于颗粒阻力的形状系数,将沙粒转化为球形的修正方法得到结果更接近实际。     

风沙跃移运动中的Magnus效应

给出了计及Magnus效应的风一沙耦合跃移运动模型,计算了当风沙流达到自平衡状态后Magnus效应对贴地风速廓线、沙粒跃移轨迹的影响,指出在Msagnus效应影响下,沙粒跃移运动对风速廓线的影响类似于粗糙床面对沙粒跃移的影响,其粗糙度约为沙粒跃移高度的0．8倍。     

跃移雪颗粒运动特性的数值模拟研究

运用拉格朗日方法,根据雪颗粒的运动机理,对单个跃移颗粒运行轨迹进行了数值计算.同时,结合颗粒的受力特征,对影响其运动的一些参数,如粒子直径、摩擦速度和颗粒密度等进行了分析.结果表明,随着雪颗粒直径或雪粒密度的增大,粒子竖向受力与重力之比减小,粒子的跃移高度则随之增大;而当粒子起跳速度不同时,水平方向受力与重力之比的变化规律也表现不同,进而影响其水平方向的跃移长度.基于经验公式,建立了雪面侵蚀质量通量与跃移粒子数目之间的关系,计算了雪颗粒质量传输率,并与经典文献的计算方法进行了比较.结果显示,随着雪颗粒直径和摩擦速度的增大,粒子质量传输率随之增大;随着阈值摩擦速度增大,粒子质量传输率随之减小.所得计算结果与Iversen经验公式吻合较好.     

风速脉动对跃移运动的影响

对风沙流中随机风场模型进行了描述,然后基于高斯概率分布,利用随机Fourier级数生成脉动风场,进而对混合粒径颗粒在湍流随机风场中的运动进行数值计算,给出了颗粒的垂直质量通量和水平质量通量以及沙粒粒径分布的变化情况。实例分析表明,所提出的算法在给定条件下,数值模拟结果与实际观测结果符合良好。     

沙粒的跃移与悬移

通过考虑沙粒在大气边界层中的湍流随机运动,分析了影响沙粒运动轨迹的各种因素．由于存在湍流脉动速度,因此在相同粒径、起跳初速度和摩阻风速情况下,沙粒运动轨迹会出现较大差异．沙粒的垂直波动速度可以作为一个反映其运动性质的参数,不仅与粒径、起跳初速度,而且与摩阻风速有关．利用其作为判断沙粒运动形式的判据,计算了给定沙粒粒径和摩阻风速条件下沙粒进入跃移和悬移的临界起跳初速度及相应运动形式的组分,并与考虑风沙耦合作用与不考虑风沙耦合作用的轨迹统计结果进行了比较．     

风沙电多场耦合中沙粒跃移运动的受力分析

基于风沙电多场耦合跃移运动的理论模型,通过对稳定风沙运动的数值模拟,将沙粒在运动过程中所受各种作用力,如空气阻力、Magnus力、Saffman力、静电力与重力进行了对比研究,讨论了沙粒所受各种作用力随摩阻风速、沙粒粒径、沙粒初始旋转角速度和初始起跳速度的变化规律．     

一种简便的静电场数值模拟方法

利用两个实例,我们阐述了对静电场做数值模拟的一种简便方法。     

跃移层中沙粒浓度分布特征的研究

对于不同的沙粒起跳初速度分布函数，利用数值方法模拟了沙粒浓度沿高程的分布规律，发现风沙强度沿高程具有“递增－饱和－衰减”的按层分布的规律，进一步详细研究了影响风沙流结构的主要因素，如风速分布、地表粗糙度、沙粒粒径等对风沙流结构的影响方式和规律。     

油藏水平井数值模拟研究中一种新的求解方法

提出了底水油藏水平井油水两相渗流问题的物理模型和数学模型，并建立了差分方程，然后利用ＩＭＰＥＳ求解方法的思路，并对其不足之处进行了改进，从而形成了一种新的求解方法．这种方法是采用上游网格饱和度插值计算相邻网格间的水量交换，用本阶段时间中点的饱和度对压力方程进行再迭代．保持了ＩＭＰＥＳ方法简单方便、占计算机内存少、计算速度快的优点，而精度明显提高，结果更趋合理．实际计算表明，新方法效果良好，可广泛用于底水、边水或人工注水油藏油水两相渗流问题的数值模型研究，对水平井和垂直井均适用．     

三维油资源渗流力学运移聚集的大规模数值模拟和应用

盆地发育的运移聚集史数值模拟,其功能是重建油气盆地的运移聚集演化史,它对于油资源的勘探,确定油藏位置和计算油藏贮量,寻找新的油田,具有极其重要的价值．从渗流力学和地质科学实际出发,研究了三维（带断层和通道）多层问题的渗流力学和石油地质非线性特征,提出一类迎风分数步隐式迭代格式,对胜利油田东营凹陷、惠民凹陷、滩海地区和阳信洼陷的实际问题进行数值模拟试验,结果和实际地质情况吻合．对模型问题,得到数值方法收敛性,成功解决了这一重要问题．     

塑料异型材挤出流动的一种数值模拟方法

提出一种基于横截面法和流动路径法的异型材挤出成型中熔体在模头内流动的数值模拟方法。此方法的主要特点是三维流动问题转变为一个一维和一个二维问题求解,它可预测截面复杂的异型材挤出模头压力降和具有多个分支流道的挤出模头出口处物料的平均流速,并用实际算例说明了数值模拟的特点和用途。     

三维风沙运动的CFD-DEM数值模拟

采用计算流体力学（Computational Fluid Dynamics,CFD）和颗粒离散元（Discrete Element Method,DEM）耦合的方法模拟三维风沙运动,并且将三维模拟结果和二维模拟结果以及实验结果进行了对比.计算结果表明：沙粒水平速度随着高度按幂函数规律增加,沙床表面附近沙粒撞击和起跳速度的概率分布均可用对数正态函数描述,沙粒撞击和起跳角度的概率分布均可用指数函数描述,沙粒水平速度、展向速度和垂直速度在不同高度处的概率分布可分别用对数正态分布、正态分布和正态分布表示.与二维计算结果的分析对比表明：二维计算得到的颗粒速度的分布规律和三维计算结果类似,但二维计算的颗粒表观密度明显偏大,由此导致输沙量计算偏大.和实验结果的对比表明：三维计算得到的颗粒速度概率分布与实验基本保持一致.     

三维潮流泥沙运动的一种数值模拟

采用σ坐标变换,将具有自由面的三维非线性潮流微分方程转换为模型所采用的模式方程;运用三角网络有限差分格式,将该模式方程和泥沙扩散方程分别离散,提出了一种新的三维潮流和泥沙运动的数值模拟方法,并与理论解进行了比较,对渤海湾湾顶典型玻潮流流与泥沙运动进行了数值模拟,以较高的辨率揭示了该海域潮流与泥沙运动的典型特征,计算结果与实际吻合良好。     

一种模拟覆冰双分裂导线舞动的数值方法

给出一种覆冰双分裂导线舞动的数值模拟方法。用Hamilton变分原理建立系统的动力学平衡方程,利用罚函数法引入子导线上间隔棒连接点的运动约束条件。采用Newmark法进行时间积分,Newton-Raphson迭代求解非线性方程,实现对双分裂导线舞动的数值模拟。数值模拟了垂直双分裂和水平双分裂覆冰导线的舞动,验证了方法的正确性。利用该方法能考虑由于迎风侧子导线尾流对背风侧子导线的影响而出现的作用在两子导线上的空气动力载荷不同的情况,为分裂导线舞动的深入研究提供了一种有效的途径。     

含夹矸厚煤层综放开采顶煤运移规律数值模拟

根据鹤岗富力矿含中位厚夹矸煤层地质条件及煤的赋存状况,应用RFPA2D数值模拟软件分别对-450南18-2工作面含0.5、1.0、1.5 m厚夹矸开采煤体的运移情况进行数值模拟。结果表明,采动过程中,上部顶煤垂直方向的位移量均大于下部,而下部顶煤始动点距工作面的距离均大于上部。夹矸层强度与煤层强度既定情况下,工作面顶煤的运移规律及破碎程度均与含夹矸层的厚度相关。含1.0 m厚夹矸顶煤总位移量大于含0.5、1.5 m厚夹矸顶煤总位移量;含1.0 m厚夹矸顶煤比含0.5、1.5 m厚夹矸顶煤破碎充分。该研究为含夹矸条件下厚煤层综放开采冒放工艺的确定提供了依据。     

污染物在受潮汐影响的地下水中输移的模拟

模拟沿海地区地下水及污染物的输移需要考虑海洋潮汐的影响,对于变密度污染物,考虑潮汐的影响会大大增加模型计算量,模拟此类问题要求模型具备优良的计算效果.为了获得适用于此类问题的模型,在比较各类模型的基础上,认为2DFEMFAT是合适的模型.对该模型作了修改完善,增加了潮汐模拟条件,使其能模拟潮汐影响下地下水中污染物的输移过程.为了验证修改后模型的正确性,作了潮汐条件下污染物输移的室内实验.数值模拟与实验数据的比较发 现,两者较为吻合,验证了模型的正确性和有效性.该模型可以推广应用到大尺度的实际问题.     

一种非线性动力学系统数值模拟的新方法

给出了应用ＭＡＴＬＡＢ语言动态仿真工具－ＳＩＭＵＬＩＮＫ对非线性动力学系统进行数值模拟的方法,以典型的非线性动力学系统－ｖａｎ ｄｅｒ Ｐｏｌ方程、Ｄｕｆｆｉｎｇ方程和Ｌｏｒｅｎｚ方程为例进行了说明。     

群桩挤土效应的数值模拟

以上海市江宁路551号大楼的打桩工程为实例,以有限元方法为主要分析手段,对群桩地表的隆起、桩周土体的侧移、挤土产生的应力及其对周围桩体的影响等挤土效应的变化规律进行了详细研究．     

复杂边界非饱和模型有限元数值模拟

利用有限元集中质量法对一维复杂边界土壤水分运动数学模型（非饱和Richards模型）进行了数值模拟,并进行了模型验证．模型验证结果表明,该方法计算精度很高,符合实际工作要求．