Simulation of viscous water column collapse using adapting hierarchical grids

An adaptive hierarchical grid‐based method for predicting complex free surface flows is used to simulate collapse of a water column. Adapting quadtree grids are combined with a high‐resolution interface‐capturing approach and pressure‐based coupling of the Navier–Stokes equations. The Navier–Stokes flow solution scheme is verified for simulation of flow in a lid‐driven cavity at Re=1000. Two approaches to the coupling of the Navier–Stokes equations are investigated as are alternative face velocity and hanging node interpolations. Collapse of a water column as well as collapse of a water column and its subsequent interaction with an obstacle are simulated. The calculations are made on uniform and adapting quadtree grids, and the accuracy of the quadtree calculations is shown to be the same as those made on the equivalent uniform grids. Results are in excellent agreement with experimental and other numerical data. A sharp interface is maintained at the free surface. The new adapting quadtree‐based method achieves a considerable saving in the size of the computational grid and CPU time in comparison with calculations made on equivalent uniform grids. Copyright © 2005 John Wiley & Sons, Ltd.     

A subdomain boundary element method for high‐Reynolds laminar flow using stream function‐vorticity formulation

The paper presents a new formulation of the integral boundary element method (BEM) using subdomain technique. A continuous approximation of the function and the function derivative in the direction normal to the boundary element (further ‘normal flux’) is introduced for solving the general form of a parabolic diffusion‐convective equation. Double nodes for normal flux approximation are used. The gradient continuity is required at the interior subdomain corners where compatibility and equilibrium interface conditions are prescribed. The obtained system matrix with more equations than unknowns is solved using the fast iterative linear least squares based solver. The robustness and stability of the developed formulation is shown on the cases of a backward‐facing step flow and a square‐driven cavity flow up to the Reynolds number value 50 000. Copyright © 2004 John Wiley & Sons, Ltd.     

叶栅全三维粘性反问题的数值解

本文发展了一种解叶栅全三维粘性反问题的新的数值方法.基于非正交曲线坐标与相应的非正交速度分量下完全守恒型的Navier－Stokes方程,全三维反问题规定叶片表面的无量纲压力分布反求叶型。计算中叶片表面的边界条件采用一种特殊的方式来处理，即一方面强加给定的压力分布条件,另方面叶面的几何位置在迭代过程中又是可移动的，其移动速度将与Navier—Stokes方程在当地的解联系起来，从而形成一种解定常问题的新的不定常过程.试算证明了本文方法的可行性。     

Influence of Relative Permeability on the Stability Characteristics of Immiscible Flow in Porous Media

Relative permeability functions for immiscible displacements in porous media show a wide range of profiles. Although, this behavior is well known, its impact on the stability of the displacement process is unexplored. Our analysis clearly demonstrates for the first time that the viscous instability characteristics of two-phase flows are governed not only by their end point values, but are strongly dependent on the actual profile of relative permeability functions. Linear stability analysis predicts the capacity of the flow to develop large scale fingers which can result in substantial bypassing of the resident fluid. It is observed that relative permeability functions attributed to drainage processes yield a more unstable displacement as compared to functions related to imbibition processes. Moreover, instability is observed to increase for those relative permeability functions which result from increased wettability of the wetting fluid. High accuracy numerical simulations show agreement with these predictions and demonstrate how large amplitude viscous fingers result in significant bypassing for certain relative permeability functions. In the nonlinear regime, the finger amplitude grows at a rate ∝ t^1/2 initially, drops to t^1/4 at a later time and finally grows ∝ t. The basic mechanisms of finger interaction, however, are not substantially influenced by relative permeability functions.     

Nonlinear Modeling Study of Aerodynamic Characteristics of an X38-like Vehicle at Strong Viscous Interaction Regions

Strong viscous interaction and multiple flow regimes exist when vehicles fly at high altitude and high Mach number conditions. The Navier–Stokes(NS) solver is no longer applicable in the above situation. Instead, the direct simulation Monte Carlo (DSMC) method or Boltzmann model equation solvers are usually needed. However, they are computationally more expensive than the NS solver. Therefore, it is of great engineering value to establish the aerodynamic prediction model of vehicles at high altitude and high Mach number conditions. In this paper, the hypersonic aerodynamic characteristics of an X38-like vehicle in typical conditions from 70 km to 110 km are simulated using the unified gas kinetic scheme (UGKS), which is applicable for all flow regimes. The contributions of pressure and viscous stress on the force coefficients are analyzed. The viscous interaction parameters, Mach number, and angle of attack are used as independent variables, and the difference between the force coefficients calculated by UGKS and the Euler solver is used as a dependent variable to establish a nonlinear viscous interaction model between them in the range of 70–110 km. The evaluation of the model is completed using the correlation coefficient and the relative orthogonal distance. The conventional viscous interaction effect and rarefied effect are both taken into account in the model. The model can be used to quickly obtain the hypersonic aerodynamic characteristics of X38-like vehicle in a wide range, which is meaningful for engineering design.     

