Ⅰ型拉伸试样金属注射成形多相流动的三维数值模拟

将粉末注射成形充模过程视为粘结剂、粉末和空气的三相流动过程,基于多相流理论给出金属注射成形的多相流动控制方程。根据金属注射成形的工艺特点确定多相流动控制方程的初边界条件,用多相流数值分析软件CFX对Ⅰ型拉伸试样的铁粉注射成形充模流动过程实现了三维数值模拟,分析了模腔中不同位置处粘结剂一粉末流动的速度曲线。数值模拟结果表明在Ⅰ型拉伸试样金属注射成形充模流动过程中发生过两相分离的现象,数值模拟的瞬态信息可用于金属注射成形产品缺陷的分析与控制。     

粉末注射成形动力方程的混沌特征分析及研究

粉末注射成形（PIM）充模过程具有混沌的特征．本文利用软件ANSYS对PIM充模过程的动力学方程数值求解。通过吸引子形态描述法分析对比两种不同注射速度下的动态吸引子图,计算混沌吸引子形态特征量,定量的表述混沌吸引子变化规律．研究表明粉末注射成形充模过程中存在混沌现象．     

求解二维不可压缩Navier-Stokes方程的混合型微分求积法

微分求积法（ＤＱＭ）能以较少的网格点求得微分方程的高精度数值解,但采用单纯的微分求积法求解二维不可压缩Ｎａｖｉｅｒ＿Ｓｔｏｋｅｓ 方程时,只能对低雷诺数流动获得较好的数值解,当雷诺数较高时会导致数值解不收敛·　为此,提出了一种微分求积法与迎风差分法混合求解二维不可压缩Ｎａｖｉｅｒ＿Ｓｔｏｋｅｓ 方程的预估＿校正数值格式,用伪时间相关算法以较少的网格点获得了较高雷诺数流动的数值解·　作为算例,对１∶１ 和１∶２ 驱动方腔内的流动进行了计算,得到了较好的数值结果·     

求解二维浅水方程的一种高分辨率有限体积法

本文将一种van Albada型可微的限制器函数引入到二维浅水方程的求解中，发展了一种求解二维浅水方程的有限体积法．数值实验结果表明，该方法不仅计算精度高，而且较其它求解二维浅水方程的高精度有限体积法，在数值解的收敛性能方面大有改善．     

二维Volterra积分方程数值解的渐近展开及其外推

本文讨论了求解二维Ｖｏｌｔｅｒｒａ积分方程的Ｎｙｓｔｒｏｍ方法，得到了数值解的逐项渐近展开，从而可进行Ｒｉｃｈａｒｄｓｏｎ外推，提高数值解的精度。     

二维Fredholm积分方程Nystrom方法的渐近展开及其外推

本文讨论了求解二维第二类Ｆｒｅｄｈｏｌｍ积分方程的Ｎｙｓｔｒｏｍ方法，得到了数值解的逐项渐近展开，从而可进行Ｒｉｃｈａｒｄｓｏｎ外推，提高数值解的精度。     

基于Picard迭代的PN×PN-2谱元法求解定常不可压缩Navier-Stokes方程

该文给出了一种求解二维定常不可压缩Navier-Stokes方程的基于Picard线性化迭代的P_N×P_N-2谱元法.通过Picard线性化将不可压缩Navier-Stokes方程的求解转化为一系列线性的Stokes-type方程,再利用非交错网格的P_N×P_N-2谱元法计算每个迭代步的Stokes-type方程.为了消除伪压力模,压力离散比速度离散低两阶,非交错网格的应用使得方程的离散方便且不会带来相应的插值误差,从而保证了谱精度.通过此方法数值计算了有精确解的Stokes流动、Kovasznay流动和方腔顶盖驱动流,结果表明,迭代收敛非常快,误差收敛达到了谱精度收敛,并且避免了压力震荡的出现,表明了该文方法准确可靠.     

二维非线性Volterra积分方程数值解的渐近展开及其外推

本文讨论了求解二维非线性Ｖｏｌｔｅｒｒａ积分方程的Ｎｙｓｔｒｏｍ方法，得到了数值解的逐项渐近展开。从而可进行Ｒｉｃｈａｒｄｓｏｎ外推，提高数值解的精度。     

直圆管突扩通道内宾汉流体湍流流场的数值研究

本文依据牛顿流体中建立的标准ｋ＿ε湍流模型这一基本思想，考虑宾汉流体的本构方程，建立了适用于求解宾汉流体湍流流动的控制方程·采用压力修正算法，实现了宾汉流体速度场与压力场的关联·在理论研究基础上，对直圆管突扩通道内宾汉流体湍流流动进行了数值研究，并探讨了直圆管突扩通道内宾汉流体湍流流动机理·     

