In the present work concentrated on the two-dimensional problem of generalized thermoelasticity for a fiber-reinforced anisotropic thick plate under initial stress. Using generalized thermoelasticity theory with fractional order heat conduction, the problem has been solved by a normal mode analysis. The effect of hydrostatic initial stresses and fractional order parameter is shown graphically on the distributions of the temperature, displacement and thermal stress components. It is found from the graphs that the initial stress and the fractional parameter significantly influences the varieties of field amounts.

In this paper, we propose a finite element method for the elasticity problems which have displacement discontinuity along the material interface using uniform grids. We modify the immersed finite element method introduced recently for the computation of interface problems having homogeneous jumps [20, 22]. Since the interface is allowed to cut through the element, we modify the standard Crouzeix-Raviart basis functions so that along the interface, the normal stress is continuous and the jump of the displacement vector is proportional to the normal stress. We construct the broken piecewise linear basis functions which are uniquely determined by these conditions. The unknowns are only associated with the edges of element, except the intersection points. Thus our scheme has fewer degrees of freedom than most of the XFEM type of methods in the existing literature [1,8,13]. Finally, we present numerical results which show optimal orders of convergence rates.

In this paper we extend the idea of interpolated coefficients for a semilinear problem to the quadratic triangular finite volume element method. At first we introduce quadratic triangular finite volume element method with interpolated coefficients for a boundary value problem of semilinear elliptic equation. Next we derive convergence estimate in H 1-norm, L 2-norm and L ∞-norm, respectively. Finally an example is given to illustrate the effectiveness of the proposed method.

In this work, element free Galerkin (EFG) method is posed for solving nonlinear, reaction-diffusion systems which are often employed in mathematical modeling in developmental biology. A predicator-corrector scheme is applied, to avoid directly solving of coupled nonlinear systems. The EFG method employs the moving least squares (MLS) approximation to construct shape functions. This method uses only a set of nodal points and a geometrical description of the body to discretize the governing equation. No mesh in the classical sense is needed. However a background mesh is used for integration purpose. Numerical solutions for two cases of interest, the Schnakenberg model and the Gierer-Meinhardt model, in various regions is presented to demonstrate the effects of various domain geometries on the resulting biological patterns.

The nonlinear Dirac equation is an important model in quantum physics with a set of conservation laws and a multi-symplectic formulation. In this paper, we propose energy-preserving and multi-symplectic wavelet algorithms for this model. Meanwhile, we evidently improve the efficiency of these algorithms in computations via splitting technique and explicit strategy. Numerical experiments are conducted during long-term simulations to show the excellent performances of the proposed algorithms and verify our theoretical analysis.

The nonlinear transverse vibrations of ordered and disordered two-dimensional (2D) two-span composite laminated plates are studied. Based on the von Karman's large deformation theory, the equations of motion of each-span composite laminated plate are formulated using Hamilton's principle, and the partial differential equations are discretized into nonlinear ordinary ones through the Galerkin's method. The primary resonance and 1/3 sub-harmonic resonance are investigated by using the method of multiple scales. The amplitude-frequency relations of the steady-state responses and their stability analyses in each kind of resonance are carried out. The effects of the disorder ratio and ply angle on the two different resonances are analyzed. From the numerical results, it can be concluded that disorder in the length of the two-span 2D composite laminated plate will cause the nonlinear vibration localization phenomenon, and with the increase of the disorder ratio, the vibration localization phenomenon will become more obvious. Moreover, the amplitude-frequency curves for both primary resonance and 1/3 sub-harmonic resonance obtained by the present analytical method are compared with those by the numerical integration, and satisfactory precision can be obtained for engineering applications and the results certify the correctness of the present approximately analytical solutions.

In this paper, we use moving mesh finite element method based upon 4P 1–P 1 element to solve the time-dependent Navier-Stokes equations in 2D. Two-layer nested meshes are used including velocity mesh and pressure mesh, and velocity mesh can be obtained by globally refining pressure mesh. We use hierarchy geometry tree to store the nested meshes. This data structure make convienence for adaptive mesh method and the construction of multigrid preconditioning. Several numerical problems are used to show the effect of moving mesh.

In this paper, a Chebyshev-collocation spectral method is developed for Volterra integral equations (VIEs) of second kind with weakly singular kernel. We first change the equation into an equivalent VIE so that the solution of the new equation possesses better regularity. The integral term in the resulting VIE is approximated by Gauss quadrature formulas using the Chebyshev collocation points. The convergence analysis of this method is based on the Lebesgue constant for the Lagrange interpolation polynomials, approximation theory for orthogonal polynomials, and the operator theory. The spectral rate of convergence for the proposed method is established in the L∞-norm and weighted L2-norm. Numerical results are presented to demonstrate the effectiveness of the proposed method.

