The second order weighted and shifted Grünwald difference (WSGD) operators are developed in [Tian, Zhou and Deng, Math. Comput., 84 (2015), pp. 1703–1727] to solve space fractional partial differential equations. Along this direction, we further design a new family of second order WSGD operators; by properly choosing the weighted parameters, they can be effectively used to discretize space (Riemann-Liouville) fractional derivatives. Based on the new second order WSGD operators, we derive a family of difference schemes for the space fractional advection diffusion equation. By von Neumann stability analysis, it is proved that the obtained schemes are unconditionally stable. Finally, extensive numerical experiments are performed to demonstrate the performance of the schemes and confirm the convergence orders.

In this paper, the magnetic field effects on natural convection of power-law nanofluids in rectangular enclosures are investigated numerically with the lattice Boltzmann method. The fluid in the cavity is a water-based nanofluid containing Cu nanoparticles and the investigations are carried out for different governing parameters including Hartmann number (0.0≤Ha≤20.0), Rayleigh number (104≤Ra≤106), power-law index (0.5≤n≤1.0), nanopartical volume fraction (0.0≤ϕ≤0.1) and aspect ratio (0.125≤AR≤8.0). The results reveal that the flow oscillations can be suppressed effectively by imposing an external magnetic field and the augmentation of Hartmann number and power-law index generally decreases the heat transfer rate. Additionally, it is observed that the average Nusselt number is increased with the increase of Rayleigh number and nanoparticle volume fraction. Moreover, the present results also indicate that there is a critical value for aspect ratio at which the impact on heat transfer is the most pronounced.

In this paper, we consider the numerical method that is proposed and analyzed in [J. Cao and C. Xu, J. Comput. Phys., 238 (2013), pp. 154–168] for the fractional ordinary differential equations. It is based on the so-called block-by-block approach, which is a common method for the integral equations. We extend the technique to solve the nonlinear system of fractional ordinary differential equations (FODEs) and present a general technique to construct high order schemes for the numerical solution of the nonlinear coupled system of fractional ordinary differential equations (FODEs). By using the present method, we are able to construct a high order schema for nonlinear system of FODEs of the order α, α>0. The stability and convergence of the schema is rigorously established. Under the smoothness assumption f,g ∈ C 4[0,T], we prove that the numerical solution converges to the exact solution with order 3+α for 0<α≤1 and order 4 for α>1. Some numerical examples are provided to confirm the theoretical claims.

In this paper, we extend using the Runge-Kutta discontinuous Galerkin method together with the front tracking method to simulate the compressible two-medium flow on unstructured meshes. A Riemann problem is constructed in the normal direction in the material interfacial region, with the goal of obtaining a compact, robust and efficient procedure to track the explicit sharp interface precisely. Extensive numerical tests including the gas-gas and gas-liquid flows are provided to show the proposed methodologies possess the capability of enhancing the resolutions nearby the discontinuities inside of the single medium flow and the interfacial vicinities of the two-medium flow in many occasions.

A new Semi-Lagrangian scheme is proposed to discretize the surface convection-diffusion equation. The other involved equations including the the level-set convection equation, the re-initialization equation and the extension equation are also solved by S-L schemes. The S-L method removes both the CFL condition and the stiffness caused by the surface Laplacian, allowing larger time step than the Eulerian method. The method is extended to the block-structured adaptive mesh. Numerical examples are given to demonstrate the efficiency of the S-L method.

Inspired by the behavior of the blind for hill-climbing using a stick to detect a higher place by drawing a circle, we propose a heuristic direct search method to solve the unconstrained optimization problems. Instead of searching a neighbourhood of the current point as done in the traditional hill-climbing, or along specified search directions in standard direct search methods, the new algorithm searches on a surface with radius determined by the motion of the stick. The significant feature of the proposed algorithm is that it only has one parameter, the search radius, which makes the algorithm convenient in practical implementation. The developed method can shrink the search space to a closed ball, or seek for the final optimal point by adjusting search radius. Furthermore our algorithm possesses multi-resolution feature to distinguish the local and global optimum points with different search radii. Therefore, it can be used by itself or integrated with other optimization methods flexibly as a mathematical optimization technique. A series of numerical tests, including high-dimensional problems, have been well designed to demonstrate its performance.

