Numerical Modeling of Interior Ballistics Processes in Light Gas Guns

An improved procedure for modeling the operation of a light-gas gun is proposed. The motion of working bodies in both the firing chamber and the light-gas chamber is studied within the framework of the mechanics of heterogeneous media. The problem of barrel heating taking into account its melting and removal of thermal ablation products into the medium inside the bore is solved in a coupled formulation. Heat and mass transfer and friction on the barrel surface are calculated using empirical dependences. The deformable piston is considered compressible and elastoviscoplastic. Allowance is made for the presence of a clearance between the lateral surface of the piston and the barrel bore walls and the associated gas flow between the firing and the light-gas chamber. Calculation results are given.     

Optimum Gas Journal Bearing with the Flexible Surface

The optimum shape of the cylindrical internal deformable (flexible) surface of a gas journal bearing of infinite length is designed. The variational problem of determining the clearance shape giving the maximum bearing load capacity is formulated and solved.     

Effect of Insoluble Surfactants on Drainage and Rupture of a Film between Drops Interacting under a Constant Force

The deformation, drainage, and rupture of an axisymmetrical film between colliding drops in the presence of insoluble surfactants under the influence of van der Waals forces is studied numerically at small capillary and Reynolds numbers and small surfactant concentrations. Constant-force collisions of Newtonian drops in another Newtonian fluid are considered. The mathematical model is based on the lubrication equations in the gap between drops and the creeping flow approximation of Navier–Stokes equations in the drops, coupled with velocity and stress boundary conditions at the interfaces. A nonuniform surfactant concentration on the interfaces, governed by a convection–diffusion equation, leads to a gradient of the interfacial tension which in turn leads to additional tangential stress on the interfaces (Marangoni effects). The mathematical problem is solved by a finite-difference method on a nonuniform mesh at the interfaces and a boundary-integral method in the drops. The whole range of the dispersed to continuous-phase viscosity ratios is investigated for a range of values of the dimensionless surfactant concentration, Peclét number, and dimensionless Hamaker constant (covering both “nose” and “rim” rupture). In the limit of the large Peclét number and the small dimensionless Hamaker constant (characteristic of drops in the millimeter size range) a fair approximation to the results is provided by a simple expression for the critical surfactant concentration, drainage being virtually uninfluenced by the surfactant for concentrations below the critical surfactant concentration and corresponding to that for immobile interfaces for concentrations above it.     

A coupled immersed boundary‐lattice Boltzmann method with smoothed point interpolation method for fluid‐structure interaction problems

The immersed boundary‐lattice Boltzmann method has been verified to be an effective tool for fluid‐structure interaction simulation associated with thin and flexible bodies. The newly developed smoothed point interpolation method (S‐PIM) can handle the largely deformable solids owing to its softened model stiffness and insensitivity to mesh distortion. In this work, a novel coupled method has been proposed by combining the immersed boundary‐lattice Boltzmann method with the S‐PIM for fluid‐structure interaction problems with large‐displacement solids. The proposed method preserves the simplicity of the lattice Boltzmann method for fluid solvers, utilizes the S‐PIM to establish the realistic constitutive laws for nonlinear solids, and avoids mesh regeneration based on the frame of the immersed boundary method. Both two‐ and three‐dimensional numerical examples have been carried out to validate the accuracy, convergence, and stability of the proposed method in consideration of comparative results with referenced solutions.     

可变形多孔介质渗透系数的测定方法

在Ｂｉｏｔ理论基础上给出可变形多孔介质耦合渗流基本方程；求出小试件一维定常耦合渗流问题的解答；表明在一维流固耦合情况下试件内部压力梯度有明显的非均匀性。因此通过实验确定可变形多孔介质渗透系数在数学上可归结为微分方程的反问题，传统的测试渗透系数的方法需要改进。介绍了可变形多孔介质渗透系数的测试原理和测试方法。对粒状多孔材料实验的结果表明，传统实验方法得到的渗透系数误差较大。     

A study of gridless method with dynamic clouds of points for solving unsteady CFD problems in aerodynamics