Similarity solution is investigated for the synchronous grouting of shield tunnel under the vertical non-axisymmetric displacement boundary condition in the paper. The synchronous grouting process of shield tunnel was simplified as the cylindrical expansion problem, which was based on the mechanism between the slurry and stratum of the synchronous grouting. The stress harmonic function on the horizontal and vertical ground surfaces is improved. Based on the virtual image technique, stress function solutions and Boussinesq's solution, elastic solution under the vertical non-axisymmetric displacement boundary condition on the vertical surface was proposed for synchronous grouting problems of shield tunnel. In addition, the maximum grouting pressure was also obtained to control the vertical displacement of horizontal ground surface. The validity of the proposed approach was proved by the numerical method. It can be known from the parameter analysis that larger vertical displacement of the horizontal ground surface was induced by smaller tunnel depth, smaller tunnel excavation radius, shorter limb distance, larger expansion pressure and smaller elastic modulus of soils.

Mass migration of cells (via wave motion) plays an important role in many biological processes, particularly chemotaxis. We study the existence of travelling wave solutions for a chemotaxis model on a microscopic scale. The interaction between nutrients and chemoattractants are considered. Unlike previous approaches, we allow for diffusion of substrates, degradation of chemoattractants and cell growth (constant and linear growth rate). We apply asymptotic methods to investigate the behaviour of the solutions when cells are highly sensitive to extracellular signalling. Explicit solutions are demonstrated, and their biological implications are presented. The results presented here extend and generalize known results.

In this paper the space variable-order fractional Schrödinger equation (VOFSE) is studied numerically, where the variable-order fractional derivative is described here in the sense of the quantum Riesz-Feller definition. The proposed numerical method is the weighted average non-standard finite difference method (WANSFDM). Special attention is given to study the stability analysis and the convergence of the proposed method. Finally, two numerical examples are provided to show that this method is reliable and computationally efficient.

The beginning of the transition from the laminar to a turbulent flow is usually the generation of instability Tollmien-Schlichting (T-S) waves in the boundary layer. Previously, most numerical and experimental researches focused on generating instability T-S waves through the external disturbances such as acoustic waves and vortical disturbances interacting with wall roughness or at the leading-edge of flatplate, whereas only a few paid attention to the excitation of the T-S waves directly by free-stream turbulence (FST). In this study, the generating mechanism of the temporal mode T-S waves under free-stream turbulence is investigated by using direct numerical simulation (DNS) and fast Fourier transform. Wave packets superposed by a group of stability, neutral and instability T-S waves are discovered in the boundary layer. In addition, the relation between the amplitude of the imposed free-stream turbulence and the amplitude of the excited T-S wave is also obtained.

In the paper, we present an efficient two-grid method for the approximation of two-dimensional nonlinear reaction-diffusion equations using a expanded mixed finite-element method. We transfer the nonlinear reaction diffusion equation into first order nonlinear equations. The solution of the nonlinear system on the fine space is reduced to the solutions of two small (one linear and one non-linear) systems on the coarse space and a linear system on the fine space. Moreover, we obtain the error estimation for the two-grid algorithm. It is showed that coarse space can be extremely coarse and achieve asymptotically optimal approximation as long as the mesh sizes satisfy . An numerical example is also given to illustrate the effectiveness of the algorithm.

The finite difference (FD) method is popular in the computational fluid dynamics and widely used in various flow simulations. Most of the FD schemes are developed on the uniform Cartesian grids; however, the use of nonuniform or curvilinear grids is inevitable for adapting to the complex configurations and the coordinate transformation is usually adopted. Therefore the question that whether the characteristics of the numerical schemes evaluated on the uniform grids can be preserved on the nonuniform grids arises, which is seldom discussed. Based on the one-dimensional wave equation, this paper systematically studies the characteristics of the high-order FD schemes on nonuniform grids, including the order of accuracy, resolution characteristics and the numerical stability. Especially, the Fourier analysis involving the metrics is presented for the first time and the relation between the resolution of numerical schemes and the stretching ratio of grids is discussed. Analysis shows that for smooth varying grids, these characteristics can be generally preserved after the coordinate transformation. Numerical tests also validate our conclusions.

