饱和砂土中泥浆渗透的变形-渗流-扩散耦合计算模型

传统的泥浆渗透计算中没有考虑土体变形和浆液流速的影响.根据泥浆颗粒的质量守恒定律推导了耦合流速的浓度扩散方程,并通过在浓度方程中引入沉积系数进一步计算得到沉积颗粒的质量;同时,以沉积量作为耦合项对毕奥固结方程中的水量连续方程进行了修正,在此基础上建立了变形-渗流-扩散耦合的控制方程及其变分原理. 采用有限单元法求解基本方程,运用了时间增量法与直接迭代法,并利用一维试验验证计算方法的可靠性,并与赫齐格的经典模型的计算结果进行了比较,结果表明,本文建立的模型的计算结果可以较好地预测各组试验中颗粒的沉积规律,且吻合程度优于仅考虑颗粒对流和扩散的传统计算方法. 最后,将泥浆在槽壁中的渗透简化为二维问题并进行了计算,计算结果与工程认识相符合,泥浆的沉积填充效应随深度的增加而增大,施工时需要严格控制浅层作业段的机械垂直度;成槽机的下斗抓挖时机可以根据地层填充的致密程度进行计算,对现场施工具有一定的指导意义.      

谱元方法求解环形空间内自然对流换热

提出了将谱元方法应用到极坐标系下，利用极坐标系下的谱元方法求解环形空间内自然对流问题。具体求解了原始变量速度和压力的不可压缩Navier-Stokes方程和能量方程，通过在时间方向采用时间分裂方法和空间采用谱元方法对方程进行离散求解，取得了与基准解较一致的计算结果。     

A comparative study of high‐order variable‐property segregated algorithms for unsteady low Mach number flows

In this paper, five different algorithms are presented for the simulation of low Mach flows with large temperature variations, based on second‐order central‐difference or fourth‐order compact spatial discretization and a pressure projection‐type method. A semi‐implicit three‐step Runge–Kutta/Crank–Nicolson or second‐order iterative scheme is used for time integration. The different algorithms solve the coupled set of governing scalar equations in a decoupled segregate manner. In the first algorithm, a temperature equation is solved and density is calculated from the equation of state, while the second algorithm advances the density using the differential form of the equation of state. The third algorithm solves the continuity equation and the fourth algorithm solves both the continuity and enthalpy equation in conservative form. An iterative decoupled algorithm is also proposed, which allows the computation of the fully coupled solution. All five algorithms solve the momentum equation in conservative form and use a constant‐ or variable‐coefficient Poisson equation for the pressure. The efficiency of the fourth‐order compact scheme and the performances of the decoupling algorithms are demonstrated in three flow problems with large temperature variations: non‐Boussinesq natural convection, channel flow instability, flame–vortex interaction. Copyright © 2010 John Wiley & Sons, Ltd.     

The solution of a convection–diffusion–reaction equation arising from cathodic reduction in a pulsating crack

This paper derives the convection–diffusion-reaction equation governing the reaction between the dissolved oxygen in sea-water and the steel walls of a pulsating crack. By the neglect of the diffusion term it is shown that an exact solution of the convection-reaction equation can be obtained. A numerical method for the solution of the complete convection–diffusion-reaction equation is derived by the use of finite differences. The numerical computation of the initial transient and the final periodic steady-state values is also discussed.     

DETERMINATION OF DISPLACEMENT OF A BEAM WITH ONE END OVERHANGING IN SYMMETRIC BENDING WITH CONVERSION METHOD

置换法应用于求解一端外伸梁,在对称弯曲的条件下,根据直梁挠曲线所在平面内其与切线所成图形的边角几何关系,推导出求解该形梁的挠度和转角的置换法位移方程,其变量是相应的置换梁自由端的挠度、梁长、梁轴线位置坐标等. 对具体载荷梁的求解过程是：先以具体量值填充左、右置换梁自由端的挠度,再将其代入该置换法位移方程的统一表达式,即得到所求梁段的挠度、转角的方程全解. 所用的计算为代数方程的分式四则运算,只需挠曲线和叠加原理概念,无需积分,一般无需查挠度表,结果精确. 给出工程背景的算例.     

混合层流动拟序结构的大涡模拟

