In this work we utilize the boundary integral equation and the Dual Reciprocity Boundary Element Method (DRBEM) for the solution of the steady state convection-diffusion-reaction equations with variable convective coefficients in two-dimension. The DRBEM is a numerical method to transform the domain integrals into the boundary only integrals by using the fundamental solution of Helmholtz equation. Some examples are calculated to confirm the accuracy of the approach. The results obtained by the analytic solutions are in good agreement with ones provided by the DRBEM technique.

In this paper, a Legendre-collocation spectral method is developed for the second order Volterra integro-differential equation with pantograph delay. We provide a rigorous error analysis for the proposed method. The spectral rate of convergence for the proposed method is established in both L2-norm and L∞-norm.

The significance of flow optimization utilizing the lattice Boltzmann (LB) method becomes obvious regarding its advantages as a novel flow field solution method compared to the other conventional computational fluid dynamics techniques. These unique characteristics of the LB method form the main idea of its application to optimization problems. In this research, for the first time, both continuous and discrete adjoint equations were extracted based on the LB method using a general procedure with low implementation cost. The proposed approach could be performed similarly for any optimization problem with the corresponding cost function and design variables vector. Moreover, this approach was not limited to flow fields and could be employed for steady as well as unsteady flows. Initially, the continuous and discrete adjoint LB equations and the cost function gradient vector were derived mathematically in detail using the continuous and discrete LB equations in space and time, respectively. Meanwhile, new adjoint concepts in lattice space were introduced. Finally, the analytical evaluation of the adjoint distribution functions and the cost function gradients was carried out.

Magnetohydrodynamic natural convection heat transfer in a rotating, differentially heated enclosure is studied numerically in this article. The governing equations are in velocity, pressure and temperature formulation and solved using the staggered grid arrangement together with MAC method. The governing parameters considered are the Hartmann number, 0≤Ha≤70, the inclination angle of the magnetic field, 0°≤θ 90°, the Taylor number, 8.9 x 104≤Ta≤1.1 x 106 and the centrifugal force is smaller than the Coriolis force and the both forces were kept below the buoyancy force. It is found that a sufficiently large Lorentz force neutralizes the effect of buoyancy, inertial and Coriolis forces. Horizontal or vertical direction of the magnetic field was most effective in reducing the global heat transfer.

Aiming at complex configuration and large deformation, an efficient hybrid method for dynamic mesh generation is presented in this paper, which is based on Radial Basis Functions (RBFs) and Delaunay graph mapping. Based on the computational mesh, a set of very coarse grid named as background grid is generated firstly, and then the computational mesh can be located at the background grid by Delaunay graph mapping technique. After that, the RBFs method is applied to deform the background grid by choosing partial mesh points on the boundary as the control points. Finally, Delaunay graph mapping method is used to relocate the computational mesh by employing area or volume weight coefficients. By applying different dynamic mesh methods to a moving NACA0012 airfoil, it can be found that the RBFs-Delaunay graph mapping hybrid method is as accurate as RBFs and is as efficient as Delaunay graph mapping technique. Numerical results show that the dynamic meshes for all test cases including one two-dimensional (2D) and two three-dimensional (3D) problems with different complexities, can be generated in an accurate and efficient manner by using the present hybrid method.

We develop a numerical method to simulate a two-phase compressible flow with sharp phase interface on Eulerian grids. The scheme makes use of a levelset to depict the phase interface numerically. The overall scheme is basically a finite volume scheme. By approximately solving a two-phase Riemann problem on the phase interface, the normal phase interface velocity and the pressure are obtained, which is used to update the phase interface and calculate the numerical flux between the flows of two different phases. We adopt an aggregation algorithm to build cell patches around the phase interface to remove the numerical instability due to the breakdown of the CFL constraint by the cell fragments given by the phase interface depicted using the levelset function. The proposed scheme can handle problems with tangential sliping on the phase interface, topological change of the phase interface and extreme contrast in material parameters in a natural way. Though the perfect conservation of the mass, momentum and energy in global is not achieved, it can be quantitatively identified in what extent the global conservation is spoiled. Some numerical examples are presented to validate the numerical method developed.

