MHD free convection flow over an isothermal vertical cone with temperature dependent viscosity

Free convection flow over an isothermal vertical cone immersed in a fluid with variable viscosity and MHD is studied in this paper. Using appropriate variables, the basic equations are transformed into the non-dimensional boundary-layer equations. These equations are then solved numerically using a very efficient implicit finite-difference method known as Crankl-Nicolson scheme. Detailed results for the velocity, temperature, skin friction, and heat transfer rates for a selection of parameter sets consisting of the viscosity parameter, magnetic field parameter, and Prandtl number are discussed. In order to validate our numerical results, the present results are compared with the available work in the literature and are found to be in an excellent agreement.     

一种减少壁热误差的自适应热粘性

使用单元中心型拉氏方法,研究一维球、柱坐标等径向对称流体力学方程组的数值格式和减少壁热误差的方法.简要分析壁热误差与差分格式修正方程的关系.通过对Riemann问题声波和HLL近似解的比较,提出减少壁热误差的一种自适应的热通量粘性.多项数值实验表明该方法可以获得令人满意的计算结果.     

Immersed boundary method for the simulation of 2D viscous flow based on vorticity–velocity formulations

A new immersed boundary method based on vorticity–velocity formulations for the simulation of 2D incompressible viscous flow is proposed in present paper. The velocity and vorticity are respectively divided into two parts: one is the velocity and vorticity without the influence of the immersed boundary, and the other is the corrected velocity and the corrected vorticity derived from the influence of the immersed boundary. The corrected velocity is obtained from the multi-direct forcing to ensure the well satisfaction of the no-slip boundary condition at the immersed boundary. The corrected vorticity is derived from the vorticity transport equation. The third-order Runge–Kutta for time stepping, the fourth-order finite difference scheme for spatial derivatives and the fourth-order discretized Poisson for solving velocity are applied in present flow solver. Three cases including decaying vortices, flow past a stationary circular cylinder and an in-line oscillating cylinder in a fluid at rest are conducted to validate the method proposed in this paper. And the results of the simulations show good agreements with previous numerical and experimental results. This indicates the validity and the accuracy of present immersed boundary method based on vorticity–velocity formulations.     

A differentially interpolated direct forcing immersed boundary method for predicting incompressible Navier–Stokes equations in time-varying complex geometries

A dispersion-relation-preserving dual-compact scheme developed in Cartesian grids is applied together with the immersed boundary method to solve the flow equations in irregular and time-varying domains. The artificial momentum forcing term applied at certain points in cells containing fluid and solid allows an imposition of velocity condition to account for the motion of solid body. We develop in this study a differential-based interpolation scheme which can be easily extended to three-dimensional simulation. The results simulated from the proposed immersed boundary method agree well with other numerical and experimental results for the chosen benchmark problems. The accuracy and fidelity of the IB flow solver developed to predict flows with irregular boundaries are therefore demonstrated.     

一种基于多矩的有限体积浸入边界法

基于多矩VSIAM3格式及浸入边界法,提出一套在复杂计算区域内求解不可压缩流动的数值格式.不可压N-S方程使用VSIAM3格式进行离散,引入浸入边界法处理复杂、移动边界,使用虚拟网格方法计算动量方程修正项,同时还考虑了对连续方程的修正.使用标准算例对数值模式进行验证.     

Method of adaptive artificial viscosity

A new finite-difference method for the numerical solution of gas dynamics equations is proposed. This method is a uniform monotonous finite-difference scheme of second-order approximation on time and space outside of domains of shock and compression waves. This method is based on inputting adaptive artificial viscosity (AAV) into gas dynamics equations. In this paper, this method is analyzed for 2D geometry. The testing computations of the movement of contact discontinuities and shock waves and the breakup of discontinuities are demonstrated.     

非结构网格下非定常流场双时间步长的加权ENO-强紧致杂交高分辨率格式

