On the stability of the finite element immersed boundary method

The immersed boundary (IB) method is a mathematical formulation for fluid–structure interaction problems, where immersed incompressible visco-elastic bodies or boundaries interact with an incompressible fluid.The original numerical scheme associated to the IB method requires a smoothed approximation of the Dirac delta distribution to link the moving Lagrangian domain with the fixed Eulerian one.We present a stability analysis of the finite element immersed boundary method, where the Dirac delta distribution is treated variationally, in a generalized visco-elastic framework and for two different time-stepping schemes.     

Incompressible flow computations over moving boundary using a novel upwind method

In this paper a novel method for simulating unsteady incompressible viscous flow over a moving boundary is described. The numerical model is based on a 2D Navier–Stokes incompressible flow in artificial compressibility formulation with Arbitrary Lagrangian Eulerian approach for moving grid and dual time stepping approach for time accurate discretization. A higher order unstructured finite volume scheme, based on a Harten Lax and van Leer with Contact (HLLC) type Riemann solver for convective fluxes, developed for steady incompressible flow in artificial compressibility formulation by Mandal and Iyer (AIAA paper 2009-3541), is extended to solve unsteady flows over moving boundary. Viscous fluxes are discretized in a central differencing manner based on Coirier’s diamond path. An algorithm based on interpolation with radial basis functions is used for grid movements. The present numerical scheme is validated for an unsteady channel flow with a moving indentation. The present numerical results are found to agree well with experimental results reported in literature.     

A novel immersed boundary procedure for flow and heat simulations with moving boundary

This article describes a novel immersed boundary procedure for computing the flow and heat transfer problems with moving and complex boundary. Although the immersed boundary techniques have been successfully implemented to these flow and heat simulations, a frequently encountered drawback of this method is the relatively low accuracy proximate to the boundary due to the spreading of forcing function or the interpolation scheme. In this study, we propose a moving-grid process under the arbitrary Lagrangian-Eulerian framework to reduce the numerical diffusion near the immersed boundary. The incompressible Navier-Stokes equations are discretized spatially using unstructured finite element method, and advanced temporally by an operator-splitting scheme. The methodology is validated by the simulations of flow induced by an oscillating cylinder in a free stream. The capability of the proposed method is further demonstrated by good predictions of flow passing the rotating fan in a channel and also flow driven by two independent rotating fans in a circular cavity.     

Determination of a control parameter in a one-dimensional parabolic equation using the moving least-square approximation

In this paper the approximation of moving least-square (MLS) is used for finding the solution of a one-dimensional parabolic inverse problem with source control parameter. Comparing with other numerical methods based on meshes such as finite difference method, finite element method and boundary element method, etc. the MLS approximation has merits of simpler numerical procedures, lower computation cost and arbitrary nodes. The result of a numerical example is presented.     

Large-eddy simulation in the near-field of a transient multi-component gas jet with density gradients

Large-eddy simulation (LES) is a research tool that is increasingly being used to study practical engineering flows because of continuous improvements in computational power. This paper outlines an LES model developed for the study of multi-component transient gas jets with density gradients. The compressible LES formulation together with the numerical model, boundary conditions, perturbation, and parallelization are discussed. A non-dissipative sixth-order finite difference scheme is used to discretize the governing equations, and a low-pass sixth-order spatial filtering scheme is employed to avoid the growth of high-frequency modes. The conditions at the boundaries are implemented using Navier–Stokes characteristic boundary conditions. In addition to code performance, results are presented from a study of an impulsively started jet at high pressure and temperature with an injected to ambient gas density ratio of approximately 3.5.     

Modelling and simulation of moving contact line problems with wetting effects

