This paper presents a least-squares finite element method for linear Phan-Thien–Tanner(PTT) viscoelastic fluid flows. We consider stabilized weights in the LS method for the viscoelastic model and prove that the LS approximation converges to the linearized solutions of the linear PTT model; the convergence is at the optimal rate for the velocity in the H1-norm and at suboptimal rates for the stress and pressure in the L2-norm, respectively. For numerical experiments, we first consider the flow through a planar channel to illustrate our theoretical results. The LS method is then applied to a flow through the slot channel with two depth ratios and the effects of physical parameters are discussed. Numerical solutions of the channel problem indicate that flow characteristics of the viscoelastic polymer solution are described by the results obtained using the method. Furthermore, we present the hole pressure for various Weissenberg numbers, and compare with that derived from the Higashitani–Pritchard (HP) theory.

We investigate the singularity formation of a 3D model that was recently proposed by Hou and Lei (2009) in [15] for axisymmetric 3D incompressible Navier–Stokes equations with swirl. The main difference between the 3D model of Hou and Lei and the reformulated 3D Navier–Stokes equations is that the convection term is neglected in the 3D model. This model shares many properties of the 3D incompressible Navier– Stokes equations. One of the main results of this paper is that we prove rigorously the finite time singularity formation of the 3D inviscid model for a class of initial boundary value problems with smooth initial data of finite energy. We also prove the global regularity of the 3D inviscid model for a class of small smooth initial data.

No abstract uploaded!

We study a fourth order finite difference method for the unsteady incompressible Navier-Stokes equations in vorticity formulation. The scheme is essentially compact and can be implemented very efficiently. Either Briley’s formula, or a new higher order formula, which will be derived in this paper, can be chosen as the vorticity boundary condition. By formal Taylor expansion, the new formula for the vorticity on the boundary gives 4th order accuracy; while Briley’s formula provides only 3rd order accuracy. However, the use of either formula results in a stable method and achieves full 4th order accuracy. The convergence analysis of the scheme with our new formula will be given in this paper, while that with Briley’s formula has been established in earlier literature. The consistency analysis is easier than that of Briley’s formula, no Strang type analysis is needed. In the stability analysis part, we adopt the technique of controlling some local terms by the diffusion term via discrete elliptic regularity. Physical no-slip boundary conditions are used throughout.

In this paper, we study the superconvergence behavior of discontinuous Galerkin methods using upwind numerical fluxes for one-dimensional linear hyperbolic equations with degenerate variable coefficients. The study establishes superconvergence results for the flux function approximation as well as for the DG solution itself. To be more precise, we first prove that the DG flux function is superconvergent towards a particular flux function of the exact solution, with an order of $O(h^{k+2})$, when piecewise polynomials of degree $k$ are used. We then prove that the highest superconvergence rate of the DG solution itself is $O(h^{k+\frac 32})$ as the variable coefficient degenerates or achieves the value zero in the domain. As byproducts, we obtain superconvergence properties for the DG solution and the DG flux function at special points and for cell averages. All theoretical findings are confirmed by numerical experiments.