针对三维非定常、可压缩流场的Navier-Stokes方程组,本文提出一种新的双时间步长高精度快速迭代格式。该格式在时间上具有二阶精度,在空间离散上不低于三阶。在对流项与粘性项的处理上,本格式分别采用了加权ENO-强紧致格式与紧致四阶精度格式的思想。几个典型算例的实践表明:计算结果与相关实验数据比较吻合,初步表明了该算法可以在非结构网格下具有高效率与高分辨率的特征。     

A momentum exchange-based immersed boundary-lattice Boltzmann method for simulating incompressible viscous flows

A momentum exchange-based immersed boundary-lattice Boltzmann method is presented in this Letter for simulating incompressible viscous flows. This method combines the good features of the lattice Boltzmann method (LBM) and the immersed boundary method (IBM) by using two unrelated computational meshes, an Eulerian mesh for the flow domain and a Lagrangian mesh for the solid boundaries in the flow. In this method, the non-slip boundary condition is enforced by introducing a forcing term into the lattice Boltzmann equation (LBE). Unlike the conventional IBM using the penalty method with a user-defined parameter or the direct forcing scheme based on the Navier–Stokes (NS) equations, the forcing term is simply calculated by the momentum exchange of the boundary particle density distribution functions, which are interpolated by the Lagrangian polynomials from the underlying Eulerian mesh. Numerical examples show that the present method can provide very accurate numerical results.     

流体力学拉氏守恒滑移线算法设计

本文在Wilkins滑移线算法的基础上, 设计了从区网格边人工黏性对主点施加算法. 结合拉氏“相容性格式”, 通过定义“接触力”和“接触力做功”, 利用局部修正内能的方法设计了总能量守恒的流体力学拉氏格式. 该修正方法可保证对称性、守恒性; 提高数值模拟分辨率.     

Nested Cartesian grid method in incompressible viscous fluid flow

In this work, the local grid refinement procedure is focused by using a nested Cartesian grid formulation. The method is developed for simulating unsteady viscous incompressible flows with complex immersed boundaries. A finite-volume formulation based on globally second-order accurate central-difference schemes is adopted here in conjunction with a two-step fractional-step procedure. The key aspects that needed to be considered in developing such a nested grid solver are proper imposition of interface conditions on the nested-block boundaries, and accurate discretization of the governing equations in cells that are with block-interface as a control-surface. The interpolation procedure adopted in the study allows systematic development of a discretization scheme that preserves global second-order spatial accuracy of the underlying solver, and as a result high efficiency/accuracy nested grid discretization method is developed. Herein the proposed nested grid method has been widely tested through effective simulation of four different classes of unsteady incompressible viscous flows, thereby demonstrating its performance in the solution of various complex flow–structure interactions. The numerical examples include a lid-driven cavity flow and Pearson vortex problems, flow past a circular cylinder symmetrically installed in a channel, flow past an elliptic cylinder at an angle of attack, and flow past two tandem circular cylinders of unequal diameters. For the numerical simulations of flows past bluff bodies an immersed boundary (IB) method has been implemented in which the solid object is represented by a distributed body force in the Navier–Stokes equations. The main advantages of the implemented immersed boundary method are that the simulations could be performed on a regular Cartesian grid and applied to multiple nested-block (Cartesian) structured grids without any difficulty. Through the numerical experiments the strength of the solver in effectively/accurately simulating various complex flows past different forms of immersed boundaries is extensively demonstrated, in which the nested Cartesian grid method was suitably combined together with the fractional-step algorithm to speed up the solution procedure.     

圆筒形半透明介质内非稳态复合导热与辐射的研究

本文研究细长电加热体在圆筒形半透明介质中的温度响应，通过对控制非稳态导热与辐射复合传热过程的积分-微分方程直接进行数值求解，分析了热辐射对内部径向热流及温度变化的影响．模拟计算结果对热线法测量半透明介质导热系数的研究具有理论指导意义。     

A lattice Boltzmann method for immiscible multiphase flow simulations using the level set method