This paper presents a numerical scheme for computing moving contact line flows with wetting effects. The numerical scheme is based on Arbitrary Lagrangian Eulerian (ALE) finite elements on moving meshes. In the computations, the wetting effects are taken into account through a weak enforcement of the prescribed equilibrium contact angle into the model equations. The equilibrium contact angle is included in the variational form of the model by replacing the curvature with Laplace Beltrami operator and integration by parts. This weak implementation allows that the contact angle determined by the numerical scheme differs from the equilibrium value and develops a certain dynamics. The Laplace Beltrami operator technique with an interface/boundary resolved mesh is well-suited for describing the dynamic contact angle observed in experiments. We consider the spreading and the pendant liquid droplets to investigate this implementation of the contact angle. It is shown that the dynamic contact angle tends to the prescribed equilibrium contact angle when time goes to infinity. However, the dynamics of the contact angle is influenced by the slip at the moving contact line. This work has been partially supported by the German Research Foundation (DFG) through the grant To143/9.     

A Galerkin collocation method for some integral equations of the first kind

Here we present a certain modified collocation method which is a fully discretized numerical method for the solution of Fredholm integral equations of the first kind with logarithmic kernel as principal part. The scheme combines high accuracy from Galerkin's method with the high speed of collocation methods. The corresponding asymptotic error analysis shows optimal order of convergence in the sense of finite element approximation. The whole method is an improved boundary integral method for a wide class of plane boundary value problems involving finite element approximations on the boundary curve. The numerical experiments reveal both, high speed and high accuracy.     

A new fourth-order difference scheme for solving an N-carrier system with Neumann boundary conditions

In this paper, a numerical method is developed to solve an N-carrier system with Neumann boundary conditions. First, we apply the compact finite difference scheme of fourth order for discretizing spatial derivatives at the interior points. Then, we develop a new combined compact finite difference scheme for the boundary, which also has fourth-order accuracy. Lastly, by using a Padé approximation method for the resulting linear system of ordinary differential equations, a new compact finite difference scheme is obtained. The present scheme has second-order accuracy in time direction and fourth-order accuracy in space direction. It is shown that the scheme is unconditionally stable. The present scheme is tested by two numerical examples, which show that the convergence rate with respect to the spatial variable from the new scheme is higher and the solution is much more accurate when compared with those obtained by using other previous methods.     

An Adaptive Moving Mesh Finite Element Solution of the Regularized Long Wave Equation

An adaptive moving mesh finite element method is proposed for the numerical solution of the regularized long wave (RLW) equation. A moving mesh strategy based on the so-called moving mesh PDE is used to adaptively move the mesh to improve computational accuracy and efficiency. The RLW equation represents a class of partial differential equations containing spatial-time mixed derivatives. For the numerical solution of those equations, a \(C^0\) finite element method cannot apply directly on a moving mesh since the mixed derivatives of the finite element approximation may not be defined. To avoid this difficulty, a new variable is introduced and the RLW equation is rewritten into a system of two coupled equations. The system is then discretized using linear finite elements in space and the fifth-order Radau IIA scheme in time. A range of numerical examples in one and two dimensions, including the RLW equation with one or two solitary waves and special initial conditions that lead to the undular bore and solitary train solutions, are presented. Numerical results demonstrate that the method has a second order convergence and is able to move and adapt the mesh to the evolving features in the solution.     

Simultaneous Refinement and Coarsening for Adaptive Meshing

In the numerical simulation of the combustion process and microstructural evolution, we need to consider the adaptive meshing problem for a domain that has a moving boundary. During the simulation, the region ahead of the moving boundary needs to be refined (to satisfy stronger numerical conditions), and the submesh in the region behind the moving boundary should be coarsened (to reduce the mesh size). We present a unified scheme for simultaneously refining and coarsening a mesh. Our method uses sphere packings and guarantees that the resulting mesh is well-shaped and is within a constant factor of the optimal possible in the number of mesh elements. We also present several practical variations of our provably good algorithm.     

An energy-preserving algorithm for nonlinear Hamiltonian wave equations with Neumann boundary conditions