Direct numerical simulations of the transition process from steady laminar to chaotic flow are considered in this study with the relatively new incompressible lattice Boltzmann equation. Numerically, a multiple relaxation time fully incompressible lattice Boltzmann equation is implemented in a 2D driven cavity. Spatial discretization is 2nd-order accurate, and the Kolmogorov length scale estimation based on Reynolds number (Re) dictates grid resolution. Initial simulations show the method to be accurate for steady laminar flows, while higher Re simulations reveal periodic flow behavior consistent with an initial Hopf bifurcation at Re 7,988. Non-repeating flow behavior is observed in the phase space trajectories above Re 13,063, and is evidence of the transition to a chaotic flow regime. Finally, flows at Reynolds numbers above the chaotic transition point are simulated and found with statistical properties in good agreement with literature.

This article deals with numerical inversion for the initial distribution in the multi-term time-fractional diffusion equation using final observations. The inversion problem is of instability, but it is uniquely solvable based on the solution's expression for the forward problem and estimation to the multivariate Mittag-Leffler function. From view point of optimality, solving the inversion problem is transformed to minimizing a cost functional, and existence of a minimum is proved by the weakly lower semi-continuity of the functional. Furthermore, the homotopy regularization algorithm is introduced based on the minimization problem to perform numerical inversions, and the inversion solutions with noisy data give good approximations to the exact initial distribution demonstrating the efficiency of the inversion algorithm.

Multi-field coupling simulation method based on the physical principles is used to simulate the discharge characteristics of nanosecond pulsed plasma synthetic jet actuator. Considering the effect of the energy transferring for air, the flow characteristics of nanosecond pulsed plasma synthetic jet actuator are simulated. The elastic heating sources and ion joule heating sources are the two main sources of energy. Through the collisions, the energy of ions is transferred to the neutral gas quickly. The flow characteristics of a series of blast waves and the synthetic jet which erupt from the plasma synthetic jet (PSJ) actuator are simulated. The blast wave not only promotes outward, but also accelerates the gas mixing the inhaled gas from the outside cavity with the residual gas inside the cavity. The performances of PSJ actuator fluctuate in the first three incentive cycles and become stable after that.

Effects of viscous dissipation and radiation on MHD natural convection in oblique porous cavity with constant heat flux is studied numerically in the present article. The right inclined wall is maintained at a constant cold temperature Tc and the left inclined wall has a constant heat flux q with length S, while the remainder of the left wall is adiabatic. The horizontal walls are assumed to be adiabatic. The governing equations are obtained by applying the Darcy model and Boussinesq approximations. COMSOL's finite element method is used to solve the non-dimensional governing equations together with specified boundary conditions. The governing parameters of this study are Rayleigh number (Ra=10,100,200,250,500 and 1000), Hartmann number (0≤ Ha ≤20), inclination angle of the magnetic field (0° ≤ ω ≤ π/2), Radiation (0≤R≤15), the heater flux length (0.1≤H≤1) and inclination angle of the sloping wall (–π/3≤ϕ≤π/3). The results are considered for various values of the governing parameters in terms of streamlines, isotherms and averageNusselt number. It is found that the intensity of the streamlines and the isotherm patterns decrease with an increment in Hartmann number. The overall heat transfer is significantly increased with the increment of the viscous dissipation and the radiation parameters.

Jet impingement induced heat transfer is an important issue in engineering science. This paper presents results of large eddy simulation (LES) of normally impinging elliptic air-jet heat transfer at a Reynolds number of 4400, with orifice-to-plate distance fixed to be 5 in the unit of jet nozzle effective diameter . The elliptic aspect ratio (a/b) is 3/2. While the target wall is heated under some condition of constant heat flux. The LES are carried out using dynamic subgrid model and Open-FOAM. The distributions ofmean velocity components, velocity fluctuations, and subgrid stresses in vertical and radial directions, and the Nusselt numbers involving heat transfer through the target wall are discussed. The comparison with existing experimental and numerical results shows good agreement.

A high-efficient algorithm to solve Crank-Nicolson scheme for variable coefficient parabolic problems is studied in this paper, which consists of the Function Time-Extrapolation Algorithm (FTEA) and Matrix Time-Extrapolation Algorithm (MTEA). First, FTEA takes a linear combination of previous l level solutions as good initial value of Un(see Time-extrapolation algorithm (TEA) for linear parabolic problems, J. Comput. Math., 32(2) (2014), pp. 183–194), so that Conjugate Gradient (CG)-iteration counts decrease to 1/3~1/4 of direct CG. Second, MTEA uses a linear combination of exact matrix values in level L, L+s, L+2s to predict matrix values in the following s–1 levels, and the coefficients of the linear combination is deduced by the quadric interpolation formula, then fully recalculate the matrix values at time level L+3s, and continue like this iteratively. Therefore, the number of computing the full matrix decreases by a factor 1/s. Last, the MTEA is analyzed in detail and the effectiveness of new method is verified by numerical experiments.