An immersed-boundary method for compressible viscous flows and its application in the gas-kinetic BGK scheme

An immersed-boundary (IB) method is proposed and applied in the gas-kinetic BGK scheme to simulate incompressible and compressible viscous flows with complex stationary and moving boundaries on stationary Cartesian grids. In this method the ghost-cell technique is used to satisfy the boundary condition on the immersed boundary. A novel idea, “local boundary determination”, is put forward to identify the ghost cells, each of which may have several different ghost-cell constructions corresponding to different boundary segments. Thus, the singular behavior of the ghost cell is eliminated. Furthermore, the so-called “fresh-cell” problem that occurs when implementing the IB method in a moving-boundary simulation is resolved by a simple temporal extrapolation. The method is first applied in the gas-kinetic BGK scheme to simulate the Taylor–Couette flow, wherein the second-order spatial accuracy of the method is validated and the “super-convergence” of the BGK scheme is observed. After that the flow between a circular cylinder and a square cylinder is used as a test case to showcase the advantage of this method in resolving the singularity problem. Then the supersonic flow around a stationary cylinder, the incompressible flow around an oscillating cylinder and the compressible flow around a moving airfoil are simulated to verify that this method can be used to simulate compressible flows and handle moving boundaries. These numerical tests demonstrate the good performance of the proposed immersed-boundary method for the study of incompressible/compressible flow problems with complex stationary/moving boundaries.     

Application of the multigrid approach for solving 3D Navier-Stokes equations on hexahedral grids using the discontinuous Galerkin method

The Galerkin method with discontinuous basis functions is adapted for solving the Euler and Navier-Stokes equations on unstructured hexahedral grids. A hybrid multigrid algorithm involving the finite element and grid stages is used as an iterative solution method. Numerical results of calculating the sphere inviscid flow, viscous flow in a bent pipe, and turbulent flow past a wing are presented. The numerical results and the computational cost are compared with those obtained using the finite volume method.     

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

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

An adaptive integration scheme for heat conduction in additive manufacturing

Solidification dynamics are important for determining final microstructure in additively manufactured parts. Recently, researchers have adopted semi-analytical approaches for predicting heat conduction effects at length and time scales not accessible to complex multi-physics numerical models. The present work focuses on improving a semi-analytical heat conduction model for additive manufacturing by designing an adaptive integration technique. The proposed scheme considers material properties, process conditions, and the inherent physical behavior of the transient heat conduction around both stationary and moving heat sources. Both the adaptive integration scheme and a technique for calculating only the points within the melt pools are described in detail. The full algorithm is then implemented and compared against a simple Riemann sum integration scheme for a variety of cases that highlight process and material variations relevant to additive manufacturing. The new scheme is shown to have significant improvements in computational efficiency, solution accuracy, and usability.     

圆柱形平头弹体对薄靶板的垂直冲击

本文提出了弹体冲击靶板时弹靶接触面的运动速度和接触面上所受法向应力的解析表达式.这些解析式是著名的Hopkins-Kolsky理论的推广.由于弹道极限速度在工程实际中的重要性,本文也给出了弹道极限速度的解析表达式.本文并证明了,作用在与靶板相接触的冲塞圆柱面上的沿板厚方向的剪应力,与板厚方向的坐标无关.     

Front tracking for two-phase flow in reservoir simulation by adaptive mesh

In this article, an algorithm for the numerical approximation of two-phase flow in porous media by adaptive mesh is presented. A convergent and conservative finite volume scheme for an elliptic equation is proposed, together with the finite difference schemes, upwind and MUSCL, for a hyperbolic equation on grids with local refinement. Hence, an IMPES method is applied in an adaptive composite grid to track the front of a moving solution. An object-oriented programmation technique is used. The computational results for different examples illustrate the efficiency of the proposed algorithm. © 1997 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 13: 673–697, 1997     

基于Roe格式的可压与不可压流的统一计算方法

摘要：以Navier-Stokes方程为基础，基于有限体积的时间推进的预处理技术．提出了一个可以用来求解可压与不可压流场的统一的计算方法，原始变量选用压力、速度与温度，通过矩阵变换与重构，使得对流项系数矩阵在可压与小可压条件下都不会奇异．将可压与不可压流场的计算方法统一起来。采用Roe格式计算对流通量，采用中心差分格式计算扩散通量．算例表明，该方法可以进行高Mach数、中等Mach数、低Mach数及不可压流场的计算。由于采用了Roe格式，该方法还可以捕获不连续流场的间断面。     

The inviscid nature of the asymmetry of the separating flow of a uniform stream around symmetric bodies

The “inviscid” nature of the asymmetry is demonstrated using the example of the separating unsteady flow of an ideal incompressible fluid around a cylinder which is expanding at a constant velocity, that is, a non-steady-state analogue of steady-state flow around a cone at an angle of attack. An asymmetric flow structure is realized for a symmetrical positioning of the points of separation of the vortex sheets. This is evidence of the secondary role of viscosity, which can manifest itself through an “inverse” effect on the position of the points of separation. New asymmetric solutions and processes by which they arise, which are different from the classical bifurcation of the symmetric solution, are found. Together with an investigation of stability, an analysis of the global pattern of “self-similar” streamlines is carried out in the selection of the “realizable” solutions. The global pattern must correspond to the scheme adopted when constructing the theoretical model.     