采用大涡模拟方法对空间发展的二维平面混合层进行了数值模拟 ,动量方程采用分步投影法求解 ,亚格子项采用标准Smagorinsky亚格子模式模拟 ,压力泊松方程采用修正的循环消去法快速求解 ,同时求解了标志物输运方程以实现数值流场显示。模拟结果给出了混合层流动的瞬态发展过程以及流动中拟序结构的发展演变过程 ,成功地模拟了混合层发展中的各种瞬态细节过程 ,如涡的卷起、增长 ,涡与涡之间的配对、合并过程 ,以及大涡破碎为小涡的级联过程 ,为各种以混合层流动为原型流动的射流、尾流等工业流动的控制和优化提供了理论基础。     

Numerical solution of the incompressible Navier–Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier–Stokes equations is presented. The scheme is constructed on a staggered‐mesh grid system. The convection terms are discretized with a fifth‐order‐accurate upwind compact difference approximation, the viscous terms are discretized with a sixth‐order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth‐order difference approximation on a cell‐centered mesh. Time advancement uses a three‐stage Runge–Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented. Copyright © 1999 John Wiley & Sons, Ltd.     

ANALYTICAL SOLUTION OF THERMAL CONSOLIDATION OF SATURATED SOILS WITH TWO PARALLEL CYLINDRICAL HEAT SOURCES1）

基于应力平衡条件、渗流连续方程、能量守恒方程,考虑土颗粒和孔隙水热膨胀系数的不同,建立考虑热水力耦合的饱和土体三维热固结控制方程。利用傅里叶变换和拉普拉斯变换导出变换域上的控制方程,解得点热源在变换域和实数域上的解析解,再利用区域积分给出两平行圆柱形热源热固结土体温度、孔压、位移的解析解,并对其进行分析,发现径距比增大会导致两热源温度相互影响程度减弱,热固结系数减小会导致孔压和位移的峰值增大。     

SOME PROBLEMS WITH THE METHOD OF FUNDAMENTAL SOLUTION USING RADIAL BASIS FUNCTIONS

The present work describes the application of the method of fundamental solutions (MFS) along with the analog equation method (AEM) and radial basis function (RBF) approximation for solving the 2D isotropic and anisotropic Helmholtz problems with different wave numbers. The AEM is used to convert the original governing equation into the classical Poisson's equation, and the MFS and RBF approximations are used to derive the homogeneous and particular solutions, respectively. Finally, the satisfaction of the solution consisting of the homogeneous and particular parts to the related governing equation and boundary conditions can produce a system of linear equations, which can be solved with the singular value decomposition (SVD) technique. In the computation, such crucial factors related to the MFS-RBF as the location of the virtual boundary, the differential and integrating strategies, and the variation of shape parameters in multi-quadric (MQ) are fully analyzed to provide useful reference.     

NUMERICAL SIMULATION OF BROWNIAN COAGULATION AND MIXING OF NANOPARTICLES IN 2D RAYLEIGH-B?NARD CONVECTION

采用泰勒展开矩方法对二维瑞利-贝纳德热对流系统(1×10⁶ ≤Ra ≤1 ×10⁸) 中纳米颗粒群的混合和凝并特性进行了数值模拟. 结果显示颗粒群随时间演化经历了扩散阶段、混合阶段、充分混合阶段3 个阶段, 随着颗粒群混合和凝并的进行, 颗粒数目浓度减少, 颗粒群的平均体积增大; 得到了颗粒分布函数各特征量与温度相关系数以及各特征量的空间分布标准偏差在3 个阶段的不同特征; 得到了颗粒分布函数各阶矩以及平均体积长时间演化的渐近行为, 结果与零维渐近解析解一致. 最后, 本文进一步研究了无量纲数(包括瑞利数Ra, 斯密特数Sc_M, 达姆科勒数Da) 对颗粒群达到自保持分布时间的影响, 发现该时间随着Ra和Sc_M的增大呈对数率减小, 随着Da的增大呈线性增大     

Large eddy simulation of natural convection along a vertical isothermal surface

Large eddy simulations of natural convection along a vertical isothermal surface have been carried out using a parallel CFD code SMAFS (Smoke Movement And Flame Spread) developed by the first author to study the dynamics of the natural convection flow and the associated convective heat transfer, with sub-grid scale turbulence modeled using the Smagorinsky model. In the computation, the filtered governing equations are discretized using finite volume method, with the variables at the cell faces in the finite volume discrete equations approximated by a second order bounded QUICK scheme and the diffusion term computed based on central difference scheme. The computation was time marched explicitly, with momentum equations solved based on a second order fractional-step Adams–Bashford scheme and enthalpy computed using a second order Runge–Kutta scheme. The Poisson equation for pressure from the continuity equation was solved using a multi-grid solver. The results including the temperature and velocity profiles of the boundary layer and the local heat transfer rate are analyzed. Comparison is made with experimental data and good agreement is found.     