We consider the lattice Boltzmann method for immiscible multiphase flow simulations. Classical lattice Boltzmann methods for this problem, e.g. the colour gradient method or the free energy approach, can only be applied when density and viscosity ratios are small. Moreover, they use additional fields defined on the whole domain to describe the different phases and model phase separation by special interactions at each node. In contrast, our approach simulates the flow using a single field and separates the fluid phases by a free moving interface. The scheme is based on the lattice Boltzmann method and uses the level set method to compute the evolution of the interface. To couple the fluid phases, we develop new boundary conditions which realise the macroscopic jump conditions at the interface and incorporate surface tension in the lattice Boltzmann framework. Various simulations are presented to validate the numerical scheme, e.g. two-phase channel flows, the Young–Laplace law for a bubble and viscous fingering in a Hele-Shaw cell. The results show that the method is feasible over a wide range of density and viscosity differences.     

Numerical study on transient response of droplet deformation in a steady electric field

The transient response of droplet deformation in a steady electric field is investigated by the numerical simulation and the motion of interface is captured by level-set method. The numerical scheme is validated and found to be in good agreement with classic analytical solutions. The effects of electric field intensity, interfacial tension, oil viscosity and droplet size on the transient deformation process are systematically discussed. The numerical results show that electric field intensity can accelerate the deformation of the droplet, while interfacial tension and oil viscosity damp it. Furthermore, the relation between electric capillary number and dimensionless deformation time is obtained.     

二维卡门涡街的格子Boltzmann仿真

格子气自动机和格子Boltzmann方法的迅速发展提供了一类求解流体力学问题的新的方法。本文中，我们介绍了Boltzmann方法，解决了格子气方法中的缺点，通过选择适当平衡分布及其参数，导出了Navier－Stokers方程，并得到了声速和粘性系数。最后在微机上模拟了在无限长平板流动问题及绕单分离板的流动问题，得到了卡门涡街。结果表该模型有格子气方法及其它的数值方法所没有的优点，计算更精确、更直观、更有效。     

Mechanical behavior of pathological and normal red blood cells in microvascular flow based on modified level-set method

The research of the motion and deformation of the RBCs is important to reveal the mechanism of blood diseases. A numerical method has been developed with level set formulation for elastic membrane immersed in incompressible fluid. The numerical model satisfies mass and energy conservation without the leaking problems in classical Immersed Boundary Method(IBM), at the same time, computing grid we used can be much smaller than the general literatures. The motion and deformation of a red blood cell(including pathological normal status) in microvascular flow are simulated. It is found that the Reynolds number and membrane's stiffness play an important role in the transmutation and oscillation of the elastic membrane. The normal biconcave shape of the RBC is propitious to create high deformation than other pathological shapes. With reduced viscosity of the interior fluid both the velocity of the blood and the deformability of the cell reduced. With increased viscosity of the plasma both the velocity of the blood and the deformability of the cell reduced. The tank treading of the RBC membrane is observed at low enough viscosity contrast in shear flow. The tank tread fixed inclination angle of the cell depends on the shear ratio and viscosity contrast, which can be compared with the experimental observation well.     

Assessment of localized artificial diffusivity scheme for large-eddy simulation of compressible turbulent flows

The localized artificial diffusivity method is investigated in the context of large-eddy simulation of compressible turbulent flows. The performance of different artificial bulk viscosity models are evaluated through detailed results from the evolution of decaying compressible isotropic turbulence with eddy shocklets and supersonic turbulent boundary layer. Effects of subgrid-scale (SGS) models and implicit time-integration scheme/time-step size are also investigated within the framework of the numerical scheme used. The use of a shock sensor along with artificial bulk viscosity significantly improves the scheme for simulating turbulent flows involving shocks while retaining the shock-capturing capability. The proposed combination of Ducros-type sensor with a negative dilatation sensor removes unnecessary bulk viscosity within expansion and weakly compressible turbulence regions without shocks and allows it to localize near the shocks. It also eliminates the need for a wall-damping function for the bulk viscosity while simulating wall-bounded turbulent flows. For the numerical schemes used, better results are obtained without adding an explicit SGS model than with SGS model at moderate Reynolds number. Inclusion of a SGS model in addition to the low-pass filtering and artificial bulk viscosity results in additional damping of the resolved turbulence. However, investigations at higher Reynolds numbers suggest the need for an explicit SGS model. The flow statistics obtained using the second-order implicit time-integration scheme with three sub-iterations closely agrees with the explicit scheme if the maximum Courant–Friedrichs–Lewy is kept near unity.     