In this paper, a novel energy-preserving numerical scheme for nonlinear Hamiltonian wave equations with Neumann boundary conditions is proposed and analyzed based on the blend of spatial discretization by finite element method (FEM) and time discretization by Average Vector Field (AVF) approach. We first use the finite element discretization in space, which leads to a system of Hamiltonian ODEs whose Hamiltonian can be thought of as the semi-discrete energy of the original continuous system. The stability of the semi-discrete finite element scheme is analyzed. We then apply the AVF approach to the Hamiltonian ODEs to yield a new and efficient fully discrete scheme, which can preserve exactly (machine precision) the semi-discrete energy. The blend of FEM and AVF approach derives a new and efficient numerical scheme for nonlinear Hamiltonian wave equations. The numerical results on a single-soliton problem and a sine-Gordon equation are presented to demonstrate the remarkable energy-preserving property of the proposed numerical scheme.     

The Local Artificial Boundary Conditions for Numerical Simulations of the Flow Around a Submerged Body

We consider in this work the numerical approximations of the two-dimensional steady potential flow around a body moving in a liquid of finite constant depth at constant speed and distance below a free surface. Several vertical segments are introduced as the upstream and the downstream artificial boundaries, where a sequence of high-order local artificial boundary conditions are proposed. Then the original problem is solved in a finite computational domain, which is equivalent to a variational problem. The numerical approximations for the original problem are obtained by solving the variational problem with the finite element method. The numerical examples show that the artificial boundary conditions given in this work are very effective.     

Hybrid Spectral Difference Methods for Elliptic Equations on Exterior Domains with the Discrete Radial Absorbing Boundary Condition

The hybrid spectral difference methods (HSD) for the Laplace and Helmholtz equations in exterior domains are proposed. We consider the fictitious domain method with the absorbing boundary conditions (ABCs). The HSD method is a finite difference version of the hybridized Galerkin method, and it consists of two types of finite difference approximations; the cell finite difference and the interface finite difference. The fictitious domain is composed of two subregions; the Cartesian grid region and the boundary layer region in which the radial grid is imposed. The boundary layer region with the radial grid makes it easy to implement the discrete radial ABC. The discrete radial ABC is a discrete version of the Bayliss–Gunzburger–Turkel ABC without pertaining any radial derivatives. Numerical experiments confirming efficiency of our numerical scheme are provided.     

A capillary free boundary problem governed by the Navier-Stokes equations

In this paper an unsteady capillary free boundary in a viscous imcompressible fluid is studied from the numerical point of view. The linearized problem is formulated, and existence and uniqueness of a solution is proved. The linearized problem is discretized by a finite element method in space and a finite difference method in time. The numerical scheme is proved to be convergent and stable. The method is applied to the computation of free boundaries, eigenfrequencies, damping factors and vibration modes. The influence of the initial conditions, the boundary conditions, the contact angle, the Ohnesorge number and the Bond number is studied.     

A new ADI technique for two-dimensional parabolic equation with an integral condition

Nonclassical parabolic initial-boundary value problems arise in the study of several important physical phenomena. This paper presents a new approach to treat complicated boundary conditions appearing in the parabolic partial differential equations with nonclassical boundary conditions. A new fourth-order finite difference technique, based upon the Noye and Hayman (N-H) alternating direction implicit (ADI) scheme, is used as the basis to solve the two-dimensional time dependent diffusion equation with an integral condition replacing one boundary condition. This scheme uses less central processor time (CPU) than a second-order fully implicit scheme based on the classical backward time centered space (BTCS) method for two-dimensional diffusion. It also has a larger range of stability than a second-order fully explicit scheme based on the classical forward time centered space (FTCS) method. The basis of the analysis of the finite difference equations considered here is the modified equivalent partial differential equation approach, developed from the 1974 work of Warming and Hyeet. This allows direct and simple comparison of the errors associated with the equations as well as providing a means to develop more accurate finite difference methods. The results of numerical experiments for the new method are presented. The central processor times needed are also reported. Error estimates derived in the maximum norm are tabulated.     

A numerical method for a kinetic equation and its application to propagating shock waves