Instability of fully developed mixed convection with viscous dissipation in a vertical porous channel

Flow driven by an externally imposed pressure gradient in a vertical porous channel is analysed. The combined effects of viscous dissipation and thermal buoyancy are taken into account. These effects yield a basic mixed convection regime given by dual flow branches. Duality of flow emerges for a given vertical pressure gradient. In the case of downward pressure gradient, i.e. upward mean flow, dual solutions coincide when the intensity of the downward pressure gradient attains a maximum. Above this maximum no stationary and parallel flow solution exists. A nonlinear stability analysis of the dual solution branches is carried out limited to parallel flow perturbations. This analysis is sufficient to prove that one of the dual solution branches is unstable. The evolution in time of a solution in the unstable branch is also studied by a direct numerical solution of the governing equation.     

一种新型高鲁棒性动网格技术及其应用

首先对四元数进行李代数空间指数映射,解决了多个四元数插值问题,并结合距离倒数插值方法实现网格边界扰动向空间网格的传播,建立了新型高鲁棒性的四元数变形网格技术. 针对该型动网格技术中由于大型矩阵运算量引起的运算效率低问题,同时利用四元数方法在动网格变形中具备与物面边界高阶一致性的特点,提出了分层次变形策略,避免了面向全流场网格节点的大型矩阵运算;进一步基于无限插值技术较强的逻辑保持能力,建立了面向结构网格分层混合变形方法. 充分利用多区域重叠、对接网格变形技术中隐含的并行性,基于对等式编程思想及MPI 库函数对动网格程序进行并行化编程,建立了高效高鲁棒性的变形网格技术. 以某型客机翼身组合体气动弹性分析为范例,研究了不同方法之间的计算效率以及鲁棒性,进一步将分层混合变形网格技术应用于某型支线客机全机型架外形设计与修正,验证了所建立的新型动网格技术的高效性与鲁棒性.      

气体动理学格式与多尺度流动模拟

介绍了气体动理学格式(GKS)的基本构造原理及其在两种典型多尺度流动模拟中的应用。GKS利用介观BGK方程的跨尺度演化解来构造网格界面上的数值通量,从而发展出能随计算网格尺度变化自动切换物理模型的多尺度方法。对湍流这种宏观多尺度流动,发展了高精度GKS方法并成功用于低雷诺数湍流的直接数值模拟;为实现对高雷诺数湍流的高效精细模拟,基于拓展BGK方程和已有的RANS,LES模型建立了新型多尺度模拟框架。对跨流域稀薄流动,发展了适合大规模并行的三维统一气体动理学格式(UGKS),并建立了适合轴对称稀薄流动的UGKS。研究表明,GKS在多尺度流动高效模拟中的优异性能,具有很好的发展前景。     

Nodal operator splitting adaptive finite element algorithms for the Navier–Stokes equations

Solution algorithms for solving the Navier–Stokes equations without storing equation matrices are developed. The algorithms operate on a nodal basis, where the finite element information is stored as the co-ordinates of the nodes and the nodes in each element. Temporary storage is needed, such as the search vectors, correction vectors and right hand side vectors in the conjugate gradient algorithms which are limited to one-dimensional vectors. The nodal solution algorithms consist of splitting the Navier–Stokes equations into equation systems which are solved sequencially. In the pressure split algorithm, the velocities are found from the diffusion–convection equation and the pressure is computed from these velocities. The computed velocities are then corrected with the pressure gradient. In the velocity–pressure split algorithm, a velocity approximation is first found from the diffusion equation. This velocity is corrected by solving the convection equation. The pressure is then found from these velocities. Finally, the velocities are corrected by the pressure gradient. The nodal algorithms are compared by solving the original Navier–Stokes equations. The pressure split and velocity–pressure split equation systems are solved using ILU preconditioned conjugate gradient methods where the equation matrices are stored, and by using diagonal preconditioned conjugate gradient methods without storing the equation matrices. © 1998 John Wiley & Sons, Ltd.     

