In this paper, the (G'/G)-expansion method is suggested to establish new exact solutions for fractional differential-difference equations in the sense of modified Riemann-Liouville derivative. The fractional complex transform is proposed to convert a fractional partial differential difference equation into its differential difference equation of integer order. With the aid of symbolic computation, we choose nonlinear lattice equations to illustrate the validity and advantages of the algorithm. It is shown that the proposed algorithm is effective and can be used for many other nonlinear lattice equations in mathematical physics and applied mathematics.

In this paper, a DG (Discontinuous Galerkin) method which has been widely employed in CFD (Computational Fluid Dynamics) is used to solve the two-dimensional time-domain Maxwell’s equations for complex geometries on unstructured mesh. The element interfaces on solid boundary are treated in both curved way and straight way. Numerical tests are performed for both benchmark problems and complex cases with varying orders on a series of grids, where the high-order convergence in accuracy can be observed. Both the curved and the straight solid boundary implementation can give accurate RCS (Radar Cross-Section) results with sufficiently small mesh size, but the curved solid boundary implementation can significantly improve the accuracy when using relatively large mesh size. More importantly, this CFD-based high-order DG method for the Maxwell’s equations is very suitable for complex geometries.

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.

In this paper, we proposes and analyzes the mixed 4th-order Runge-Kutta scheme of conditional nonlinear perturbation (CNOP) approach for the EI Niño-Southern Oscillation (ENSO) model. This method consists of solving the ENSO model by using a mixed 4th-order Runge-Kutta method. Convergence, the local and global truncation error of this mixed 4th-order Runge-Kutta method are proved. Furthermore, optimal control problem is developed and the gradient of the cost function is determined.

The lattice Boltzmann equation (LBE) is considered as a promising approach for simulating flows of liquid and gas. Most of LBE studies have been devoted to regular square LBE and few works have focused on the rectangular LBE in the simulation of gaseous microscale flows. In fact, the rectangular LBE, as an alternative and efficient method, has some advantages over the square LBE in simulating flows with certain computational domains of large aspect ratio (e.g., long micro channels). Therefore, in this paper we expand the application scopes of the rectangular LBE to gaseous microscale flow. The kinetic boundary conditions for the rectangular LBE with a multiple-relaxation-time (MRT) collision operator, i.e., the combined bounce-back/specular-reflection (CBBSR) boundary condition and the discrete Maxwell's diffuse-reflection (DMDR) boundary condition, are studied in detail. We observe some discrete effects in both the CBBSR and DMDR boundary conditions for the rectangular LBE and present a reasonable approach to overcome these discrete effects in the two boundary conditions. It is found that the DMDR boundary condition for the square MRT-LBE can not realize the real fully diffusive boundary condition, while the DMDR boundary condition for the rectangular MRT-LBE with the grid aspect ratio a≠1 can do it well. Some numerical tests are implemented to validate the presented theoretical analysis. In addition, the computational efficiency and relative difference between the rectangular LBE and the square LBE are analyzed in detail. The rectangular LBE is found to be an efficient method for simulating the gaseous microscale flows in domains with large aspect ratios.

This article addresses the resolution of the inverse problem for the parameter identification in orthotropic materials with a number of measurements merely on the boundaries. The inverse problem is formulated as an optimization problem of a residual functional which evaluates the differences between the experimental and predicted displacements. The singular boundary method, an integration-free, mathematically simple and boundary-only meshless method, is employed to numerically determine the predicted displacements. The residual functional is minimized by the Levenberg-Marquardt method. Three numerical examples are carried out to illustrate the robustness, efficiency, and accuracy of the proposed scheme. In addition, different levels of noise are added into the boundary conditions to verify the stability of the present methodology.

The maximum entropy method for the Hausdorff moment problem suffers from ill conditioning as it uses monomial basis {1,x,x2,...,xn}. The maximum entropy method for the Chebyshev moment probelm was studied to overcome this drawback in. In this paper we review and modify the maximum entropy method for the Hausdorff and Chebyshev moment problems studied in and present the maximum entropy method for the Legendre moment problem. We also give the algorithms of converting the Hausdorff moments into the Chebyshev and Lengendre moments, respectively, and utilizing the corresponding maximum entropy method.

A spectral Jacobi-collocation approximation is proposed for Volterra delay integro-differential equations with weakly singular kernels. In this paper, we consider the special case that the underlying solutions of equations are sufficiently smooth. We provide a rigorous error analysis for the proposed method, which shows that both the errors of approximate solutions and the errors of approximate derivatives decay exponentially in L∞ norm and weighted L2 norm. Finally, two numerical examples are presented to demonstrate our error analysis.

In this paper, the idea of a combination of variable separation approach and the extended homoclinic test approach is proposed to seek non-travelling wave solutions of Calogero equation. The equation is reduced to some (1+1)-dimensional nonlinear equations by applying the variable separation approach and solves reduced equations with the extended homoclinic test technique. Based on this idea and with the aid of symbolic computation, some new explicit solutions can be obtained.

We compare in this paper the properties of Osher flux with O-variant and P-variant (Osher-O flux and Osher-P flux) in finite volume methods for the two-dimensional Euler equations and propose an entropy fix technique to improve their robustness. We consider both first-order and second-order reconstructions. For inviscid hypersonic flow past a circular cylinder, we observe different problems for different schemes: A first-order Osher-O scheme on quadrangular grids yields a carbuncle shock, while a first-order Osher-P scheme results in a dislocation shock for high Mach number cases. In addition, a second-order Osher scheme can also yield a carbuncle shock or be unstable. To improve the robustness of these schemes we propose an entropy fix technique, and then present numerical results to show the effectiveness of the proposed method. In addition, the influence of grid aspects ratio, relative shock position to the grid and Mach number on shock stability are tested. Viscous heating problem and double Mach reflection problem are simulated to test the influence of the entropy fix on contact resolution and boundary layer resolution.

The elastic-plastic response of the functionally graded thick-walled tube subjected to internal pressure is investigated by using the relation of the volume average stresses of constituents and the macroscopic stress of composite material in micromechanics. The tube consists of two idealized isotropic elastic-plastic materials whose volume fractions are power functions of the radius. As the internal pressure increases, the deformations of one phase and two phases from elastic to plastic are analyzed. In order to simplify the calculations we assume both materials with the same Poisson's ratio. By using the assumption of a uniform strain field within the representative volume element and the Tresca yield criterion, the theoretical solutions are obtained for the case of two elastic phases and the case of two plastic phases, and the function of the radial displacement is presented for the case with both elastic and plastic phases. The yield criterion of functionally graded material is given in terms of the yield stresses and volume fractions of constituents rather than Young's modulus and yield stress with different unknown parameters of the whole material in the existing papers. Finally we also discuss the position where the plastic deformation first occurs and the conditions for which material first yields in the tube.

The differential quadrature method has been widely used in scientific and engineering computation. However, for the basic characteristics of time domain differential quadrature method, such as numerical stability and calculation accuracy or order, it is still lack of systematic analysis conclusions. In this paper, according to the principle of differential quadrature method, it has been derived and proved that the weighting coefficients matrix of differential quadrature method meets the important V-transformation feature. Through the equivalence of the differential quadrature method and the implicit Runge-Kutta method, it has been proved that the differential quadrature method is A-stable and s-stage s-order method. On this basis, in order to further improve the accuracy of the time domain differential quadrature method, a class of improved differential quadrature method of s-stage 2s-order has been proposed by using undetermined coefficients method and Padé approximations. The numerical results show that the proposed differential quadrature method is more precise than the traditional differential quadrature method.

In this paper, we study carbon emission trading whose market is gaining in popularity as a policy instrument for global climate change. The mathematical model is presented for pricing options on CO2 emission allowance futures with jump diffusion processes, and a so-called fitted finite volume method is proposed to solve the pricing model for the spatial discretization, in which the Crank-Nicolson is employed for time stepping. In addition, the stability and the convergence of the fully discrete scheme are given, and some numerical results, which are compared with the closed form solution and the Monte Carlo simulation solution, are provided to demonstrate the rates of convergence and the robustness of the numerical method.

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 work, the hybrid solution reconstruction formulation proposed by Luo et al. [H. Luo, H. Dai, P. F. de Sousa and B. Yin, On the numerical oscillation of the direct-forcing immersed-boundary method for moving boundaries, Computers & Fluids, 56 (2012), pp. 61–76] for the finite-difference discretization on Cartesian meshes is implemented in the finite-element framework of the local domain-free discretization (DFD) method to reduce the numerical oscillations in the simulation of moving-boundary flows. The reconstruction formulation is applied at fluid nodes in the immediate vicinity of the immersed boundary, which combines weightly the local DFD solution with the specific values obtained via an approximation of quadratic polynomial in the normal direction to the wall. The quadratic approximation is associated with the no-slip boundary condition and the local simplified momentum equation. The weighted factor suitable for unstructured triangular and tetrahedral meshes is constructed, which is related to the local mesh intervals near the immersed boundary and the distances from exterior dependent nodes to the boundary. Therefore, the reconstructed solution can account for the smooth movement of the immersed boundary. Several numerical experiments have been conducted for two- and three-dimensional moving-boundary flows. It is shown that the hybrid reconstruction approach can work well in the finite-element context and effectively reduce the numerical oscillations with little additional computational cost, and the spatial accuracy of the original local DFD method can also be preserved.

This paper extends an algorithm of P1-conservative interpolation on triangular meshes to tetrahedral meshes and thus constructs an approach of solution reconstruction for three-dimensional problems. The conservation property is achieved by local mesh intersection and the mass of a tetrahedron of the current mesh is calculated by the integral on its intersection with the background mesh. For each current tetrahedron, the overlapped background tetrahedrons are detected efficiently. A mesh intersection algorithm is proposed to construct the intersection of a current tetrahedron with the overlapped background tetrahedron and mesh the intersection region by tetrahedrons. A localization algorithm is employed to search the host units in background mesh for each vertex of the current mesh. In order to enforce the maximum principle and avoid the loss of monotonicity, correction of nodal interpolated solution on tetrahedral meshes is given. The performance of the present solution reconstruction method is verified by numerical experiments on several analytic functions and the solution of the flow around a sphere.

In this paper, a high-order curved mesh generation method for Discontinuous Galerkin methods is introduced. First, a regular mesh is generated. Second, the solid surface is re-constructed using cubic polynomial. Third, the elastic governing equations are solved using high-order finite element method to provide a fully or partly curved grid. Numerical tests indicate that the intersection between element boundaries can be avoided by carefully defining the elasticity modulus.

The paper makes a comparative study of the finite element method (FEM) and the finite difference method (FDM) for two-dimensional fractional advection-dispersion equation (FADE) which has recently been considered a promising tool in modeling non-Fickian solute transport in groundwater. Due to the non-local property of integro-differential operator of the space-fractional derivative, numerical solution of FADE is very challenging and little has been reported in literature, especially for high-dimensional case. In order to effectively apply the FEM and the FDM to the FADE on a rectangular domain, a backward-distance algorithm is presented to extend the triangular elements to generic polygon elements in the finite element analysis, and a variable-step vector Grünwald formula is proposed to improve the solution accuracy of the conventional finite difference scheme. Numerical investigation shows that the FEM compares favorably with the FDM in terms of accuracy and convergence rate whereas the latter enjoys less computational effort.

In this paper, we consider a least squares nonconforming finite element of low order for solving the transport equations. We give a detailed overview on the stability and the convergence properties of our considered methods in the stability norm. Moreover, we derive residual type a posteriori error estimates for the least squares nonconforming finite element methods under H–1-norm, which can be used as the error indicators to guide the mesh refinement procedure in the adaptive finite element method. The theoretical results are supported by a series of numerical experiments.

This paper provides a proof of robustness of the restricted additive Schwarz preconditioner with harmonic overlap (RASHO) for the second order elliptic problems with jump coefficients. By analyzing the eigenvalue distribution of the RASHO preconditioner, we prove that the convergence rate of preconditioned conjugate gradient method with RASHO preconditioner is uniform with respect to the large jump and meshsize.