In this work, we study numerical methods for a coupled fluid-porous media flow model. The model consists of Stokes equations and Darcy's equations in two neighboring subdomains, coupling together through certain interface conditions. The weak form for the coupled model is of saddle point type. A mortar finite element method is proposed to approximate the weak form of the coupled problem. In our method, nonconforming Crouzeix-Raviart elements are applied in the fluid subdomain and the lowest order Raviart-Thomas elements are applied in the porous media subdomain; Meshes in different subdomains are allowed to be nonmatching on the common interface; Interface conditions are weakly imposed via adding constraint in the definition of the finite element space. The well-posedness of the discrete problem and the optimal error estimate for the proposed method are established. Numerical experiments are also given to confirm the theoretical results.

In this article a comparison study of the numerical methods for compressible two-phase flows is presented. Although many numerical methods have been developed in recent years to deal with the jump conditions at the fluid-fluid interfaces in compressible multiphase flows, there is a lack of a detailed comparison of these methods. With this regard, the transport five equation model, the modified ghost fluid method and the cut-cell method are investigated here as the typical methods in this field. A variety of numerical experiments are conducted to examine their performance in simulating inviscid compressible two-phase flows. Numerical experiments include Richtmyer-Meshkov instability, interaction between a shock and a rectangle SF6 bubble, Rayleigh collapse of a cylindrical gas bubble in water and shock-induced bubble collapse, involving fluids with small or large density difference. Based on the numerical results, the performance of the method is assessed by the convergence order of the method with respect to interface position, mass conservation, interface resolution and computational efficiency.

In this paper, we analyze a nonconforming finite element method for the computation of transmission eigenvalues and the corresponding eigenfunctions. The error estimates of the eigenvalue and eigenfunction approximation are given, respectively. Finally, some numerical examples are provided to validate the theoretical results.

An inverse problem of determining unknown source parameter in a parabolic equation is considered. The variational iteration method (VIM) is presented to solve inverse problems. The solution gives good approximations by VIM. A numerical example shows that the VIM works effectively for an inverse problem.

In this paper, we introduce the Hamiltonian boundary value method (HBVM) to solve nonlinear Hamiltonian PDEs. We use the idea of Fourier pseudospectral method in spatial direction, which leads to the finite-dimensional Hamiltonian system. The HBVM, which can preserve the Hamiltonian effectively, is applied in time direction. Then the nonlinear Schrödinger (NLS) equation and the Korteweg-de Vries (KdV) equation are taken as examples to show the validity of the proposed method. Numerical results confirm that the proposed method can simulate the propagation and collision of different solitons well. Meanwhile the corresponding errors in Hamiltonian and other intrinsic invariants are presented to show the good preservation property of the proposed method during long-time numerical calculation.

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.

The equation of state (EOS) plays a crucial role in hyperbolic conservation laws for the compressible fluid. Whereas, the solid constitutive model with elastic-plastic phase transition makes the analysis of the solid Riemann problem more difficult. In this paper, one-dimensional elastic-perfectly plastic solid Riemann problem is investigated and its exact Riemann solver is proposed. Different from previous works treating the elastic and plastic phases integrally, we resolve the elastic wave and plastic wave separately to understand the complicate nonlinear waves within the solid and then assemble them together to construct the exact Riemann solver for the elastic-perfectly plastic solid. After that, the exact solid Riemann solver is associated with the fluid Riemann solver to decouple the fluid-solid multi-material interaction. Numerical tests, including gas-solid, water-solid high-speed impact problems are simulated by utilizing the modified ghost fluid method (MGFM).

