Numerical simulation two phase flows of casting filling process using SOLA particle level set method

Mould filling process is a typical gas–liquid metal two phase flow phenomenon. Numerical simulation of the two phase flows of mould filling process can be used to properly predicate the back pressure effect, the gas entrapment defects, and better understand the complex motions of the gas phase and the liquid phase. In this paper, a novel sharp interface incompressible two phase numerical model for mould filling process is presented. A simple ghost fluid method like discretization method and a density evaluation method at face centers of finite difference staggered grid are proposed to overcome the difficulties when solving two phase Navier–Stokes equations with large-density ratio and large-viscosity ratio. A new mass conservation particle level set method is developed to capture the gas–liquid metal phase interface. The classical pressure-correction based SOLA algorithm is modified to solve the two phase Navier–Stokes equations. Two numerical tests including the Zalesak disk problem and the broken dam problem are used to demonstrate the accuracy of the present method. The numerical method is then adopted to simulate three mould filling examples including two high speed CCD camera imaging water filling experiments and an in situ X-ray imaging experiment of pure aluminum filling. The simulation results are in good agreement with the experiments.     

On a singularly perturbed mixed problem for a linear parabolic equation with nonlinear boundary conditions

A mixed problem for a linear parabolic equation with a small parameter multiplying the time derivative and with nonlinear boundary conditions is solved. Such boundary conditions arise in some heat and mass transfer problems, for example, in cooling thin spherical gas-filled shells or in the case of a gas filling such shells with gas-permeable walls.     

实际气体的幂律状态方程

为克服Onnes(昂内斯)气体状态方程参数多的缺点, 提出了描述实际气体的幂律状态方程.该模型仅包含两个参数,其中幂律函数的阶数可以为任意实数,刻画了实际气体偏离理想气体的程度.满足幂律状态方程的实际气体称为幂律气体.应用于描述氮气（N₂）和四氟甲烷（CF₄）两种实际气体的研究表明,与Onnes气体状态方程相比,幂律气体状态方程可以用较少的参数,准确地描述气体状态方程中压强和体积的幂律关系.此外,温度越低,幂律函数的阶数越小,反映了气体的实际状态越偏离理想气体.     

一类发烟装液弹装填率的数学模型

装液弹的装填率对弹的安全可靠性极其重要 .如果装液弹的装填率不合理 ,不是造成弹腔容积的浪费 ,就会造成弹内压过大 ,易于破坏密封性能引起渗漏甚至破裂而造成损失和危险 .本文通过对一类发烟装液炮弹弹内压力分析 ,得出了弹内压力计算公式 .在分析压力诸因素时 ,我们运用了固体热膨胀理论 ;考虑了液体的压缩性 .这是本文区别于以往压力计算的两个特点 .通过这类发烟装液炮弹内压压力曲线的分析 ,导出了装液 (炮 )弹弹腔空隙率合理选定最后的数学模型 (PLT方程——可由此确定合理的装填率 ) ,由此所计算的空隙率数据与国外文献值相符 ,并引入了“预极限温度”的概念 .可以相信 ,所谓“预极限温度”,将是装液 (炮 )弹 (或任何装液容器 )设计者必须认真考虑的问题     

An explicit staggered-grid method for numerical simulation of large-scale natural gas pipeline networks

We present an explicit second order staggered finite difference (FD) discretization scheme for forward simulation of natural gas transport in pipeline networks. By construction, this discretization approach guarantees that the conservation of mass condition is satisfied exactly. The mathematical model is formulated in terms of density, pressure, and mass flux variables, and as a result permits the use of a general equation of state to define the relation between the gas density and pressure for a given temperature. In a single pipe, the model represents the dynamics of the density by propagation of a non-linear wave according to a variable wave speed. We derive compatibility conditions for linking domain boundary values to enable efficient, explicit simulation of gas flows propagating through a network with pressure changes created by gas compressors. We compare our staggered grid method with an explicit operator splitting method and a lumped element scheme, and perform numerical experiments to validate the convergence order of the new discretization approach. In addition, we perform several computations to investigate the influence of non-ideal equation of state models and temperature effects on pipeline simulations with boundary conditions on various time and space scales.     

Simulation of eye deformation in the measurement of intraocular pressure