A numerical scheme for moving bottlenecks in traffic flow

We consider a stronglycoupled PDE-ODE systemthat describes themoving bottlenecks created by several buses on a road. The model consists of a scalar conservation law modeling the main traffic evolution and a series of ODEs accounting for the trajectories of the buses, which depend on the surrounding traffic density. The moving bottlenecks are expressed by inequality constraints on the flux.We generalize a conservative finite volume scheme for the constrained hyperbolic PDE using a reconstruction technique and a tracking algorithm for the ODEs. Numerical simulations illustrate the impact of the buses on traffic flow.     

A composite Chebyshev finite difference method for nonlinear optimal control problems

In this paper, a composite Chebyshev finite difference method is introduced and is successfully employed for solving nonlinear optimal control problems. The proposed method is an extension of the Chebyshev finite difference scheme. This method can be regarded as a non-uniform finite difference scheme and is based on a hybrid of block-pulse functions and Chebyshev polynomials using the well-known Chebyshev–Gauss–Lobatto points. The convergence of the method is established. The nice properties of hybrid functions are then used to convert the nonlinear optimal control problem into a nonlinear mathematical programming one that can be solved efficiently by a globally convergent algorithm. The validity and applicability of the proposed method are demonstrated through some numerical examples. The method is simple, easy to implement and yields very accurate results.     

Chaotic invasive weed optimization algorithm with application to parameter estimation of chaotic systems

This paper introduces a novel hybrid optimization algorithm by taking advantage of the stochastic properties of chaotic search and the invasive weed optimization (IWO) method. In order to deal with the weaknesses associated with the conventional method, the proposed chaotic invasive weed optimization (CIWO) algorithm is presented which incorporates the capabilities of chaotic search methods. The functionality of the proposed optimization algorithm is investigated through several benchmark multi-dimensional functions. Furthermore, an identification technique for chaotic systems based on the CIWO algorithm is outlined and validated by several examples. The results established upon the proposed scheme are also supplemented which demonstrate superior performance with respect to other conventional methods.     

A Hybrid WENO Method with Modified Ghost Fluid Method for Compressible Two-Medium Flow Problems

In this paper, we develop a simplified hybrid weighted essentially non-oscillatory (WENO) method combined with the modified ghost fluid method (MGFM) [31] to simulate the compressible two-medium flow problems. The MGFM can turn the two-medium flow problems into two single-medium cases by defining the ghost fluids state in terms of the predicted the interface state, which makes the material interface “invisible”. For the single medium flow case, we adapt between the linear upwind scheme and the WENO scheme automatically by identifying the regions of the extreme points for the reconstruction polynomial as same as the hybrid WENO scheme [55]. Instead of calculating their exact locations, we only need to know the regions of the extreme points based on the zero point existence theorem, which is simpler for implementation and saves computation time. Meanwhile, it still keeps the robustness and has high efficiency. Extensive numerical results for both one and two dimensional two-medium flow problems are performed to demonstrate the good performances of the proposed method.     

Time‐dependent algorithms for viscoelastic flow: Finite element/volume schemes

Hybrid finite volume/element methods are investigated within the context of transient viscoelastic flows. A finite volume algorithm is proposed for the hyperbolic constitutive equation, of Oldroyd‐form, whereas the continuity/momentum balance is accommodated through a Taylor‐Galerkin finite element method. Various finite volume combinations are considered to derive accurate and stable implementations. Consistency of formulation is key, embracing fluctuation distribution and median‐dual‐cell constructs, within a cell‐vertex discretisation on triangles. In addition, we investigate the effect of treating the time‐term in a finite element fashion, using mass‐matrix iteration instead of the standard finite volume mass‐lumping approach. We devise an accurate transient scheme that captures the analytical solution at short and long time, both in core flow and near shear boundaries. In this respect, some difficulties are highlighted. A new method emerges, with the Low Diffusion B (LDB, with or without mass‐matrix iteration) as the optimal choice. We progress to a complex flow application and demonstrate some provocative features due to the influence of true transient boundary conditions on evolutionary flow‐structure in a 4:1 start‐up rounded‐corner contraction problem. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Finite Volume Methods for Multi-Symplectic PDES

We investigate the application of a cell-vertex finite volume discretization to multi-symplectic PDEs. The investigated discretization reduces to the Preissman box scheme when used on a rectangular grid. Concerning arbitrary quadrilateral grids, we show that only methods with parallelogram-like finite volume cells lead to a multi-symplectic discretization; i.e., to a method that preserves a discrete conservation law of symplecticity. One of the advantages of finite volume methods is that they can be easily adjusted to variable meshes. But, although the implementation of moving mesh finite volume methods for multi-symplectic PDEs is rather straightforward, the restriction to parallelogram-like cells implies that only meshes moving with a constant speed are multi-symplectic. To overcome this restriction, we suggest the implementation of reversible moving mesh methods based on a semi-Lagrangian approach. Numerical experiments are presented for a one dimensional dispersive shallow-water system.     

