We consider inviscid, incompressible shear flows with variable density and variable cross section. For this problem, we derived a new estimate for the growth rate of an unstable mode and a parabolic instability region which intersects semiellipse instability region under some condition.

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.

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.

The Stokes axisymmetric flow of an incompressible viscous fluid past a micropolar fluid spheroid whose shape deviates slightly from that of a sphere is studied analytically. The boundary conditions used are the vanishing of the normal velocities, the continuity of the tangential velocities, continuity of shear stresses and spin-vorticity relation at the surface of the spheroid. The hydrodynamic drag force acting on the fluid spheroid is calculated. An exact solution of the problem is obtained to the first order in the small parameter characterizing the deformation. It is observed that due to increase spin parameter value, the drag coefficient decreases. Well known results are deduced and comparisons are made with classical viscous fluid and micropolar fluid.

In this paper, we present a superconvergence result for the bi-k degree time-space fully discontinuous finite element of first-order hyperbolic problems. Based on the element orthogonality analysis (EOA), we first obtain the optimal convergence order of discontinuous Galerkin finite element solution. Then we use orthogonality correction technique to prove a superconvergence result at right Radau points, which is higher one order than the optimal convergence rate. Finally, numerical results are presented to illustrate the theoretical analysis.

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.

In this work, vortex convection is simulated using a dynamic mesh adaptation procedure. In each adaptation period, the mesh is refined in the regions where the phenomena evolve and is coarsened in the regions where the phenomena deviate since the last adaptation. A simple indicator of mesh adaptation that accounts for the solution progression is defined. The generation of dynamic adaptive meshes is based on multilevel refinement/coarsening. The efficiency and accuracy of the present procedure are validated by simulating vortex convection in a uniform flow. Two unsteady compressible turbulent flows involving blade-vortex interactions are investigated to demonstrate further the applicability of the procedure. Computed results agree well with the published experimental data or numerical results.

Understanding blood flow in human body’s cerebral arterial system is of both fundamental and practical significance for prevention and treatment of vascular diseases. The mechanism and treatment for the growth of daughter aneurysm on its mother aneurysm are not yet fully understood. Themain purpose of the present paper is to elucidate the relationships between hemodynamics and the genesis, growth, subsequent rupture of the mother and daughter aneurysm on the cerebral vascular. The intensified stents with different porosities and structures are investigated to reduce the wall shear stress and pressure of mother and daughter aneurysm. The simulation is based on a lattice Boltzmann modeling of non-Newtonian blood flow. A novel stent structurewith “dense in front and sparse in rear” is proposed,which is verified to have good potential to reduce the wall shear stress of both mother and daughter aneurysm. The simulation is based on a lattice Boltzmann modeling of non-Newtonian blood flow. A novel stent structurewith “dense in front and sparse in rear” is proposed,which is verified to have good potential to reduce the wall shear stress of both mother and daughter aneurysm.

In this paper, we propose an iterative two-grid method for the edge finite element discretizations (a saddle-point system) of Perfectly Matched Layer(PML) equations to the Maxwell scattering problem in two dimensions. Firstly, we use a fine space to solve a discrete saddle-point system of H(grad) variational problems, denoted by auxiliary system 1. Secondly, we use a coarse space to solve the original saddle-point system. Then, we use a fine space again to solve a discrete H(curl)-elliptic variational problems, denoted by auxiliary system 2. Furthermore, we develop a regularization diagonal block preconditioner for auxiliary system 1 and use H-X preconditioner for auxiliary system 2. Hence we essentially transform the original problem in a fine space to a corresponding (but much smaller) problem on a coarse space, due to the fact that the above two preconditioners are efficient and stable. Compared with some existing iterative methods for solving saddle-point systems, such as PMinres, numerical experiments show the competitive performance of our iterative two-grid method.

In this paper, we propose two new explicit multi-symplectic splitting methods for the nonlinear Dirac (NLD) equation. Based on its multi-symplectic formulation, the NLD equation is split into one linear multi-symplectic system and one nonlinear infinite Hamiltonian system. Then multi-symplectic Fourier pseudospectral method and multi-symplectic Preissmann scheme are employed to discretize the linear subproblem, respectively. And the nonlinear subsystem is solved by a symplectic scheme. Finally, a composition method is applied to obtain the final schemes for the NLD equation. We find that the two proposed schemes preserve the total symplecticity and can be solved explicitly. Numerical experiments are presented to show the effectiveness of the proposed methods.