The procedure of measuring the intraocular pressure by an optical analyzer is numerically simulated. The cornea and the sclera are considered as axisymmetrically deformable shells of revolution with fixed boundaries; the space between these shells is filled with incompressible fluid. Nonlinear shell theory is used to describe the stressed and strained state of the cornea and sclera. The optical system is calculated from the viewpoint of the geometrical optics. Dependences between the pressure in the air jet and the area of the surface reflecting the light into a photodetector are obtained. The shapes of the regions on the cornea surface are found from which the reflected light falls on the photodetector. First, the light is reflected from the center of the cornea, but then, as the cornea deforms, the light is reflected from its periphery. The numerical results make it possible to better interpret the measurement data.     

Dynamics of a rotor partially filled with a viscous incompressible fluid

A rotor partially filled with a viscous incompressible fluid is modeled as planar system. Its structural part, i. e. the rotor, is assumed to be rigid, circular, elastically supported and running with a prescribed time-dependent angular velocity. Both parts, structure and fluid, interact via the no-slip condition and the pressure. The point of departure for the mathematical formulation of the fluid filling is the Navier-Stokes equation, which is complemented by an additional equation for the evolution of its free inner boundary. Further, rotor and fluid are subjected to volume forces, namely gravitation. Trial functions are chosen for the fluid velocity field, the pressure field and the moving boundary, which fulfill the incompressibility constraint as well as the boundary conditions. Inserting these trial functions into the partial differential equations of the fluid motion, and applying the method of weighted residuals yields equations with time derivatives only. Finally, in combination with the rotor equations, a nonlinear system of 12 differential-algebraic equations results, which sufficiently describes solutions near the circular symmetric state and which may indicate the loss of its stability. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

充气薄膜球囊测压法的数学方程

把流动坐标原点设在球囊中心,探讨了用充气薄膜球囊测量流体压力的数学方程,提出了应优选圆球囊测压,并分析了测压误差及误差处理方法,绘出了仿真曲线与实验曲线,通过实验验证了仿真的正确性.     

周向加肋非圆柱壳谐振分析的一个新矩阵方法

基于一类柱壳谐振控制方程呈一阶常微分矩阵方程形式以及傅立叶级数展开,提出了一种新矩阵方法,求解两端简支具有环肋加强非圆柱壳在谐外压作用下的稳态响应．该方法和以往同类方法相比,有两个突出的优点:1) 矩阵微分方程的解采用齐次扩容精细积分法替代龙格-库塔法,提高了精度;其中传递矩阵能实现计算机精确计算．2) 环肋作用力借助Dirac-δ函数和三角级数逼近可以解析求出;除法向作用力外,还考虑了切向作用力．通过数值计算,还研究了外激励频率对壳体位移和应力的影响规律．对比有限元分析与其它方法的计算结果,表明了该方法的正确性和有效性．     

Analysis of the Stress State of the Surface Layer of a Polymeric Mass upon Filling of Bulky Compression Molds

The stress state of the surface layer of a polymeric mass during filling of bulky compression molds is analyzed. It is shown that, at particular rheological characteristics of the mass, temperature, and filling rates, cracking of the surface layer occurs, which leads to defects in the finished products. A physical analysis of this process makes it possible to conclude that the cracks arise due to the normal stresses operating in the front region of the moving polymeric mass. It is found that, under certain flow conditions, areas with a pressure lower than the atmospheric one appear on the surface of the polymer. If the tensile stresses arising in these local regions are higher than the tensile strength of the mass, the continuity of the composition is broken in the direction determined by the greatest rate of the normal deformation. To confirm the reliability of the crack-formation mechanism proposed, the distribution of the pressure and normal stresses over the free surface is calculated based on a numerical method. These calculations show that, by comparing the stress level achieved in the front region with the tensile-strength characteristics of the polymeric composition, it is possible to predict, with a sufficient accuracy, the possibility of crack formation in the surface layer of such a mass under given flow conditions and thus to solve the question on flawless manufacturing of products.     

Analysis of equations of state for a real gas based on the notion of division of strain into free and elastic strains