Semi-Lagrangian (S-L) methods have no CFL stability constraint, and are more stable than the Eulerian methods. In the literature, the S-L method for the level-set re-initialization equation was complicated, which may be unnecessary. Since the re-initialization procedure is auxiliary, we propose to use the first-order S-L scheme coupled with a projection technique to improve the accuracy at the grid points just adjacent to the interface. Standard second-order S-L method is used for evolving the level-set convection equation. The implementation is simple, including on the block-structured adaptive mesh. The efficiency of the S-L method is demonstrated by extensive numerical examples including passive convection of interfaces with corners/kinks/large deformation under given velocity fields, a geometrical flow with topological changes, simulations of bubble/ droplet dynamics in incompressible two-phase flows. In terms of accuracy it is comparable to the other existing methods.

This paper analyzes semi-discrete and fully discrete hybrid stress quadrilateral finite element methods for 2-dimensional linear elastodynamic problems. The methods use a 4 node hybrid stress quadrilateral element in the space discretization. In the fully discrete scheme, an implicit second-order scheme is adopted in the time discretization. We derive optimal a priori error estimates for the two schemes and an unconditional stability result for the fully discrete scheme. Numerical experiments confirm the theoretical results.

This paper introduces two novel conformal structure-preserving algorithms for solving the coupled damped nonlinear Schrödinger (CDNLS) system, which are based on the conformal multi-symplectic Hamiltonian formulation and its conformal conservation laws. The proposed algorithms can preserve corresponding conformal multi-symplectic conservation law and conformal momentum conservation law in any local time-space region, respectively. Moreover, it is further shown that the algorithms admit the conformal charge conservation law, and exactly preserve the dissipation rate of charge under appropriate boundary conditions. Numerical experiments are presented to demonstrate the conformal properties and effectiveness of the proposed algorithms during long-time numerical simulations and validate the analysis.

In present work, we investigate numerical simulation of steady natural convection flow in the presence of weak magnetic Prandtl number and strong magnetic field by involving algebraic decay in mainstream velocity. Before passing to the numerical simulation, we formulate the set of boundary layer equations with the inclusion of the effects of algebraic decay velocity, aligned magnetic field and buoyant body force in the momentum equation. Later, finite difference method with primitive variable formulation is employed in the physical domain to compute the numerical solutions of the flow field. Graphical results for the velocity, temperature and transverse component of magnetic field as well as surface friction, rate of heat transfer and current density are presented and discussed. It is pertinent to mention that the simulation is performed for different values of algebraic decay parameter α, Prandtl number Pr, magnetic Prandtl number Pm and magnetic force parameter S.

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.

The vectorial form of the Wave Propagation Method (VWM), regarding the dispersion of harmonic plain (elasto-dynamic) waves within certain wave-guides, is developed for the vibration analysis of circular cylindrical shells. To obtain this goal, all plain waves are divided into positive-negative going wave vectors along with the shell axis. Based on the Flügge thin shell theory, the shell continuity as well as boundary conditions are well satisfied by introducing the propagation and reflection matrices. Furthermore, all elements of the reflection matrix are derived for certain classical supports. As an example, for demonstrating the feasibility of VWM in the shell vibration analysis, a circular cylindrical shell with two ended flexible support is adopted. The natural frequencies of the systemaswell asmode shapes are obtained using VWM. The aquired results are compared with those of the previous works and found in excellent agreement. It is also found that VWM could mathematically provide a reduced dimensional matrix (dominant matrix) to calculate the natural frequencies of the system. Accordingly, the proposed method can provide high computational efficiency and remarkable accuracy, simultaneously.

In this paper, a hybrid dynamic mesh generation method for multi-block structured grid is presented based on inverse distance weighting (IDW) interpolation and transfinite interpolation (TFI). The major advantage of the algorithm is that it maintains the effectiveness of TFI, while possessing the ability to deal with multi-block structured grid from the IDW method. In this approach, dynamic mesh generation is made in two steps. At first, all domain vertexes with known deformation are selected as sample points and IDW interpolation is applied to get the grid deformation on domain edges. Then, an arc-length-based TFI is employed to efficiently calculate the grid deformation on block faces and inside each block. The present approach can be well applied to both two-dimensional (2D) and three-dimensional (3D) problems. The proposed method has been well-validated by several test cases. Numerical results show that dynamic meshes with high quality can be generated in an accurate and efficient manner.