A moving particle semi‐implicit method for free surface flow: Improvement in inter‐particle force stabilization and consistency restoring

The moving particle semi‐implicit (MPS) method has been widely applied in free surface flows. However, the implementation of MPS remains limited because of compressive instability occurred when the particles are under compressive stress states. This study proposed an inter‐particle force stabilization and consistency restoring MPS (IFS‐CR‐MPS) method to overcome this numerical instability. For inter‐particle force stabilization, a hyperbolic‐shaped quintic kernel function is developed with a non‐negative and smooth second order derivative to satisfy the stability criterion under compressive stress state. Then, a contrastive study is conducted on the contradiction between the common understanding of the conventional MPS hyperbolic‐shaped kernel function and its performance. The result shows that the conventional MPS hyperbolic‐shaped kernel function can easily cause violent repulsive inter‐particle force and then lead to the compressive instability. Therefore, the first order derivative of the modified hyperbolic‐shaped quintic kernel function is recommended as the form of the contribution of the neighbor particles to achieve a more stable inter‐particle repulsive force. For consistency restoring, the Taylor series expansion and the hyperbolic‐shaped quintic kernel are combined to improve the accuracy of the viscosity and pressure calculation. The IFS‐CR‐MPS algorithm is subsequently verified by the inviscid hydrostatic pressure, jet impacting, and viscous droplet impacting problems. These results can be used for choosing kernel function and the contribution of neighbor particles in particle methods. Copyright © 2016 John Wiley & Sons, Ltd.     

Analytical and numerical solution for viscous dissipation effect on the heat transfer in a deformable channel

In this paper, shooting method and homotopy perturbation technique are applied for the flow analysis of temporal energy transport in a deformation channel with isothermal walls. An incompressible viscous fluid fills the space inside the channel. Analytical and numerical solutions are developed for the momentum and energy equations. The viscous dissipation effects are taken into account. Graphs for pertinent flow parameters are sketched and discussed. Comparison between the analytical and numerical solutions indicates an excellent agreement. It is noticed that behaviors of Prandtl and Eckert numbers on the temperature are qualitatively similar. Copyright © 2011 John Wiley & Sons, Ltd.     

Mixed convection unsteady stagnation point flow over a stretching sheet with heat transfer in the presence of variable free stream

An analysis is performed for the unsteady mixed convection flow of an incompressible viscous fluid about a stagnation point on a stretching sheet in the presence of a variable free stream. The equations of motion and energy are transformed into the ordinary differential equations by using similarity transformations. Homotopy analysis method is used for the solution of the governing problem. The results have been discussed by plots. The present values of the function are shown very close to the previous limiting solutions. Copyright © 2011 John Wiley & Sons, Ltd.     

A vectorial algorithm with finite element method for prediction of powder segregation in metal injection molding

The design of the mold and the choice of the injection parameters for metal injection molding (MIM) is required to maintain homogeneity of the filled mixture. However, powder segregation is unavoidable in MIM because of the significant difference in densities of the metallic powder and the polymer binder. To achieve an effective prediction of segregation effect, a biphasic model based on mixture theory is employed. The viscous behaviors of each phase and the interaction coefficient between the flows of the two phases should be determined. The solution of two coupled Navier–Stokes equations results in a tremendous computation effort. The previous development of an explicit algorithm makes the biphasic simulation much faster than that of the classic methods. However, it is strongly desired to reduce or even eliminate the numerous global solutions for pressure fields at each time step. Hence, a new vectorial algorithm is proposed and developed to perform the simulation only by vectorial operations. It provides the anticipated efficiency in the simulation of biphasic modeling, and the advantage to use the classic elements of equal‐order interpolations. Some results produced by the two algorithms are compared with the experimental values to validate the new vectorial algorithm. Copyright © 2012 John Wiley & Sons, Ltd.