Fluid-structure-interaction problems are ubiquitous, complicated, and not yet well understood. In this paper we investigate the interaction of a leading rigid circular cylinder and a trailing compliant filament and analyze the dynamic responses of the filament in the wake of the cylinder. It is revealed that there exist two flapping states of the filament depending on the cylinder-filament separation distance and the relevant critical distance distinguishing the two states is associated with the Reynolds number and the filament length. It is also found that the drag coefficient of the cylinder is reduced but that of the filament may be increased or decreased depending on its length. Compared with a single filament in a uniform flow, the filament of the same mechanical properties flapping in the wake of the cylinder has a lower frequency and a greater amplitude.

Numerical simulations of unsteady gas flows are studied on the basis of Gas-Kinetic Unified Algorithm (GKUA) from rarefied transition to continuum flow regimes. Several typical examples are adopted. An unsteady flow solver is developed by solving the Boltzmann model equations, including the Shakhov model and the Rykov model etc. The Rykov kinetic equation involving the effect of rotational energy can be transformed into two kinetic governing equations with inelastic and elastic collisions by integrating the molecular velocity distribution function with the weight factor on the energy of rotational motion. Then, the reduced velocity distribution functions are devised to further simplify the governing equation for one- and two-dimensional flows. The simultaneous equations are numerically solved by the discrete velocity ordinate (DVO) method in velocity space and the finite-difference schemes in physical space. The time-explicit operator-splitting scheme is constructed, and numerical stability conditions to ascertain the time step are discussed. As the application of the newly developed GKUA, several unsteady varying processes of one- and two-dimensional flows with different Knudsen number are simulated, and the unsteady transport phenomena and rarefied effects are revealed and analyzed. It is validated that the GKUA solver is competent for simulations of unsteady gas dynamics covering various flow regimes.

In this paper, a new combined method is presented to simulate saltwater intrusion problem. A splitting positive definite mixed element method is used to solve the water head equation, and a symmetric discontinuous Galerkin (DG) finite element method is used to solve the concentration equation. The introduction of these two numerical methods not only makes the coefficient matrixes symmetric positive definite, but also does well with the discontinuous problem. The convergence of this method is considered and the optimal L2-norm error estimate is also derived.

In this paper, we propose a new energy-preserving scheme and a new momentum-preserving scheme for the modified regularized long wave equation. The proposed schemes are designed by using the discrete variational derivative method and the finite volume element method. For comparison, we also propose a finite volume element scheme. The conservation properties of the proposed schemes are analyzed and we find that the energy-preserving scheme can precisely conserve the discrete total mass and total energy, the momentum-preserving scheme can precisely conserve the discrete total mass and total momentum, while the finite volume element scheme merely conserve the discrete total mass. We also analyze their linear stability property using the Von Neumann theory and find that the proposed schemes are unconditionally linear stable. Finally, we present some numerical examples to illustrate the effectiveness of the proposed schemes.

In this work, we introduce a mathematical model for the theory of generalized thermoelasticity with fractional heat conduction equation. The presented model will be applied to an infinitely long hollow cylinder whose inner surface is traction free and subjected to a thermal and mechanical shocks, while the external surface is traction free and subjected to a constant heat flux. Some theories of thermoelasticity can extracted as limited cases from our model. Laplace transform methods are utilized to solve the problem and the inverse of the Laplace transform is done numerically using the Fourier expansion techniques. The results for the temperature, the thermal stresses and the displacement components are illustrated graphically for various values of fractional order parameter. Moreover, some particular cases of interest have also been discussed.

This paper gives analysis of a semi-discrete scheme using equal order interpolation to solve unsteady Navier-Stokes equations. A unified pressure stabilized term is added to our scheme. We proved the uniform error estimates with respect to the Reynolds number, provided the exact solution is smooth.

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 apply an a posteriori error control theory that we develop in a very recent paper to three families of the discontinuous Galerkin methods for the Reissner-Mindlin plate problem. We derive robust a posteriori error estimators for them and prove their reliability and efficiency.

In this paper, we apply the collocation methods to a class of Volterra integral functional equations with multiple proportional delays (VIFEMPDs). We shall present the existence, uniqueness and regularity properties of analytic solutions for this type of equations, and then analyze the convergence orders of the collocation solutions and give corresponding error estimates. The numerical results verify our theoretical analysis.

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.