A boundary condition-implemented immersed boundary-lattice Boltzmann method (IB-LBM) is presented in this work. The present approach is an improvement to the conventional IB-LBM. In the conventional IB-LBM, the no-slip boundary condition is only approximately satisfied. As a result, there is flow penetration to the solid boundary. Another drawback of conventional IB-LBM is the use of Dirac delta function interpolation, which only has the first order of accuracy. In this work, the no-slip boundary condition is directly implemented, and used to correct the velocity at two adjacent mesh points from both sides of the boundary point. The velocity correction is made through the second-order polynomial interpolation rather than the first-order delta function interpolation. Obviously, the two drawbacks of conventional IB-LBM are removed in the present study. Another important contribution of this paper is to present a simple way to compute the hydrodynamic forces on the boundary from Newton’s second law. To validate the proposed method, the two-dimensional vortex decaying problem and incompressible flow over a circular cylinder are simulated. As shown in the present results, the flow penetration problem is eliminated, and the obtained results compare very well with available data in the literature.

This paper presents an analytical solution for the interaction of electric potentials, electric displacement, elastic deformations, and describes hygrothermal effect responses in hollow and solid cylinders, subjected to mechanical load and electric potential. Exact solutions for displacement, stresses and electric potentials in functionally graded piezoelectric material are determined using the infinitesimal theory. The material properties coefficients of the present cylinder are assumed to be graded in the radial direction by a power law distribution. Numerical examples display the significant of influence of material inhomogeneity. It is interesting to note that selecting a specific value of inhomogeneity parameter can optimize the piezoelectric hollow and solid cylinders responses, which will be of particular importance in modern engineering designs.

It is very important to predict the bypass transition in the simulation of flows through turbomachinery. This paper presents a four-equation eddy-viscosity turbulence transition model for prediction of bypass transition. It is based on the SST turbulence model and the laminar kinetic energy concept. A transport equation for the non-turbulent viscosity is proposed to predict the development of the laminar kinetic energy in the pre-transitional boundary layer flow which has been observed in experiments. The turbulence breakdown process is then captured with an intermittency transport equation in the transitional region. The performance of this new transition model is validated through the experimental cases of T3AM, T3A and T3B. Results in this paper show that the new transition model can reach good agreement in predicting bypass transition, and is compatible with modern CFD software by using local variables.

In this paper, we consider a two-scale stabilized finite volume method for the two-dimensional stationary incompressible flow approximated by the lowest equal-order element pair P1 — P1 which do not satisfy the inf-sup condition. The two-scale method consist of solving a small non-linear system on the coarse mesh and then solving a linear Stokes equations on the fine mesh. Convergence of the optimal order in the H1-norm for velocity and the L2-norm for pressure are obtained. The error analysis shows there is the same convergence rate between the two-scale stabilized finite volume solution and the usual stabilized finite volume solution on a fine mesh with relation h = 𝾪 (H2). Numerical experiments completely confirm theoretic results. Therefore, this method presented in this paper is of practical importance in scientific computation.

In this paper, a conservative fifth-order upwind compact scheme using centered stencil is introduced. This scheme uses asymmetric coefficients to achieve the upwind property since the stencil is symmetric. Theoretical analysis shows that the proposed scheme is low-dissipative and has a relatively large stability range. To maintain the convergence rate of the whole spatial discretization, a proper non-periodic boundary scheme is also proposed. A detailed analysis shows that the spatial discretization implemented with the boundary scheme proposed by Pirozzoli [J. Comput. Phys., 178 (2001), pp. 81–117] is approximately fourth-order. Furthermore, a hybrid methodology, coupling the compact scheme with WENO scheme, is adopted for problems with discontinuities. Numerical results demonstrate the effectiveness of the proposed scheme.

In this paper, we develop an accurate and efficient Legendre wavelets method for numerical solution of the well known time-fractional telegraph equation. In the proposed method we have employed both of the operational matrices of fractional integration and differentiation to get numerical solution of the time-telegraph equation. The power of this manageable method is confirmed. Moreover the use of Legendre wavelet is found to be accurate, simple and fast.

In this paper, high-order Discontinuous Galerkin (DG) method is used to solve the two-dimensional Euler equations. A shock-capturing method based on the artificial viscosity technique is employed to handle physical discontinuities. Numerical tests show that the shocks can be captured within one element even on very coarse grids. The thickness of the shocks is dominated by the local mesh size and the local order of the basis functions. In order to obtain better shock resolution, a straightforward hp-adaptivity strategy is introduced, which is based on the high-order contribution calculated using hierarchical basis. Numerical results indicate that the hp-adaptivity method is easy to implement and better shock resolution can be obtained with smaller local mesh size and higher local order.

We propose two variants of tailored finite point (TFP) methods for discretizing two dimensional singular perturbed eigenvalue (SPE) problems. A continuation method and an iterative method are exploited for solving discretized systems of equations to obtain the eigen-pairs of the SPE. We study the analytical solutions of two special cases of the SPE, and provide an asymptotic analysis for the solutions. The theoretical results are verified in the numerical experiments. The numerical results demonstrate that the proposed schemes effectively resolve the delta function like of the eigenfunctions on relatively coarse grid.

In this paper, we study wave interaction with an emerged porous media. The governing equation is shallow water equations with a friction term of the linearized Dupuit-Forcheimer’s formula. From the continuity of surface and horizontal flux, we derived the wave reflection and transmission coefficient formulas. They are similar with the corresponding formulas of the submerged solid bar breakwater. We solve the equations numerically using finite volume method on a staggered grid. The numerical wave reduction in the porous media confirms the analytical wave transmission curve.

In this paper, we construct a global repair technique for the finite element scheme of anisotropic diffusion equations to enforce the repaired solutions satisfying the discrete maximum principle. It is an extension of the existing local repair technique. Both of the repair techniques preserve the total energy and are easy to be implemented. The numerical experiments show that these repair techniques do not destroy the accuracy of the finite element scheme, and the computational cost of the global repair technique is cheaper than the local repair technique when the diffusion tensors are highly anisotropic.

The Courant-Friedrichs-Lewy condition (The CFL condition) is appeared in the analysis of the finite difference method applied to linear hyperbolic partial differential equations. We give a remark on the CFL condition from a view point of stability, and we give some numerical experiments which show instability of numerical solutions even under the CFL condition. We give a mathematical model for rounding errors in order to explain the instability.