It is established that, in a general case, the equations of state of a real gas are identical to the constitutive equations of cubic strain of a deformable material. Based on this, we use the fundamental notion of deformation mechanics of the division of strain into free and elastic strain. As a result, we reduce the equation of cubic strain to the form of a generalization of the ideal gas law. We establish specific features of the gaseous state as compared with the solid and liquid states. In particular, the continuum physical meaning of an ideal gas and the universal gas constant is revealed. We construct two-parameter expressions for the compressibility coefficient for an arbitrary state of aggregation, which reveal the physical meaning of the corresponding experimental data and agree with the notions of the molecularkinetic theory of gases.     

环向和纵向加肋非均匀大变形圆柱壳的自由振动*

环向和纵向加肋非均匀圆柱壳在航空,宇航等工业广泛运用,本文使用阶梯折算法得到了环向和纵向加肋非均匀大变形圆柱壳在任意边界条件下自由振动的一般解.问题最后归结为求解一个超越代数方程,这个方程可以用一个具体的解析表达式表示出来.文中还给出了收敛性证明和算例.算例表明,利用本文的方法,可得到满意的结果.     

凝析油-气在存在裂缝的双渗介质中渗流的变产量数学模型及其应用

当凝析气藏压力低于露点压力时,有凝析油析出,致使井底附近阻力变大,影响凝析气的采出量,反凝析现象带来了一些特殊的经济问题.选择合理的开发方式,提高凝析气的采收率,成为了当前工作中迫切需要解决的难题.为了控制合理的井底流压常会引起产量的改变.主要针对凝析气在存在裂缝的双渗气藏中的渗流问题,建立气井变产量生产新模型并进行求解,同时选择相关合理参数绘制了试井理论图版,根据试井理论图版利用试井测试数据进行解释,可以得到相应动态参数为气田开发预案提供相应的理论依据.     

Heat and mass transfer in annular liquid jets: I. Formulation

Heat and mass transfer phenomena in annular liquid jets are analyzed at high Reynolds numbers by means of a model derived from the governing equations that takes into account the effects of surface tension and boundary conditions at the gas–liquid interfaces and the large differences between the thermal and mass diffusivities, densities, dynamic viscosities, and thermal conductivities between gases and liquids. The model clearly illustrates the stiffness in both space and time associated with the concentration, linear momentum and energy boundary layers, and the initial cooling of the gases enclosed by the jet when, starting from a steady state where gases are injected into the volume enclosed by the jet at a rate equal to the heat and mass absorption rates by the liquid, gas injection is stopped. It is shown that, owing to the non-linear integrodifferential coupling between the fluid dynamics and heat and mass transfer processes, the pressure of the gases enclosed by the jet may vary in either a monotonic or an oscillatory manner depending on the large number of non-dimensional parameters that govern the heat and mass transfer phenomena. For the underpressurized jets considered here, it is shown that thermal equilibrium is achieved at a much faster rate than that associated with mass transfer, double diffusive phenomena in the liquid may occur, and the mass and volume of the gases enclosed by the jet may increase or decrease as functions of time until a steady equilibrium condition is reached.     

Modelling of high-pressure gas transmission lines

Gas flow in a transmission line is described by a set of three coupled partial differential equations (PDE) expressing conservation of mass, momentum and energy; the gas properties are described by a non-ideal equation of state. A technique is introduced which reduces the energy equation into a single parameter in the mass equation without the assumption of isothermal or isentropic flow.The remaining set of PDEs is solved by two different techniques. An accurate but time-consuming technique consists of applying the method of characteristics, for which an improved representation of the friction term is presented. The second way consists of a finite-difference implementation with a second-order truncation error on an analogue computer. Both the physical assumptions and the numerical approximations are checked against data obtained from experiments in the main transport system of Gasunie.Guidelines on the analogue modelling of pipeline systems and the interpretation of simulation results conclude the paper.     

Transient modeling of non-isothermal, dispersed two-phase flow in natural gas pipelines