An efficient numerical method for kinetic equations and its application to analyses of moving shock wave problems are presented. The present study aims to give an efficient scheme for two-dimensional unsteady gas flows. An explicit MacCormack difference method is applied to solve a BGK-model equation. The efficiency and accuracy of the scheme are examined in an application to one-dimensional shock structure problems. Furthermore, the scheme is applied to a two-dimensional flow problem: nonstationary reflection of a shock wave at a wedge. The present scheme is found to be useful and efficient for the analyses of two-dimensional unsteady rarefied gas flows.     

The use of continuous boundary elements in the boundary elements method for domains with non-smooth boundaries via finite difference approach

A numerical method is presented in this article to deal with the drawback of boundary elements method (BEM) at corner points. The use of continuous elements instead of the discontinuous ones has been recommended in the BEM literature widely because of the simplicity and accuracy. However the continuous elements lead to certain difficulties for problems where their domains contain corners. In this paper the finite difference method (FDM) has been applied to obtain some constraints for boundary points near the corners to deal with this drawback. Because of its simplicity and capability, the new scheme is applicable on BEM problems for all geometries, all governing equations and general boundary conditions, easily. Since the Dirichlet boundary condition is more critical than the other ones, we will focus on it in the numerical implementation. The numerical results show that the new treatment improves the accuracy of BEM significantly.     

Numerical solution of potential flow equations with a predictor-corrector finite difference method

We develop a numerical solution algorithm of the nonlinear potential flow equations with the nonlinear free surface boundary condition.A finite difference method with a predictor-corrector method is applied to solve the nonlinear potential flow equations in a two-dimensional (2D) tank.The irregular tank is mapped onto a fixed square domain with rectangular cells through a proper mapping function.A staggered mesh system is adopted in a 2D tank to capture the wave elevation of the transient fluid.The finite difference method with a predictor-corrector scheme is applied to discretize the nonlinear dynamic boundary condition and nonlinear kinematic boundary condition.We present the numerical results of wave elevations from small to large amplitude waves with free oscillation motion,and the numerical solutions of wave elevation with horizontal excited motion.The beating period and the nonlinear phenomenon are very clear.The numerical solutions agree well with the analytical solutions and previously published results.     

A direct-forcing immersed boundary method for the thermal lattice Boltzmann method

In this study, a direct-forcing immersed boundary method (IBM) for thermal lattice Boltzmann method (TLBM) is proposed to simulate the non-isothermal flows. The direct-forcing IBM formulas for thermal equations are derived based on two TLBM models: a double-population model with a simplified thermal lattice Boltzmann equation (Model 1) and a hybrid model with an advection–diffusion equation of temperature (Model 2). As an interface scheme, which is required due to a mismatch between boundary and computational grids in the IBM, the sharp interface scheme based on second-order bilinear and linear interpolations (instead of the diffuse interface scheme, which uses discrete delta functions) is adopted to obtain the more accurate results. The proposed methods are validated through convective heat transfer problems with not only stationary but also moving boundaries – the natural convection in a square cavity with an eccentrically located cylinder and a cold particle sedimentation in an infinite channel. In terms of accuracy, the results from the IBM based on both models are comparable and show a good agreement with those from other numerical methods. In contrast, the IBM based on Model 2 is more numerically efficient than the IBM based on Model 1.     

浅水动边界问题的位移法模拟

为精确模拟浅水波非线性演化过程中的动边界,提出一种基于位移的Hamilton变分原理,并进而导出一种基于位移的浅水方程（Shallow Water Equation based on Displacement,SWE D）.SWE D以位移为基本未知量,可以精确满足动边界处的零水深要求并精确捕捉动态边界位置,且解具有协调性.在Hamilton变分原理的框架下,分别采用有限元和保辛积分算法对该浅水方程进行空间离散和时间积分,可有效地处理不平水底情况,保证对非线性演化进行长时间仿真的精度.数值算例表明该方法适用于浅水动边界问题的数值模拟.