Fault diagnosis for a class of nonlinear uncertain hybrid systems

This paper presents a fault diagnosis architecture for a class of hybrid systems with nonlinear uncertain time-driven dynamics, measurement noise, and autonomous and controlled mode transitions. The proposed approach features a hybrid estimator based on a modified hybrid automaton framework. The fault detection scheme employs a filtering approach that attenuates the effect of the measurement noise and allows tighter mode-dependent thresholds for the detection of both discrete and parametric faults while guaranteeing no false alarms due to modeling uncertainty and mode mismatches. Both the hybrid estimator and the fault detection scheme are linked with an autonomous guard events identification (AGEI) scheme that handles the effects of mode mismatches due to autonomous mode transitions and allows effective mode estimation. Finally, the fault isolation scheme anticipates which fault events may have occurred and dynamically employs the appropriate isolation estimators for isolating the fault by calculating suitable thresholds and estimating the parametric fault magnitude through adaptive approximation methods. Simulation results from a five-tank hybrid system illustrate the effectiveness of the proposed approach.     

Three-dimensional mixed finite element-finite volume approach for the solution of density-dependent flow in porous media

The density-dependent flow and transport problem in groundwater on three-dimensional triangulations is solved numerically by means of a mixed hybrid finite element scheme for the flow equation combined with a mixed hybrid finite element-finite volume (MHFE-FV) time-splitting-based technique for the transport equation. This procedure is analyzed and shown to be an effective tool in particular when the process is advection dominated or when density variations induce the formation of instabilities in the flow field. From a computational point of view, the most effective strategy turns out to be a combination of the MHFE and a spatially variable time-splitting technique in which the FV scheme is given by a second-order linear reconstruction based on the least-squares minimization and the Barth–Jespersen limiter. The recent saltpool problem introduced as a benchmark test for density-dependent solvers is used to verify the accuracy and reliability of this approach.     

Initial states iterative learning for three-dimensional ballistic endpoint control

In this paper, an initial states iterative learning control algorithm is proposed for control of the ballistic endpoint displacement in three-dimensional space, where the target is moving and the projectile experiences system uncertainties. The characteristics of the three-dimensional ballistic process are formulated and explored, and the learning algorithm is proposed in the spatial domain. The algorithm consists of two parts. First, the initial speed and angles are iteratively learned to make the projectile attain a fixed position. Second, the shooting time is learned to tune the arrival time of the projectile. Since the dimensions of the solution space are larger than that of the task space, three control manners, including shooting speed, shooting angle and their combination, are researched respectively. Through rigorously analyzed, it is proved that the algorithm is convergent and the multiple initial states can be adjusted simultaneously. Finally, an example of practical cannonball projection is presented to verify the effectiveness of the proposed algorithms.     

Numerical simulation and visualization of vortical structure transformation in the flow past a sphere at an increasing degree of stratification

Based on numerical simulation and visualization, the vortex structure of the flow past a sphere moving uniformly and horizontally in a linearly (density) stratified viscous fluid with an increasing degree of stratification (with the internal Froude number Fr decreasing from infinity to 0.005) at Re = 100 is analyzed in detail for the first time. The classification of the flow regimes is refined. The direct numerical simulation is based on the method of splitting with respect to physical factors (MERANZH) with an explicit hybrid finite-difference scheme, which is second-order accurate in space, monotone, and has a minimal numerical viscosity and dispersion.     

A numerical method for mass conservative coupling between fluid flow and solute transport

We present a new coupled discretization approach for species transport in an incompressible fluid. The Navier-Stokes equations for the flow are discretized by the divergence-free Scott-Vogelius element on barycentrically refined meshes guaranteeing LBB stability. The convection-diffusion equation for species transport is discretized by the Voronoi finite volume method. In accordance to the continuous setting, due to the exact integration of the normal component of the flow through the Voronoi surfaces, the species concentration fulfills discrete global and local maximum principles. Besides of the numerical scheme itself, we present important aspects of its implementation. Further, for the case of homogeneous Dirichlet boundary conditions, we give a convergence proof for the coupled scheme. We report results of the application of the scheme to the interpretation of limiting current measurements in an electrochemical flow cell with cylindrical shape.     

两相流体非线性渗流模型及其应用

基于三参数非线性渗流运动定律、质量守恒定律及椭圆渗流的概念,建立了低渗透介质中两相流体椭圆非线性渗流数学模型,运用有限差分法与外推法求得了其解,导出了两相流体椭圆非线性渗流条件下油井见水前后开发指标的计算公式,进行了实例分析。结果表明:非线性渗流对含水饱和度分布影响较大;非线性渗流使得水驱油推进速度比线性渗流的快,使油井见水时间提前,使得石油开发指标变差;非线性渗流使得同一时刻的压差比线性渗流的大,使石油开发难度加大。这为低渗油藏垂直裂缝井开发工程提供了科学依据。