A novel meshless scheme is proposed for inverse source identification problems of Helmholtz-type equations. It is formulated by the non-singular general solutions of the Helmholtz-type equations augmented with radial basis functions. Under this meshless scheme, we can determine smooth source terms from partially accessible boundary measurements with accurate results. Numerical examples are presented to verify validity and accuracy of the present scheme. It is demonstrated that the present scheme is simple, accurate, stable and computationally efficient for inverse smooth source identification problems.

In this paper, we are devoted to nonlinear stability and B-convergence of additive Runge-Kutta (ARK) methods for nonlinear stiff problems with multiple stiffness. The concept of -algebraic stability of ARK methods for a class of stiff problems Kστ is introduced, and it is proven that this stability implies some contractive properties of the ARK methods. Some results on optimal B-convergence of ARK methods for Kσ,0 are given. These new results extend the existing ones of RK methods and ARK methods in the references. Numerical examples test the correctness of our theoretical analysis.

The linear operator plays an important role in the computational process of Homotopy Analysis Method (HAM). In HAM frame any kind of linear operator can be chosen to develop a solution. Hence, it is easy to introduce the modified/physical parameter dependent linear operators. The effective use of these operators has been judged through solving fluid flow problems. Modification in linear operators affects the solution and improves the computational efficiency of HAM for larger values of parameters. The convergence rate of the solution is rapid and several times higher resulting in lesser computational time.

A vortex ring impacting a three-dimensional bump is studied using large eddy simulation for a Reynolds number Re = 4 × 104 based on the initial diameter and translational speed of the vortex ring. The effects of bump height and vortex core thickness for thin and thick vortex rings on the vortical flow phenomena and the underlying physical mechanisms are investigated. Based on the analysis of the evolution of vortical structures, two typical kinds of vortical structures, i.e., the wrapping vortices and the hair-pin vortices, are identified and play an important role in the flow state evolution. The boundary vorticity flux is analyzed to reveal the mechanism of the vorticity generation on the bump surface. The circulation of the primary vortex ring reasonably elucidates some typical phases of flow evolution. Further, the analysis of turbulent kinetic energy reveals the transition from laminar to turbulent state. The results obtained in this study provide physical insight into the understanding of the mechanisms relevant to the flow evolution and the flow transition to turbulent state.

Interaction of vortex rings with solid is an important research topic of hydrodynamic. In this study, a multiple-relaxation time (MRT) lattice Boltzmann method (LBM) is used to investigate the flow of a vortex ring impacting spheroidal particles. The MRT-LBM is validated through the cases of vortex ring impacting a flat wall. The vortex evolution due to particle size, the aspect ratio of a prolate particle, as well as Reynolds (Re) number are discussed in detail. When the vortex ring impacting a stationary sphere, the primary and secondary vortex rings wrap around each other, which is different from the situation of the vortex ring impacting a plate. For the vortex ring impacting with a prolate spheroid, the secondary vortex ring stretches mainly along the long axis of the ellipsoid particle. However, it is found that after the vortex wrapping stage, the primary vortex recovers along the short axis of the particle faster than that in the long axis, i.e., the primary vortex ring stretches mainly along the short axis of the particle. That has never been address in the literature.

Aerosol modeling is very important to study the behavior of aerosol dynamics in atmospheric environment. In this paper we consider numerical methods for the nonlinear aerosol dynamic equations on time and particle size. The finite volume element methods based on the linear interpolation and Hermite interpolation are provided to approximate the aerosol dynamic equation where the condensation and removal processes are considered. Numerical examples are provided to show the efficiency of these numerical methods.

In this paper we develop explicit fast exponential Runge-Kutta methods for the numerical solutions of a class of parabolic equations. By incorporating the linear splitting technique into the explicit exponential Runge-Kutta schemes, we are able to greatly improve the numerical stability. The proposed numerical methods could be fast implemented through use of decompositions of compact spatial difference operators on a regular mesh together with discrete fast Fourier transform techniques. The exponential Runge-Kutta schemes are easy to be adopted in adaptive temporal approximations with variable time step sizes, as well as applied to stiff nonlinearity and boundary conditions of different types. Linear stabilities of the proposed schemes and their comparison with other schemes are presented. We also numerically demonstrate accuracy, stability and robustness of the proposed method through some typical model problems.