Application of streamwise diffusion to time-dependent free convection of liquid metals

A numerical analysis is given for the application of streamwise diffusion to high-intensity flows with marginal spatial resolution. Terms are added to the momentum equation which are similar to those used in the Petrov-Galerkin, Taylor-Galerkin and balancing tensor diffusivity methods. Values for the streamwise viscosity are obtained from mesh refinement studies. An illustration is given for the time-dependent free convection of a liquid metal in a cavity with differentially heated sided walls. The spatial problem is solved with the Galerkin finite element method and the time integration is performed with the backward Euler method. Solution quality and computation time are compared for results with and without added streamwise diffusion. For the cases presented, streamwise diffusion eliminates spurious oscillations and saves computation time without compromising the solution.     

谐波平衡法在动导数快速预测中的应用研究

谐波平衡法以傅里叶级数展开为基础,将周期性非定常流场的非定常求解过程转化为几个定常流场的耦合求解过程,并通过重建得到整个流场的非定常过程. 建立了基于谐波平衡法的动导数快速预测方法,数值模拟了超声速带翼导弹俯仰的动态流场,并通过积分法获取了俯仰动导数,与实验结果吻合很好;且在同等计算精度下,谐波平衡法的计算效率是双时间步方法的13 倍. 应用谐波平衡法研究了较大范围内减缩频率对俯仰动导数的影响规律. 研究发现,对于本外形,当减缩频率降低到一定值后,俯仰动导数的值迅速变化,甚至发生变号;对此现象产生的原因进行了深入分析,并通过对导弹自激俯仰运动的数值模拟验证了该结果. 此外,针对大攻角条件下动态流场非线性强的特点,开展了谐波平衡法在大攻角下的适用性研究. 结果表明,谐波平衡法在大攻角下也能取得很好的计算结果.      

BIOT’S CONSOLIDATION ANALYSIS FOR CYLINDRICAL SOIL-SAMPLE BASED ON HANSBO’S FLOW1)

诸多黏性土渗流试验表明,在低水力梯度下渗流会出现明显偏离Darcy定律的现象. 为了分析渗流的非Darcy特性对固结过程 的影响,引入Hansbo渗流方程描述圆柱土样内的渗流,重新推导轴对称条件下的Biot固结方程,并给出方程的Crank--Nicolson有限 差分格式. 通过与Darcy渗流条件下轴对称Biot固结方程解析解的对比,验证计算方法的有效性. 然后分析Hansbo模型参数对圆柱 土样固结过程的影响. 计算结果表明：与Darcy渗流相比,Hansbo渗流会延缓圆柱土样的固结过程. 随着Hansbo渗流参数m或I¹的增大,在固结前期,Mandel--Cryer效应会更加显著,即孔隙水压力峰值将提高,且达到该峰值的时间 会延迟;在固结中后期,孔压消散滞后的现象也更加明显. 不过,Hansbo渗流对位移的影响很小.     

Rayleigh-Bénard湍流热对流研究的进展、现状及展望

对流现象广泛存在于恒星和行星里.在地球上,对流现象在诸如大气、海洋、地核和地幔等众多动力学系统中起着重要作用.Rayleigh-Bénard(RB)湍流热对流系统是从这些复杂的自然现象中抽象出来的研究对流问题的经典流体力学模型.本文主要从湍流传热、羽流、大尺度流动结构、速度和温度脉动的小尺度统计和非传统RB对流等几个方面着重评述近年来RB对流的若干研究新进展,并对今后的研究做出展望.     

Investigation of the turbulent model for pressure fluctuations with spectral theory

The pressure fluctuations in turbulent shear flows are investigated with the theory of spectral analysis.An expression for pressure spectra is analytically derived in terms of velocity spectra.This derivation is based on a formal solution of the Navier-Stokes equation and quasi-normal assumption to express the third and fourth order velocity correlations in terms of double velocity correlation.Then,a turbulent model for the computation of pressure fluctuation intensity with Renolds stress and mean flow velocity gradients is established.The turbulent constants in the model are calculated from the assumptions about the general behaviour of velocity spectra in high Renolds number flows.Comparison with direct simulation of turbulent boundary layer is made.It is found that the turbulent-turbulent,cross correlation,and turbulent-shear source terms for mean square value of pressure fluctuation are about the same magnitude.