基于外场实验的风力机叶片三维效应研究

针对一台33 kW水平轴风电机组开展了外场实验,得到其叶片7个断面翼型的压力分布曲线;基于求解时均N-S方程对风轮进行三维数值模拟,以及将叶片各断面作为二维翼型进行数值计算,分别得到各断面翼型的压力分布曲线及升阻力系数．通过将外场实验、三维和二维数值计算所得压力分布曲线及升阻力系数进行对比分析,研究了三维效应对风力机气动性能的影响．研究表明,从叶尖到叶根各断面翼型的压差先增大后逐渐减小,叶片表面压力分布曲线比较明显地反映了从叶尖到叶根流动分离的变化;叶片表面压力分布的三维数值计算结果较二维计算结果更加接近于外场实验值;风力机叶片表面的三维流动对叶片的气动性能影响较大,在叶尖和叶根部分尤为突出.     

Cavity flow simulation of Carreau–Yasuda non-Newtonian fluids using PIM meshfree method

Lid-driven cavity flow of a purely-viscous non-Newtonian fluid obeying Carreau–Yasuda rheological model is studied numerically using the PIM meshfree method combined with the Characteristic-Based Split-A algorithm. Results are reported for the velocity and pressure profiles at Reynolds numbers as high as 1000 for a non-Newtonian fluid obeying Carreau–Yasuda rheological model. For Newtonian fluids, results obtained from our PIM–CBS-A method show good agreement with benchmark results published in the literature and obtained using finite difference and/or Finite Element Methods. Our numerical results are also consistent with recent published results obtained using another meshfree method called LSM with the advantage that PIM needs less “points” to achieve the same degree of accuracy. Results obtained for the Carreau–Yasuda model reveals the strong effect of the shear-thinning behavior of a fluid on its flow kinematics within the cavity.     

精确计算Bingham流体圆管层流压降及速度分布规律的新方法

Bingham(宾汉)模型情况下,多采用通用公式进行圆管层流压降的解析计算,即将Bingham模型本构方程代入粘性流体圆管层流流动通用公式进行计算,仅能得到压降的解析解.新方法结合Bingham流体本构方程与运动方程,建立有关力学平衡方程,并运用代数方程的根式解理论对圆管层流流动时的非线性方程进行求解,可直接求得Bingham流体圆管层流压降及速度流核区半径的解析解,进一步可求得圆管层流速度解析解;Bingham流体圆管层流速度的直接影响因素为流量、塑性粘度和屈服值,研究发现速度流核宽度与屈服值成正比,与流量及塑性粘度成反比,且流核的宽度越大,流核区的速度越小.     

粉末注射成形坯孔隙度的分形模型

粉末注射成形坯是一种具有分形特性的典型的多孔介质,借助于多孔介质孔隙结构的分形理论,对粉末注射成形坯孔隙率的分形模型进行推导。首先分析了粉末注射成形坯孔隙结构的双重分形特性,介绍了粉末注射成形坯孔隙分布分形维数和孔隙迂曲分形维数,然后推导出粉末注射成形坯孔隙度的分形模型。     

An existence result for a coupled system modeling a fully equivalent global pressure formulation for immiscible compressible two-phase flow in porous media

A new model describing immiscible, compressible two-phase flow, such as water-gas, through heterogeneous porous media is considered. The main feature of this model is the introduction of a new global pressure and the full equivalence to the original equations. The resulting equations are written in a fractional flow formulation and lead to a coupled system which consists of a nonlinear parabolic equation (the global pressure equation) and a nonlinear diffusion-convection one (the saturation equation). Under some realistic assumptions on the data, we show an existence result with the help of appropriate regularizations and a time discretization. We use suitable test functions to get a priori estimates. In order to pass to the limit in nonlinear terms, we also obtain compactness results which are nontrivial due to the degeneracy of the system.     

Steady flow in a deformable porous medium

A nonlinear model for a steady flow in a deformable porous medium is considered. The flow is governed by the poroelasticity system consisting of an elasticity equation for the displacement of the porous medium and Darcy's equation for the pressure in the fluid. This poroelasticity system is nonlinear when the permeability in Darcy's equation is assumed to depend on the dilatation of the porous medium. Existence and uniqueness of a weak solution of this poroelasticity system is established under rather weak assumptions on the regularity of the data. Convergence of a finite element approximation is proved and verified through numerical experiments. Copyright © 2013 John Wiley & Sons, Ltd.     

Coupled lattice Boltzmann method for simulating electrokinetic flows: A localized scheme for the Nernst–Plank model