Unsteady-state or transient two-phase flow, caused by any change in rates, pressures or temperature at any location in a two-phase flow line, may last from a few seconds to several hours. In general, these changes are an order of magnitude longer than the transient encountered during single-phase flow. The primary reason for this phenomenon is that the velocity of wave propagation in a two-phase mixture is significantly slower. Interfacial transfer of mass, momentum and energy further complicate the problem. It is primarily due to the numerical difficulties anticipated in accurately modeling transient two-phase flow that the state of the art in this important area is restricted to a handful of studies with direct applicability to petroleum and gas engineering. A limited amount of information on the subject of two-phase transport phenomena is available in the petroleum engineering literature. Most of the publications for two-phase flow of gas assume that temperature is constant over the entire length of the pipeline.This study is the first effort to simulate the non-isothermal, one-dimensional, transient homogenous two-phase flow gas pipeline system using two-fluid conservation equations. The modified Peng–Robinson equation of state is used to calculate the vapor–liquid equilibrium in multi-component natural gas to find the vapor and liquid compressibility factors. Mass transfer between the gas and the liquid phases is treated rigorously through flash calculation, making the algorithm capable of handling retrograde condensation. The liquid droplets are assumed to be spheres of uniform size, evenly dispersed throughout the gas phase.The method of solution is the fully implicit finite difference method. This method is stable for gas pipeline simulations when using a large time step and therefore minimizes the computation time. The algorithm used to solve the non-linear finite difference thermo-fluid equations for two-phase flow through a pipe is based on the Newton–Raphson method.The results show that the liquid condensate holdup is a strong function of temperature, pressure, mass flow rate, and mixture composition. Also, the fully implicit method has advantages, such as the guaranteed stability for large time step, which is very useful for simulating long-term transients in natural gas pipeline systems.     

Free Cavitational Deceleration of a Circular Cylinder in a Liquid after Impact

The problem is considered about the vertical continuous impact and subsequent free deceleration of a circular cylinder semi-immersed in a liquid. The specificity of this problem is that, under certain conditions, some areas of low pressure near the body appear and the attached cavities are formed. The separation zones and the motion law of the cylinder are unknown in advance and have to be determined in solving the problem. The study of the problem is conducted by a direct asymptotic method effective for small spans of time. Some nonlinear problem with unilateral constraints is formulated that is solved together with the equation defining the law of motion of the cylinder. In the case when the space above the external free surface of a liquid is filled with a gas with low pressure (vacuum), an analytical solution of the problem is constructed. To determine the main hydrodynamic characteristics (the separation point and acceleration of the cylinder), we derive a system of transcendental equations with elementary functions. The solution of this system agrees well with the results obtained by the direct numerical method.     

A combination of energy method and spectral analysis for study of equations of gas motion

There have been extensive studies on the large time behavior of solutions to systems on gas motions, such as the Navier-Stokes equations and the Boltzmann equation. Recently, an approach is introduced by combining the energy method and the spectral analysis to the study of the optimal rates of convergence to the asymptotic profiles. In this paper, we will first illustrate this method by using some simple model and then we will present some recent results on the Navier-Stokes equations and the Boltzmann equation. Precisely, we prove the stability of the non-trivial steady state for the Navier-Stokes equations with potential forces and also obtain the optimal rate of convergence of solutions toward the steady state. The same issue was also studied for the Boltzmann equation in the presence of the general time-space dependent forces. It is expected that this approach can also be applied to other dissipative systems in fluid dynamics and kinetic models such as the model system of radiating gas and the Vlasov-Poisson-Boltzmann system.      

薄球壳在均布外压与温度耦合作用下的热屈曲研究

从张量方法推导出的轴对称薄球壳屈曲方程出发,推导出在均布外压与温度耦合作用下用位移表示的薄球壳热屈曲方程;应用虚功原理建立薄球壳屈曲最小势能泛函;进一步用Ritz(里兹)法分析了周边简支的半球壳的3种热屈曲问题.得到了: 1) 温度不超过屈曲临界温度值时,均布外压的临界载荷;2) 均布外压载荷为0时,屈曲临界温度值;3) 均布外压载荷不超过临界载荷时,屈曲临界温度值.     

扁柱面网壳的非线性动力学行为

用连续化法建立了正三角形网格的三向单层扁柱面网壳的非线性动力学方程和协调方程．在两对边简支条件下用分离变量函数法给出扁柱面网壳的横向位移．由协调方程求出张力,通过Galerkin作用得到了一个含二次、三次的非线性动力学微分方程．通过求Floquet指数讨论平衡点邻域的稳定性,用复变函数留数理论求出Melnikov函数,可得到该动力学系统发生混沌运动的临界条件．通过数值计算模拟和Poincaré映射也证明了混沌运动存在．