When solving unsteady computational fluid dynamics problems in aerodynamics with a gridless method, a cloud of points is usually required to be regenerated due to its accommodation to moving boundaries. In order to handle this problem conveniently, a fast dynamic cloud method based on Delaunay graph mapping strategy is proposed in this paper. A dynamic cloud method makes use of algebraic mapping principles and therefore points can be accurately redistributed in the flow field without any iteration. In this way, the structure of the gridless clouds is not necessarily changed so that the clouds regeneration can be avoided successfully. The spatial derivatives of the mathematical modeling of the flow are directly determined by using weighted least‐squares method in each cloud of points, and then numerical fluxes can be obtained. A dual time‐stepping method is further implemented to advance the two‐dimensional Euler equations in arbitrary Lagarangian–Eulerian formulation in time. Finally, unsteady transonic flows over two different oscillating airfoils are simulated with the above method and results obtained are in good agreement with the experimental data. Copyright © 2009 John Wiley & Sons, Ltd.     

A local domain‐free discretization method for simulation of incompressible flows over moving bodies

This paper presents a local domain‐free discretization (DFD) method for the simulation of unsteady flows over moving bodies governed by the incompressible Navier–Stokes equations. The discretization strategy of DFD is that the discrete form of partial differential equations at an interior point may involve some points outside the solution domain. All the mesh points are classified as interior points, exterior dependent points and exterior independent points. The functional values at the exterior dependent points are updated at each time step by the approximate form of solution near the boundary. When the body is moving, only the status of points is changed and the mesh can stay fixed. The issue of ‘freshly cleared nodes/cells’ encountered in usual sharp interface methods does not pose any particular difficulty in the presented method. The Galerkin finite‐element approximation is used for spatial discretization, and the discrete equations are integrated in time via a dual‐time‐stepping scheme based on artificial compressibility. In order to validate the present method for moving‐boundary flow problems, two groups of flow phenomena have been simulated: (1) flows over a fixed circular cylinder, a harmonic in‐line oscillating cylinder in fluid at rest and a transversely oscillating cylinder in uniform flow; (2) flows over a pure pitching airfoil, a heaving–pitching airfoil and a deforming airfoil. The predictions show good agreement with the published numerical results or experimental data. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical experiments with several variant WENO schemes for the Euler equations

Numerical experiments with several variants of the original weighted essentially non‐oscillatory (WENO) schemes (J. Comput. Phys. 1996; 126 :202–228) including anti‐diffusive flux corrections, the mapped WENO scheme, and modified smoothness indicator are tested for the Euler equations. The TVD Runge–Kutta explicit time‐integrating scheme is adopted for unsteady flow computations and lower–upper symmetric‐Gauss–Seidel (LU‐SGS) implicit method is employed for the computation of steady‐state solutions. A numerical flux of the variant WENO scheme in flux limiter form is presented, which consists of first‐order and high‐order fluxes and allows for a more flexible choice of low‐order schemes. Computations of unsteady oblique shock wave diffraction over a wedge and steady transonic flows over NACA 0012 and RAE 2822 airfoils are presented to test and compare the methods. Various aspects of the variant WENO methods including contact discontinuity sharpening and steady‐state convergence rate are examined. By using the WENO scheme with anti‐diffusive flux corrections, the present solutions indicate that good convergence rate can be achieved and high‐order accuracy is maintained and contact discontinuities are sharpened markedly as compared with the original WENO schemes on the same meshes. Copyright © 2008 John Wiley & Sons, Ltd.     

Study on the flow resistance of the dispersion system of deformable preformed particle gel in porous media using LBM-DEM-IMB method

After being injected into the porous media, the dispersion system of preformed particle gel (PPG) tends to enter high permeability regions and block water channeling passages, which forces the subsequent water to turn to the low permeability regions and thus increases sweep efficiency and enhances oil recovery. However, it is still unclear about the influence factors and the mechanisms how PPG increases water flow resistance, which limits the application of PPG in more oilfields. Therefore, the paper combines the lattice Boltzmann method (LBM), the discrete element method (DEM) and the improved immersed moving boundary (IMB) method to simulate the migration of deformable PPG in porous media. On the basis, the paper quantitatively analyzes the variation law of displacement pressure across the porous media and discusses the influence factors such as the PPG diameter, elasticity modulus and the number concentration. Results indicate that, because of the friction and retention of PPG in pore-throat, the displacement pressure across the porous media during PPG flooding is much higher than that during water flooding. In other words, the existence of PPG increases the flow resistance of injected water. Besides, the displacement pressure is always fluctuant resulting from the continuous process of PPG migration, retention, deformation and remigration. Influence factor analysis shows that the incremental value and fluctuation degree of flow resistance increase with the PPG diameter, elasticity modulus and the number concentration. The study not only provides useful reference for future PPG flooding, but also benefits the development of deformable particle flow theory.