In this paper, a method is proposed for extracting fracture parameters in anisotropic thermoelasticity cracking via interaction integral method within the framework of extended finite element method (XFEM). The proposed method is applied to linear thermoelastic crack problems. The numerical results of the stress intensity factors (SIFs) are presented and compared with those reported in related references. The good agreement of the results obtained by the developed method with those obtained by other numerical solutions proves the applicability of the proposed approach and confirms its capability of efficiently extracting thermoelasticity fracture parameters in anisotropic materials.

This paper presents a physical investigation and mathematical analysis on mechanical behavior of the regular jugged discontinuity. In particular, we focus on the creep property of structural plane with various slope angles under different normal stress through shear creep tests of structural plane under shear stresses. According to the test results, the shear creep property of structural plane was described and the creep velocity and long-term strength of the structural plane during shear creep were also investigated. An empirical formula is finally established to evaluate shear strength of discontinuity and a modified Burger model was proposed to represent the shear deformation property during creep.

In this paper, we investigate the error estimates of mixed finite element methods for optimal control problems governed by general elliptic equations. The state and co-state are approximated by the lowest order Raviart-Thomas mixed finite element spaces and the control variable is approximated by piecewise constant functions. We derive L2 and H–1-error estimates both for the control variable and the state variables. Finally, a numerical example is given to demonstrate the theoretical results.

In this paper, we propose a fast proximity point algorithm and apply it to total variation (TV) based image restoration. The novel method is derived from the idea of establishing a general proximity point operator framework based on which new first-order schemes for total variation (TV) based image restoration have been proposed. Many current algorithms for TV-based image restoration, such as Chambolle’s projection algorithm, the split Bregman algorithm, the Bermúdez-Moreno algorithm, the Jia-Zhao denoising algorithm, and the fixed point algorithm, can be viewed as special cases of the new first-order schemes. Moreover, the convergence of the new algorithm has been analyzed at length. Finally, we make comparisons with the split Bregman algorithm which is one of the best algorithms for solving TV-based image restoration at present. Numerical experiments illustrate the efficiency of the proposed algorithms.

In this paper, we consider the transform magnetic (TM) model of electromagnetic scattering in the cavity. By the Polynomial Preserving Recovery technique, we present superconvergence analysis for the vertex-edge-face type finite element. From the numerical example, we can see that the provided method is efficient and stable.

This study investigates the applicability of the singular boundary method (SBM), a recent developed meshless boundary collocation method, for the analysis of two-dimensional (2D) thin structural problems. The troublesome nearly-singular kernels, which are crucial in the applications of SBM to thin shapes, are dealt with efficiently by using a non-linear transformation technique. Promising SBM results with only a small number of boundary nodes are obtained for thin structures with the thickness-to-length ratio is as small as 1E-9, which is sufficient for modeling most thin layered coating systems as used in smart materials and micro-electro-mechanical systems. The advantages, disadvantages and potential applications of the proposed method, as compared with the finite element (FEM) and boundary element methods (BEM), are also discussed.

Intrigued by our recent experimental work (H. Yamaguchi and X. D. Niu, J. Fluids Eng., 133 (2011), 041302), the present study numerically investigate the flow-structure interactions (FSI) of three rigid circular particles aligned moving in an inclined channel flow at intermediate Reynolds numbers by using a momentum-exchanged immersed boundary-lattice Boltzmann method. A ”frog-leap” phenomenon observed in the experiment is successfully captured by the present simulation and flow characteristics and underlying FSI mechanisms of it are explored by examining the effects of the channel inclined angles and Reynolds numbers. It is found that the asymmetric difference of the vorticity distributions on the particle surface is the main cause of the “frog-leap” when particle moves in the boundary layer near the lower channel boundary.