A boundary integral method for simulating the dynamics of inextensible vesicles suspended in a viscous fluid in 2D

We present a new method for the evolution of inextensible vesicles immersed in a Stokesian fluid. We use a boundary integral formulation for the fluid that results in a set of nonlinear integro-differential equations for the vesicle dynamics. The motion of the vesicles is determined by balancing the non-local hydrodynamic forces with the elastic forces due to bending and tension. Numerical simulations of such vesicle motions are quite challenging. On one hand, explicit time-stepping schemes suffer from a severe stability constraint due to the stiffness related to high-order spatial derivatives and a milder constraint due to a transport-like stability condition. On the other hand, an implicit scheme can be expensive because it requires the solution of a set of nonlinear equations at each time step. We present two semi-implicit schemes that circumvent the severe stability constraints on the time step and whose computational cost per time step is comparable to that of an explicit scheme. We discretize the equations by using a spectral method in space, and a multistep third-order accurate scheme in time. We use the fast multipole method (FMM) to efficiently compute vesicle–vesicle interaction forces in a suspension with a large number of vesicles. We report results from numerical experiments that demonstrate the convergence and algorithmic complexity properties of our scheme.     

The Immersed Structural Potential Method for haemodynamic applications

In this paper, a new fluid–structure interaction immersed computational methodology, based upon the original Immersed Boundary Method (IBM) [1] is outlined with the final aim of modelling cardiovascular phenomena, specifically, heart valve related problems. The principal characteristic of such immersed techniques is the representation of any deformable or rigid body immersed within an incompressible viscous flow field as a momentum forcing source in the Navier–Stokes equations. A number of shortcomings within the immersed formulation still require further investigation and improvement, including the excessive numerical diffusion caused by the interpolation/spreading process, the need to include realistic viscoelastic composite constitutive models describing more accurately the nature of cardiovascular tissues and also the need to capture more effectively stresses developed at the fluid–structure interface. By following the same philosophy as the original IBM, a more sophisticated formulation is derived in this paper, the “Immersed Structural Potential Method (ISPM)”. The method introduced presents an alternative approach to compute the equivalent fluid–structure interaction forces at the fluid mesh, accounts for a sophisticated viscoelastic fibre-reinforced constitutive model to better describe the mechanics of cardiovascular tissues and utilises a novel time-integration methodology for the computation of the deformation gradient tensor which ensures compliance with the incompressibility constraint. A series of numerical examples will be presented in order to demonstrate the robustness and applicability of this new methodology.     

Schardin问题的数值研究

激波绕过三角楔(Schardin问题)时会产生激波马赫反射与绕射、 三角楔尾涡与涡串等复杂物理现象. 本文利用三阶精度加权基本无振荡(WENO)格式、 结构化矩形网格的自适应加密方法与沉浸边界法对Schardin问题进行了数值模拟. 数值结果清晰地显示了激波与三角楔相互作用, 在楔面发生马赫反射以及在楔角绕射诱导主涡的过程, 并与Schardin等的实验结果及相关数值结果完全符合. 另外, 数值结果还详细反映了先前实验与数值结果没有详细讨论的主涡滑移层上的涡串生成机理, 以及激波与涡串相互作用和产生声波的过程.     

Schardin问题的数值研究

激波绕过三角楔(Schardin问题)时会产生激波马赫反射与绕射、 三角楔尾涡与涡串等复杂物理现象. 本文利用三阶精度加权基本无振荡(WENO)格式、 结构化矩形网格的自适应加密方法与沉浸边界法对Schardin问题进行了数值模拟. 数值结果清晰地显示了激波与三角楔相互作用, 在楔面发生马赫反射以及在楔角绕射诱导主涡的过程, 并与Schardin等的实验结果及相关数值结果完全符合. 另外, 数值结果还详细反映了先前实验与数值结果没有详细讨论的主涡滑移层上的涡串生成机理, 以及激波与涡串相互作用和产生声波的