We present a coupled lattice Boltzmann method (LBM) to solve a set of model equations for electrokinetic flows in micro-/nano-channels. The model consists of the Poisson equation for the electrical potential, the Nernst–Planck equation for the ion concentration, and the Navier–Stokes equation for the flows of the electrolyte solution. In the proposed LBM, the electrochemical migration and the convection of the electrolyte solution contributing to the ion flux are incorporated into the collision operator, which maintains the locality of the algorithm inherent to the original LBM. Furthermore, the Neumann-type boundary condition at the solid/liquid interface is then correctly imposed. In order to validate the present LBM, we consider an electro-osmotic flow in a slit between two charged infinite parallel plates, and the results of LBM computation are compared to the analytical solutions. Good agreement is obtained in the parameter range considered herein, including the case in which the nonlinearity of the Poisson equation due to the large potential variation manifests itself. We also apply the method to a two-dimensional problem of a finite-length microchannel with an entry and an exit. The steady state, as well as the transient behavior, of the electro-osmotic flow induced in the microchannel is investigated. It is shown that, although no external pressure difference is imposed, the presence of the entry and exit results in the occurrence of the local pressure gradient that causes a flow resistance reducing the magnitude of the electro-osmotic flow.     

A computational model of a transcritical R744 ejector based on a homogeneous real fluid approach

A mathematical model of the compressible transonic single- and two-phase flow of a real fluid is discussed in this paper. The model was originally developed to simulate a refrigerant flow through a heat pump ejector. In the proposed approach, a temperature-based energy equation is replaced with an enthalpy-based formulation, in which the specific enthalpy, instead of the temperature, is an independent variable. A thermodynamic and mechanical equilibrium between gaseous and liquid phases is assumed for the two-phase flow. Consequently, real fluid properties, such as the density, the dynamic viscosity and the diffusion coefficient, are defined as functions of the pressure and the specific enthalpy. The energy equation formulation is implemented in commercial CFD software using subroutines that were developed in-house. The formulations was tested extensively for a single-phase flow of the R141b refrigerant, and for a two-phase flow of the R744 fluid (carbon dioxide) that occurred in a 3-D model of the ejector motive nozzle. In the model validation procedure, a satisfactory comparison between the experimental and computational results of the primary and secondary mass flow rates was obtained for both flow regimes. In addition, in the case of the R744 flow, the pressure distribution along the centre line of the ejector was accurately predicted as well. Furthermore, the results also shows that geometry modelling and measurement accuracy play an important in the final numerical results. As a result of the reasonable computational times, this method can be effectively used for the design of ejectors and also in geometric optimisation computations.     

Precise integration method for natural frequencies and mode shapes of ocean risers with elastic boundary conditions

The effects of the elastic constraints simulating damaged and undamaged boundaries on the natural frequencies and mode shapes of ocean risers with a variable axial tension are investigated using the precise integration method (PIM). The classical high-order variable-coefficient partial-differential governing equation of the free vibration of ocean risers is reduced to a set of first-order ordinary differential equations and efficiently solved by the PIM. The main advantages of the PIM are that the numerical results can be calculated with high accuracy even when total element number n, number of iterations N and Taylor expansion terms r are small. Moreover, the computing time is quite short. Various boundary conditions are modeled as linear elastic constraints using a pair of translational and torsional springs, and four types of boundary damage coefficients are proposed to investigate the effects of a damaged boundary on the natural frequencies. The results for specific boundary conditions show agreement with those reported in the literature, and the calculation errors are very small in comparison with the analytical solution. Overall, the methodology of PIM is applicable for the investigation of the natural frequencies and mode shapes of ocean risers with a variable axial tension and cross-section and various boundary conditions.     

凝析气藏段塞流动的初—边值问题及褶积理论的应用

推出凝析油－气两相渗流方程,引入凝析油－气两相拟压力－时间函数,使渗流方程线性化,针对低渗凝析气藏的钻杆测试问题,在研究凝析油－气的段塞流动特征的基础上,建立了描述凝析气段塞流动的初－边值问题,求出这一问题解,并把褶积理论引入凝析气段塞流动资料分析模型中,建立了短时段塞流动资料解释的新理论与新方法,作为理论的应用文中给出了分析实例。     

Study of the contraction flow using grid generation technique

We study a model procedure to solve the incompressible Navier-Stokes equations on the flow inside contraction geometry. The governing equations are expressed in the primitive variable formulation. A rectangular computational plane is arises by elliptic grid generation technique. The numerical solution is based on a technique of automatic numerical generation of a curvilinear coordinate system. By transformed the governing equation into computational plane. The time dependent momentum equations are solved explicitly for the velocity field using the explicit marching procedure, the continuity equation is applied at each grid point in the solution of pressure equation, while the successive over relaxation (SOR) method is used for the Neumann problem for pressure. We will apply the technique on several irregular-shape.